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Abstract 

An attempt has been made in this review paper to study the effect of couple-stress on some 
problems of classical elasto-dynamics. Using the linearized asymmetric theory of couple-stress, 
presented by Mindlin and Tiersten, the authors have reported the following investigations on 
elastic waves and distribution of stresses in solid elastic material. In Cosserat's pseudo- 
continuum, the propagation of surface waves is considered paying special attention to Rayleigh, 

Love and Stonele type of waves. The propagation of waves in a solid elastic layer has been 
studied thoroughly in the arena of couple-stress theory of elasticity. Axis symmetric Lamb’s 
problem is also studied in this review paper Remembering the engineering significance the 
effect of couple-stress on steady-state response to moving loads in a semi-infinite elastic 
medium is studied in details considering the three possible cases - supersonic, subsonic and 
transonic, depending on Mach numbers. Wave velocity equations are deduced and in particular 
in the last problem integral transform technique has been adopted and finally the expressions 
for displacements and stresses in improper integrals are evaluated and written in terms of real 
functions. The results obtained during the investigation in this review paper are in fair 
agreement with the results of the corresponding classical problem when the couple-stress 
parameter tends to zero. 

(Keywords : couple-stress/Cosserat's pseudo-continuum/ asymmetric theory of elasticity/ 
Galilean-transformation/Mach numbers). 
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Introduction 

The classical theory of elasticity is based on an ideal model of an elastic 
continuous medium in which the loading is transmitted through an infinitesimal area 
element in the body by means of stress vector only. The deformations thus take place 
is characterized by the symmetric strain tensors. The symmetric stress and strain 
tensors are interrelated by the generalized Hooke's law. Though the results obtained 
for homogeneous, isotropic materials with the application of the classical theory of 
elasticity are in harmony with experiments, in many cases, remarkable discrepancies 
between theory and experiments have been observed. The discrepancy between the 
classical theory of elasticity and experiment is striking where the stress concentration 
takes place in solid elastic medium i.e. in the neighbourhood of holes, notches, 
grooves, cracks and also in case of ultrasonic waves. The classical theory of elasticity 
eventually fails in the study of vibration of grain bodies and polymer'. 

It is expected from the mechanical point of view that the forces across a 
hypothetical plane within the solid should be statistically equivalent to a force and a 
couple. Such an assumption relays the fact that not only force stresses but also the 
couple stresses are transmitted through an area element, both of them obviously are 
asymmetric in nature. The corresponding deformations are also characterized by two 
asymmetric tensors namely strain tensor and curvature twist tensor. These are the 
main theme of couple-stress theory of elasticity. 

The concept of couple-stress was originally introduced by Viogt^. The complete 
theory of asymmetric elasticity was developed in 1909 by the Cosserat brothers^. In 
this theory, which wasnotvlinear in the beginning, they assumed to each molecule a 
perfectly rigid trihedron which during the process of deformation underwent not only 
the displacement but also the rotation. The material medium is termed as Cosserat 
continuum or micropolar elastic medium. In spite of the novelty of the idea, the work 
of Cosserat brothers was not duly appreciated during their life time and unnoticed for 
a good while. 

After a long time, the interests of research workers were concentrated on the 
simplified Cosserat theory of asymmetric elasticity of the so-called Cosserat Pseudo- 
continuum. By this same we understand a continuum for which the symmetric stress 
(force stress and couple-stress) may occur while the displacement of a body is 
described by a single displacement vector only. The modern derivation of the 
Cosserat's ^theory has been given by Truesdell and Toupin^ Toupin^ Aero and 
Kuvshinki Grioli*, Mindlin’, Mindlin and Triesten'*’ and Cohen". Following the 
linearized form of constitutive equations in couple-stress theory studied by Mindlin 



PROBLEMS IN COUPLE-STRESS THEORY ELASTICITY 


3 


and Tiersten some problem of elasto-dynamics have been studied in this review 
paper. 

The investigation of propagation of waves (surface waves and waves in a layer) in 
solid, elastic, homogenous and isotropic medium is one of the main feature of 
classical elasto-dynamics. Both the problems have been studied by a good number of 
investigators'^'^''. Regarding Lamb problem, it is noted that Sengupta and Ghosh^^ 
studied the axis symmetric Lamb problem in couple-stress theory of elasticity and 
Nowacki and Nowacki^^’^’ studied the plane and axis symmetric Lamb problem in 
micropolar theory of elasticity. In classical theory of elasticity, the problem of moving 
load over a solid elastic semi-space has been studied by Cole and Huth^*, Lamb^® and 
Sneddon^®. The designs of highways or airport runways, as well as the foundation 
problems in soil mechanics, particularly when the earth-mass supports a moving load 
over its free plane surface lead to the investigations of the dynamic stress distribution 
associated with the problems. It is also noted that Nath and Sengupta^' investigated 
steady-state response to moving loads in solid elastic media considering the 
supersonic-supersonic case. 

Sengupta and this research collaborators have studied different problems of elastic 
waves and vibrations in couple-stress theory of elasticity^^''"’. In this connection it is 
mentioned that another part of asymmetric elasticity, namely the mircopolar theory of 
elasticity has been developed by many authors^'^®. 

The present review paper is concerned with the general linear theory of Cosserat 
medium. We confine ourselves to the problem of elastic, homogeneous, isotropic and 
centro-symmetric bodies. Here like the classical, in the linear asymmetric theory of 
elasticity we assume that the deformations are small and square and product of the 
deformation are negligible with the respect to the linear terms. It is also assumed that 
the relation between the state of strain and stress are linear and the increase of 
temperature is inconsiderable. Lastly, we shall consider the theory of couple stress as 
macroscopic theory like the classical theory of elasticity and distance considered 
within its frame are much greater than the intermolecular distance. Hence, it may be 
assumed that the action radius of the intermolecular forces is negligible and equals to 
zero. 

Basic Equations and Relations : 

We shall consider the following basic equations and relations in studying different 
problems in this review paper. The states of stresses and couple-stresses are given in 
terms of the Cartesian tensor notation'". 



4 


P. R. SENGUPTA et al. 


x] = lixey + Xekj^ij ; = 4r]Xij + Axixjt (1) 

where the local strain Cy, rotation co/ and curvature twist tensor Xy are expressed in 
terms of the displacement m, as 

«y = {Wy + My)/2}; 00/= -{Sijt ujk)/2}; Xy = %/ (2) 

where X, p are Lame's constants, 8y is Kronecker delta, t] and t] are the material 
constraints associated with resistance to curvature, Sy* is the unit alternating tensor and 
comma denotes the partial differential with respect to the space co-ordinates. 

The deviation of the couple stress tensor p° and symmetric stress t* are given by 
P” = Py-{(p^y)/3}, t: = {(Ty + ty)/2} (3) 

The stress equation of motion, couple stress equation of motion and linearized stress 
equation of motion are given by the vector form'° as 

V.i + p7 = p#; + + (4) 

ox oy dz 

V.p + pc + Tx/ = 0 


V.t‘ + (l / 2)V X V.p" + p 7 + (l / 2)p V X c = p# (5) 

where / (sVr) is the unit spatial dyadic, f denotes the body force vector, c is the 

body couple vector and dot denotes the time derivative Inserting (1) in (5) We get the 
displacement equation of motioin'® as 

pV^u + (A, + p)W.M + qV^Vx VxM + p7 + (l/2)pVxc = i< (6) 

The equations of motion given by (4) and (6) may be written in terms of tensor 
notations as 


(L/=1, 2,3) 




^ 

dx^ dy^ dz^ 


\j + Pfi = PUi 

jij POj — 0 


( 7 ) 
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\iV\ + (X + Ji) « (u^^, - V^u)+ p/ - H , = pu- (8) 

We shall deal mainly with the following displacement equations of motion in absence 
of body forces and body couples [putting/ = c/ = 0 (i =1, 2, 3) in equiation (8)] in 
studying this review paper. 

pV'M, + (X + p)m^^, - V^M,) = pM, 2, 3) (9) 

Effect of Couple Stress on Surface Waves : 

Let us consider two homogeneous elastic medium Mi and M2 welded in contact 
(or sufficiently rough enough to prevent any sliding on common surface) at their 
common surface of separation. The two medium are separated by a plane horizontal 
boundary extended to infinity and M2 being above Mi. As a reference coordinate 
system, a set of orthogonal Cartesian axes OX1X2X3, the origin O being any point of the 
boundary and OX3 pointing normally into M2 is taken. 

Let us consider the possibility of a wave traveling in the direction Oxi in such a 
manner that (a) the disturbance is largely confined to the neighbourhood of the 
boundary, and (b) at any instance all particles in any line parallel to 0x2 have equal 
displacements. Due to (a), the wave is a surface wave; and due to (b), the case we 
have taken is analogous to the plane waves. Then the displacement components, Wi 
and W3, at any point may be expressed in the form 


Ml = O i + y 3, M3 = 4>,3 - M/.i, 


( 10 ) 


so that V^<D =0, VV = mi 3 -M3 1 s 

dx^ dxl 

where O and \|/ are functions of coordinates xi, X3 and time t. 

The displacement equations of motion are from (9) 

pV^M, +(X + p)m,,, +T1* V^(m,,, -V\) = pm, (i,J= 1, 2, 3) 


( 11 ) 


Then, in presence of couple-stresses we obtain, in view of (10) and (11), the 
following relations in M \ : 
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v^(i)-o/<=o, = (12) 

where pi, A-i, pi denote the properties of the media M\ and 

aj = {(X,, + 2p,)/p,}, p^=pi/p, and =Ti;/p, , ( 13 ) 

p' being a constant characterizing the effect of couple-stress in the medium M\. 
Similar relations hold good in M 2 with pi, A,,, pi f, , ai P, replaced by p 2 , X 2 , P 2 f j, 
(X 2 P2- 

To solve the equations (12) and similar equations for the medium we assume 

O = / 'P = (14) 

On substituting equation (14) in (12) and similar equations for the medium M 2 , and on 
simplification we get the solutions satisfying conditions for surface waves, as follows; 

In the medium Mi 

T = +Ce^>^^ 

M2 - e'iM. 

In the medium M 2 

0 = £)e'^'2''3+^-«) 

'P= 

where B, C, D, E,F,B , C', F* are all arbitrary constants and vj, r\j(J=\,2) are all 
positive imaginaries and ^ are positive real numbers given by 
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0 = 


Pj 




■sl /2 


Xj +2]Xj 


-1 


={p) 


1/2 


y 


2 tj 


^ p . 

V y 


-1 


^ q] 


2 i) 


, y 


AV 2 


+1 


( 17 ) 


The above system of equations ( 15 ) and ( 16 ) will lead us to a particular solution 
corresponding to a group of simple harmonic waves of wave-length 2 %l\ traveling 
forward with speed c. 

To obtain the velocity equations we have now to apply the following boundary 
conditions : 

(i) The displacement at the common boundary surface between M\ and M2 must 
be continuous at all times and places. 


(ii) The rotation components eoi, 02 and 03 where ©1 =--^, ©2 = ©3 = 


*2,1 


must be continuous on the common boimdary. 


(iii) The stresses a3i, 032, 032 and couple-stress p3i, 1J.32, M.33, where 


=P.“2,3-'n:k,„+W..3.j 

0,3 =X,V^(D + 2p. ((5,33 -TJ 

P 31 = -2r|>.33, P 32 = 2ti*VV.3. 


and similar expressions for M2, across the boundary surface between M\ and M2 must 
be continuous at all times and places. Using ( 10 ), ( 15 ), ( 16 ) and the condition (i), we 
obtain 


^ - KxC^ + £ + Ti2£ + K2F, 

-nA-{B + C) = r2D-{E + F), 


( 19 ) 

( 20 ) 
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b* + c*=e*+f* 

Again, substituting (15) and (16) into (19) and using condition (iii), we obtain 

-tlOh -nrfe 

= -2li,r,eD4 e -rhk fe' -rh^l-Vj f; 


(22) 

knM +^')+ wk +/teife' -C?)+C,P:)c* = 

=-knM +r)+n#j£* -/kC.Cr (23) 

+ 1 )+ 2ii,r:^^A + + 2ii, ^i;jC = 

= {x,(r/ + l)+2ii/:kD - 2 \iM.E - 2ix,^t;jF; (24) 

«ni(n? +^^K*^ + Cifc = -/ri2(ri^ + 'ek)(\iE (25) 

n* (- +CiC*)= ri;(iiij£* - C2^‘)> (26) 

and from the condition (ii), we obtain 

/il,5* + C,C’ = /n.£*-C,F-, (27) 

- (nf + + (c? - = -(n^ + f )f , (28) 

and other component <02 contributes the same equation as (21). 


From the equations (21), (23), (26) and (27) we find that only possible values of 
B*,C*, E* and are zeros. Thus there is no propagation of the displacements U 2 . 

Again, equation (25) will be satisfied if 
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K 


qXi ^ 


qXi 


(29) 


where X is a constant which depends on the parameters of the couple-stress of the 
media. From the equations (28) and (29), we obtain 


K = v^plB-vy^E, 


(30) 


wtere 


V, =(! + //<;,)/ 


M 1 ^ 
+ 


C.Ti; <;^r\ 


, v,={l + i/Qj 


M 1 ^ 


2 2 




(31) 


In view of (3 1) the equations (19), (20), (22) and (24) can be written as 




q,M--/vi;7fa 

V 


Z' V. „2 ^ 

^2-^ + lV2Pla 

I 


E = 0 


r^A + 


1 i^\P\ 

A 




+ Vipfb]B + r^D-\\- -5^ + v^plb U = 0 




2p,r,^^^ + 


-nlVz -Tl*2/22i-»fe^ +<;^K -'n24'2)-^^ + V2/2|<^ 

?292 


B-¥2p2^l% D 

{(x, + 2p, ) (r,^ + 1)- 2p, ) ^ + |2p, 

{(Xj + 2^2) (^ 2 ^ + 1)- 2^2}^ + 12^2 


£• = 0 


J 


\ ?i 

r 


+ ivipfd \ B- 


2^ 


^2^2 

tl2+-y 

92 ; 


-IV2j3|i/^£ = 0, 


( 32 ) 
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where 


11^1 1 
a = ^- r -^- r . b = - — r-r + - 


+ck -^ 0 '}.— TT + fe' +ck 

qX qX. 

(33) 

Eliminating A, B, D and F from the equations (32), we obtain the modified wave- 
velocity equation of the surface waves under the influence of couple-stresses, 






B,+~ 


f ^ 

-iv^plJ. 

q^j 

j .[ 

\ 92 j 

J 


1 -- 


2 


2p,P/,^' 


Ci^i 


'te’-c)-';9:}x2^?M;+v,p,-c 


y ^2 qi j 


!-■ 




te’+c)-<’,9;}K^+v.p.’c 


' f . X , 

f 'Kl ^ 

1 T 

P.iaK^,'+l)2p,P? 


+ ivXd-2rM 

L 1 ■ ' 

1 9 . ^ 

J JL 


■2Pfkl 


-p,{a^,(r/ +l)}2p^P 
-2P^k] 


2 

n 

+XpX 

V ql ) 

I 


(34) 


=0 


where a^, P) and t] are given by (13) andj = 1,2. 
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If the parameters of couples=stress f , and vanish, the wave-velocity equation 
(34) as obtained above is in agreement with the corresponding classical result, which 
may be written in the simplified form. 

c\(P\ - Pif + (pi^2 + P2O)] + 2 Qc^ [riTiiP2 - 2-2T12P1 + P2- Pi] + ^{rxr\\ + 1 ) 

(/•2Tl2+l) = 0, (35) 

where q = 2(p,p^ -p,p]) = 2(p, -pjand = { (c7a;)-lf , ri^ = { (c7p;)-lf 

Particular cases : 

Ralieigh waves : The particular case of the foregoing problem is to investigate the 
possibility of Rayleigh waves in an elastic solid under the influence of couple-stress. 
In this particular case the plain boundary must be a free surface so that is replaced 
by vacuum. In view of (22), (24) and (25) we obtain : 

2p,r,4^^+[{(^^ --Ti?)pi -^xqtWpx Kx(1\]b = ^, 

(36) 

{ (A., 2p.)(r,= +l)-2p.}^ 4-2 p,{ti, + 5 = 0. (37) 


Eliminating A and B from (36) and (37), we obtain ; 

4PiV,^{ni + (n,/?? / q\ ))- { {li + 2p, ) (r,^ + 1)- 2p, ) 

[{ fe ^ -■»!?) 1^1 ■^Q)px-^WxWpxKx(1x\ =0 

(38) 


Equation (38), in view of (17) becomes 


f 2 


f 2 > 

- 

f 2> 



2 A 

f 2 ^ 

^-1 

2 

^-il 

2 

1+4 

= 

4-2 

2-4 

1 

Ui J 






Ip? J 

P?J 

J 


Equation (39) determines the velocity of the Rayleigh surface wave under the 
influence of couple-stress. 

If the parameter f ^of couple-stress be zero, we obtain the classical result’® 
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{2-(cvi3;)r=4{i-(c7c;r{i-(cvp:r 

The equation (39) can be analysed to study the effect of couple-stress on the said 
classical problem by assuming the couple-stress parameter to be so small that its 
cubes and higher powers are always neglected. Thus the approximate form of equation 
(39) is 


2 



= 4 










f 

1 + 

V 


Pi 


2 

1 

J 


( 40 ) 


Squaring both the sides of (40) and re-arranging we obtain 

l-lbl 




c' c* 


^24 16 ^ 


p; 


1-^ 


V 


a 


Pf a 

-32^^4i- 


■ y 


5. + £? 
a? P: 


^2 ^ 

2 




P. 

\ 

1 



J a^P.\ 

1 J 


= 0 


("41) 


We assume X = p. (i.e. v = Vi), then ai = •\/3p, , and we obtain from equation (41) 


vp: p: 3 K 3 


16 




3-- 


pr 




= 0 


(42) 


The above equation is a cubic equation in cYPJ and has three real and positive 
roots. One of the roots lies between 0.8453 and 2. In order to satisfy the condition of 
Rayleigh waves the value of must be less than or equal to 0.09, and the 
acceptable roots for the real existence of Rayleigh waves have been computed 
numerically. It is thus observed that the velocity of Rayleigh waves increases due to 
the influence of couple-stress. 

Numerical result : 

In Table 1 the values of c/Pi have been computed for different values of in 
equatin (42), and thus it has been exhibited that the velocity of Rayleigh waves in 
presence of couple-stresses increases. 
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Table 1- Computed of c/Pi for different values of the couple-stress parameter 

10 '^ 10 '^ 10 ^ 10 -’ 10 -^ 
c/p, 0.92915 0.92288 0.92199 0.919411 0.9194025 

Love waves ; 

Let us suppose that the medium Mi is bounded by tvt'o horizontal plane surface at 
finite distance //apart, the upper plane being free, while the lower plane surface forms 
the medium Afi-It is sufficient to consider the displacement ui only. The notable fact 
here is that the displacement in M 2 may no longer dimnish with distance from the 
boundary between Mi and Mi, so that for the medium Mi we presence the full solution 
as 


)+ (43) 

and for the medium Mi, the solution is 

w, = + C*e"’ , (44) 

where ‘r]i is not necessarily imaginary, while for M , we see that r|i is an imaginary and 
^ 1,^2 are positive quantities. 

In addition to the three boundary conditions (i), (ii) and (iii) for general surface 
waves we have the conditions that there shall be no stress and couple-stress across the 
free surface X 3 = H. Hence, in view of (43), (44) the boundary conditions may be 


written as 

B* + C* = (Ai+ Ai) + (Bi + Bi), (45) 

iniB* + ; ,C* = iii2 (^1 - ^ 2 ) - 2 (Bi-Bi), (46) 

- /ri,p, (1 + ) 5 * +p,C, 

= /q,p,(l + i\p:){A - A,y p:C,(i - i]ql){B, - B,) (47) 

(-21^5- + CJC'W: = U ( 5 , + (48) 
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-5,e^'')=o 

-« 1 ^ ^A,e-‘'^-)+Cr -A ^'''')=0 


From (48) we, get 


C’=- 


K 


h^?Cf Cl 




K 


M'a'^aCi 


+ ^(4 +^2)= 

C 2 


where is an arbitrary constant which depends on the parameters. 
Equation (49) and (50) are satisfied if 




(49) 

(50) 


(51) 


(52) 


which gives 


= i tanti,i/ and Bi-B 2 = {B^ + Bt) tan h(;^H 

A + A. 


(53) 


In view of (46), (51) and the second equation of (53), we obtain 
K= (1/a) {bB* + c(Ai + A 2 )}, 


i: = (l/a)W+c(4+4)} 


where a = 


P2^2Cl 


M'1^iC2 J 


Cl 


712 




(54) 


In view of (51), (54), we obtain from (45) and (47) 


= + ApHA - A)]- ^ - ^^^Kn^/COtan KjH + cJ(A, + 4)] 


( 55 ) 



PROBLEMS IN COUPLE-STRESS THEORY OF ELASTICITY 
where 



( \ 1 ^ 

c 

( \ 1 ^ 

a 


a 
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<■! t%; m 


15 


(56) 


Equation (54), in view of the first equation of (52), becomes 


/q.(l + £]p^)+ (l - 

= + ^\P^)] tan^,/f - {(l - (il:/C') tan /?;,// + cj } 


(57) 


The above equation gives the S. H. surface waves of velocity under the inluece of 
couple-stress. If the parameters of couple-stresses, f jUnd vanish, the wave-velocity 
equation is in agreement with the classical equation of Love waves. Hence, equation 
(57) shows the remarkable effect of couple-stresses on the Love waves. 

The modified equation can be analysed to study the effect of couple-stresses on 
the classical problem by assuming the parameters of couple-stresses, ^,and to be 
small so that their cubes and higher powers are always neglected. Hence, on the above 
assumption, we obtain the following approximate form of equation (57): 


where n, =^{(cVp;)-(cVp;)^^f: -l}(/= 1,2) (59) 

Equation (58) yields a real value of c if rii are imaginary and real respectively. The 
requirement that qi, should be imaginary and qi be real is, by (59), 

P,(l + ^^ft)<c<p.(l + ^^fO 


(60) 
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which shows that for a particular ^ the range of c increase under the influence of 
couple-stresses. In other words, the velocity c of Love waves increases under the 
influence of couple-stresses. 

Stoneley waves : 

In the classical theory Stoneley waves are a generalized form of Rayleigh waves 
propagating along the common boundary of M\ and M^. Hence, Stoneley waves along 
the common boxmdary of the media M\ and M 2 under the influence of couple-stresses 
are determined by the roots of the wave equation (25). When the couple-stress 
parameters vanish, we get equation (26), which is the frequency equation of Stoneley 
waves in the classical theory. 

Effect of couple-stress on waves in a layer : 

Let us introduced a Cartesian frame of reference OX1X2X3 taking the origin in the 
middle plane of the elastic layer; the middle plane coincides with the plane 0x1X2. We 
consider the effect of couple-stresses on the propagation of waves in an elastic layer 
of thickness Ih. The planes bounding the layer are X3 = ± /j and are supposed to be free 
of stress. There exists a plane wave moving with the constant velocity c in the 
direction of Xi. Both the longitudinal and transverse waves in the infinitely extended 
layer would be propagated. The boundary surfeces of the elastic space leads evidently 
to a distortion of the state of stress, which influences the velocity of propagation of 
elastic waves. From the nature of the problem the non-zero displacement Ui and M3 at 
any point may be expressed in the form. 

Ml = <t>,l - Y,3, M 3 = (t»,3 + IK.! (61) 

where <j) and y are displacement potentials, functions of the coordinates Xi, X3 and 
time. 

The displacement equations of motion are 

+(X+p)m,^ +'n'V^(M,^ - v\)= pm; [ij= l, 2, 3] (62) 

The equation of motion (62) in view of (61), yield the following equations 

V>-c;‘$ = 0; VV-fVV = c,'i|> (63) 
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where c,* = (X + 2^)/p , < = l^/p, ^'='n/P = 

r| is a constant characterizing the existence of couple-stresses and X, p are Lame's 
elastic constants. If , the parameter of couple-stress, be zero the classical result of 
the corresponding problem follows at once. 

Oiur problem here is to seek solutions of the equations (63) subject to the 
boundary conditions 


033 = CT 31 = p .32 = 0 in the plane X 3 = ± h (64) 

The stresses a,j and couple-stresses Py are given by' 

O33 = 2 pe 33 +^(^11 +^33)5 

013+031 =4p«i3, 

®13 +5®31 =“M'13.1 +P 32.3 

On substituting for pi and pa from (61) in (65) we get 

011 = 2p((t>,ii - 033 = 2p((t) 33 - vj/ ^ J+ 

013 =p(2fi3 -V .33 +Y,ii)-nVV tJai =p(2<1>.3i -y .33 +¥.ii)+^VV 
P32 = - 2 tiV^(vp 3 ) p,2 = -2 tiV^(\Ki) 

( 66 ) 


To solve the equation (63), we assume 


<|>(^b :C 3 , 0 =fixi, X 3 )e®' ; 

\p(xi, X3, t) = g (xi, X3) 

(67) 

Owing to (67), equations (63) take the form 



vy+k:f=o ; 

V^g+k:g-rV^g = 0 

( 68 ) 

where k^=(p^/c^) , 




= 1,3 ^ 


011 — 2 peii + X(eii +£ 33 ) 
II 32 =4ri(e3i.3-e33_i) 

Pi2=4n(en.3-%.i) 


^/ = 1,3^ 
J = U 


(65) 
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The solution of the equation ( 67 ) are 


Introducing ( 69 ) into (68), we arrive at two ordinary differential equations 




r-{a^-k:)r=o 


Kdxlj 


-a‘ ^ - 


ydx] j 


-a^-^ng -=0 


( 70 ) 


the solution of which are 

/■* = A sinh v,:c3 + 5 cosh VjAIj 

g" =C sinh.V2jr3 +£)coshv2Jf3 +£'sinhv3j:3 +Fcoshv3;c3 

where 

v,=^a'-k^ , Vj v,=^a"-q\ 

={\/2e)[{\+4k:tf -l]; ={]/2e)[{\+4kief -l] 

By virtue of equations ( 71 ), ( 69 ) and ( 67 ), we finally obtain 

<t) = {^sinh^VjJCj + J 5 cosh v,X3)'^‘“'"‘^'^| 
v|; = (C sinh V2JC3 +i)cosh v2A:3 +E sin h,V3X3. + F cosh V3x:3) 

we consider the particular case 

4,, =Scosh:v,x, 

y V /(cor-Qtxii ) 

v|/, = (C sinhiV2Jc:3 + E sinh V3X3 + /^cosh V3 a:3) 




Introducing ( 73 ) into ( 61 ), it is readily observed that in this case the displacement are 
symmetric with respect to x.^ = 0, so are the stresses O33, 031 and couple-stress p.32. 
Besides, the stress Ou and the couple-stress pi2 are antisymmetric with respect to the 
plane x:3 = 0. 
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Substituting from (73) into the boundary conditions (64), we get, 

(pco^ - 2|j.a^)5cosh^V|/r + 2/|a,a (CVj coshiVj/r + Ev^ cos hv^h) = 0 
2rav,5sin hvj 2 + C {(v^ +a^)-f^{vj +a^y}sin hvjt 
+ E {(v; + a^)-t{vl + }sin hv^h = 0 
Cvj (vj - ) cos hvj 7 4- Ev^ (v, - ) cos hv^h = Q 

Since i is small, the system of equations (74) takes the form 

[a{(cYcj ) - 2}cos /zv,/r]5 + 2i{CVj cos hvji + Ev^ cos hvjj) = 0 
2iBav^ sinh v^h + C(2a^ - )sinh v^h + Ei2a^ - kl )sinh vjj = 0 '(c = a/ a) (75) 
- Cv^p^ cosh Vj/i + Ev^q^ cosh v,h = 0 



These equations lead to the transcendental equation 


tanhvjA 

tanhvj/i 


= {2 - (c V cl )\2a^ - kl j/ 4v, V 2 {l + {p^ j ql )} j^l + (v 



It is readily observed that the wave velocity equation (76), as derived above for 
the propagation of waves in an elastic layer, contains the couple-stress parameter i . 
Assuming the couple-stress parameter f to small so that is cubes and higher powers 
are always neglected. To discuss the result of the wave velocity equation (76) in the 
following cases. 

Case I : If the length of the wave is large compared with the thickness of the layer 
2 h, the quantities v\h, vjh and ah can be regarded as small and the hyperbolic tangents 
can be replaced by their arguments. But V 3 A is not small. So the equation (76) becomes 



Equation (77) determines the velocity c of the plane wave under the influece of 
couple-stress in an infinitely layer. 
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If the parameter ^ of couple-stress vanishes, we get from equation (77) 

■)} = <) 

i.e. = 4Cj(ci‘ -Cj) which is in agreement with the classical result of Rayleigh^’ 
and Lamb^*. 

Since the parameter ^of couple-stresses is small, we get the following 
approximate form of the equation (77) 

4v,^ {l + kie) = {2- I cl ) (2a^ - kl ) 

which can be rewritten as 
4(1 - c^lcl ) {l + (c ^ } = (2 - (78) 

8 r 1 

we assume X - p,;(v = 1/4); thenc,^ = 3cl, and c^/cl =- 1 + — 

3 V 6 

which shows that the symmetrical mode of vibration dcj increases due to the presence 
of couple-stress. On the other hand the velocity c of the propogation of waves in the 
elastic layer increases due to the presence of couple-stress. 

Case II ;If, again, the length of the wave is very small compared with the thickness of 
the layer 2h, the quantities v\h, v^h, v^h and ah are large and we may assume that the 
ratio of hyperbolic tangents in equation (76) approaches unity and equation (76) 
becomes 

4{\-dlclT =^{l + P^lg%2-dlciy{u{vyM)} (79) 

Equation (79) determines the velocity of the Rayleigh surface wave under the 
influence of couple-stress. If t = 0, couple-stress vanish and we get the classical 
result 

(2 - c V )^ = 4 {1 - (c V c? )P { 1 - (c V 4 r } 

The couple-stress parameter i being small we get the following approximate form of 
the equation (79), 
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('^V‘!‘)-8(c‘Ai)+‘:^f24/cf)-(l6/cf))-16jl-(cjVci'))-16{c7<;2)x 

x{l-(cV=?))-32avjl-f^ 

I V^i 


2 "*'T 

' C, 


C, Cj 


(80) 


We assume A, = |j. (i.e. v =1/4), then c, = -y/SCj and obtain from the above equation (80) 
56 32' 


'"c‘ c 


^2 3 3 


1 16 ,,, 




a i 

3 — r 


2- — 

J 3 

1 ^2 J 


^ ^2 y 


= 0 


(81) 


The above cubic equation in c^lc\ has three real and positive roots. One of the roots 

lies between 0.8453 and 2. In order to satisfy the condition of Rayleigh waves the 
value of must be less than 0.09 and the acceptable root for the real existence for 
Rayleigh wave been computed numerically. It is thus observed that the velocity of 
Rayleigh waves increases due to the influence of couple-stress. 

Numerical result : 

In the following table the values of dcj have been computed for different values 
of in equation (81) and thus it has been exhibited that the velocity of Rayleigh 
waves in presence of couple-stresses increases and the velocity of Rayleigh waves 
under the influence of couple-stress, approaches to the classical value 0.9194019 cj 
with the diminution of the value of the parameter of couple-stress. 


Table 2- Computation of dci for different values of couple-stress parameter. 



mumoni 

10’^ 

10 -^ 

10"^ 

10"^ 

dct 

0.92915 

0.92288 

0.92199 

0.919411 

0.9194025 


Let us consider another set of interesting solutions given by 


<j)2 = jfsinhviXjg'^®' 11/2 = (D cosh VjXj + jD cosh V3X3)e ^ (82) 


Introducing (82) into (61), it is readily observed that the displacement u\, 1/3, the 
stresses an, 033 and couple-stress jj.32 are antisymmetric with respect to the planers = 
0, while the stresses a^, 031 and the couple-stress jin are symmetric with respect to 
this plane. 
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Introducing (82) into the boundary condition (64) and making use of equation 
(66), we obtain the following three linear equations for^, D and F. 

(pco^ - 2fia^ )y4 sinh vji + 2/pa(Z)Vj sinh + Fv, sinh v^h) = 0 
2/av,^coshv,/; + Z)((vj +a^)-^^(vj -a^)'}coshVj/j 
+ f| (vj + ^ ^ (vj - y jcosh v^h = 0 

D (vj - a^)vj sinh v^h + F{yl - a^)v, sinh v^h = 0 


Since i is small, elimination of A, D and F from the equation (83) leads to the 
equation 


tanhv,/? 

tanhVjA 


= 4v,v,{l + (p79^)}/(2a^ -*2){2 -(cVc')} 1 


+ 


tanhVjA^ 
q\ tanhVjA^ 


(84) 


Case I. If the length of the wave is very small compared with the thickness of the 
layer, equation *84) reduces to (79) and yields the velocity of propagation of Rayleigh 
surface waves under the influence of couple-stress. If the parameter of couple-stress is 
zero, the classical result follows. 

Case II : If the length of the wave be large compared with the thickness of the layer, 
quantities vxh, vth and ah can be regarded as small and hyperbolic tangents tanh Vih, 
tanh vjA can be replaced by the first two terms of their expansions into series. Then 
equation (84) becomes 


v.{l-(v.-ftV3)} 


4v.vJl + (/?Vg7} 


I ^VjtanhVjA I 


(85) 


If i = 0 we get, on simplification’ 

(cVca ) = I {l - /c ^ )}, c = (co/a) 

which is the classical result of Rayleigh” and Lamb” since £ is small, we obtain the 
following approximation of equation (85) 



PROBLEMS IN COUPLE-STRESS THEORY OF ELASTICITY 


23 


{Aa^h^l3){\-{cHcl)]+Aa^e ( 86 ) 

whicji shows that the value of dci increase in presence of couple-stress. In other 
words the velocity c of the propagation of waves in the elastic layer increases due to 
the presence couple stress. 

Axis symmetric Lamb's problem in couple stress theory ; 

We consider a time varying loading z{r,t) = pif) exp(/cor) acting on the elastic 
semi-space bounded by the z = 0 plane z axis being pointed into the medium. The 
loading being axially symmetrical, it produces in the semi-space state of stress and 
deformation and the cylindrical co-ordinates is used to investigate the problem. 

We consider the displacement equation of motion in absence of body forces and 
body couple as 

pV^« + (X-fp)W«-t-TiWVxVxM = p« (87) 

we express the displacement u into its lameller and solenoidal components 

u = 'V^ + VxH V.H = 0 (88) 

Then the above equation of motion reduces to the following equation 

c,V> = $ ; c^l - = H (89) 

where = n/p , c' = (^ + 2li)/p , c' = p/p 

For harmonic waves we insert (j) = (j) exp(io)0 and 'H= H exp(/cot)in (89) to obtain 


where 


{v^+o])^ = 0 

CT,=a)/c, 


(w+p^)(w-K)^=o 


(90) 
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<^2 = <0/^2 ; [ (l + 4ealf + if 

Then, considering H = H' + H' the complete solution becomes 

f tf 

u = V^ + VxH +VxH 

t n 

where (j), H and H are governed by the Helmholtz equations 

(v^+a^)$ = 0 ; (V^+P^)^' = 0 ; (v^-pO^'' = 0 


Boundary conditions : The boundary conditions on z = 0 

=-j3(r)e'“' ; Z’, = /*♦ = It,* = M’s- == 0 

where'® 


Pr='C2, 


9 19/3,9 19/3 


dr 


{ 




+- 
^ 2 


r 9(i) 
1 9 


dr 




r 9({) 


/• 9(j) j 

M'ZZ 

r ^ J 


(91) 


(92) 


(93) 


we consider the particular case in which the external loading and the displacement 
vector u are independent of the body forces and body couples being discarded. In 
this case w,. = = 0, Pr = lizr = 0. 

Hence, the boundary conditions reduce to 


2 


1(9 9 

— u, + — p , + 

-^.=-p{r)e'“‘ ; p..* = 0 




= 0 


(94) 


Wkh boundary conditions (94), we seek solutions of (92). In this particular 
case = 0. Hence the equations (92) take the following forms 
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^ . 02') nr' 




72 ^ , 02'! TTir 


Applying in equations of 95), the Hankel transformation defined by 

05 00 

(j)* (a, z) = Jr<i)(r , z)j^ {ar)dr ; \]i* {a, z)= Jrvi7(r, t)Jo {ar)dr 


H'. = — 

* dr^' * 8r ' 
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(95) 


(96) 


we obtain (j)’=rle‘“, \\fl=Ce~'^‘ (97) 

where a = (a^ - a,^ f , v, = (a' - PJ f , v, = (a^ + p' f 
Hence in view of (88), (91) and (97) we obtain 

M, = exp(/©/) J{- O/fe'" + a(v,5e‘''' + v^ce'”'' )} oJ, {ar)da (98) 

0 

= exp(/oc>?)|{- aAe'^ + a (v,5e"''' + ajJ^r)do. (99) 

0 

Substituting (98) and (99) in the boundary conditions (94) we obtain 

2CTrf-{(vi‘ +a^)-^^(vi' -a^)^}a5-{(vi‘ +a^)-£^(vj -a^)^}aC = 0 

(2a^ -0j)^-2a^(v,fi + VjC) = -/j’(a)/p. ► 
av,(a’ - - 2v, (a^ - v’) C = 0 


where we have expressed the loading /i(/')by the Hankel interal 

P(r) = Jp * (a)xJ„ i<xr)da ; P*(a}= ]rP(r)y„ (ar)dr 

0 0 


( 101 ) 


Solving the system of equation (100) we obtain 
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A{a) = - 


KP* 


B{a) = - 


2cj P* 
R{a) \i, ’ 



( 102 ) 


where 



f r2) 


(la^ -a2)-4a\a 

11-% 

,Ji: = (2a^-a^) 1 + -^% 

. ^ \ V, 


Hence for a concentrated force z(r,t) 



aJ, (ar) 
aJ„(ar) 

Jj R(a) 


(103) 


z(r,f) = {PM'-)/2nr} e“' P*(a) = /2k 

Thus the displacement given by (103) under the influence of couple stress being 
know we can determine the stresses and the couple stresses from equation (1). If 
f tends to zero the classical results follow’*. 

Effect of couple-stress on steady-state response to moving loads ; 

We consider a homogeneous, isotropic elastic semi-space selecting the origin on 
the free plane boundary and ;c2-axis pointing into the medium, the Cartesian co- 
ordinates of a point of the semi-space X2 ^ 0 are supposed to be (xi, X2, X3) the 
corresponding displacement components are designated by (wi, 1/2, M3). As regards the 
line load we suppose that on plane boundary of the semi-space there acts the loading 
p5(xi+ut) where P is a constant and 5(t) is the Dirac-delta function. 

We assume that a plane strain state prevails and the elastic displacements mj, M 2 , 
(m 3 = 0) are derivable from the displacement potentials (t>(xi, X2) and X2) so that 

ui =* <l>.i - V.2, U2 = (j),2 + V.i (1 04) 

Displacement equations of motion in absence of body forces and body couples are 
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pVV,+(A, + p)p3,3,+riV^(p^^, 

,-v^p,)=PM, (/,y=l,2,3) 

(105) 


; VV-^Vvj; = c;^V 

(106) 

where cj" = (X, + 2p)/p 


e=r\/M 


Ti is a constant characterzing the existence of couple-stresses and A,, p, are Lame's 
elastic constants. If ^ = Othe couple-stress vanishes and the classical result follov/s at 
once. 

Expressions for stresses and couple-stress in terms of (j) and \[( are 

'tu=l^(2<l>,.2 + M/,u-'K.22)-'nVV ; =p(2(|),, +v(/„-i|/,J+-nVV (107) 

T,3=2riV^(vi;,),p,3=2TiV^(v|/,) 

Now, in the problem stated above, the load moves in the negative direction of xi- 
axis at a constant speed u. An observer moving with the load at the same speed would 
see the load as stationaiy. We introduce a Galilean transformation^’. 

xl = x^+ut, ^ 2=^2 i’ = t ( 108 ) 

then the boundary conditions would be independent of f. For a concentrated line load 
the boundary conditions are in the moving co-ordinates, are 

T 22 = -p5(x;), = P 3 , = 0. (109) 

Now, as the response in the elastic semi-space is in a steady state, <() and i|/ will be 
independent of t' as seen by an observer moving with the load. In other words, xj, 
and / enter <|) and v|/ in combination xi+ wt. Under this assumption, equations (106) 
become 
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V'^-^V> = c:V-^ 

^ dx[\ 

On introducing the Mach numbers^’ 


M\ - ulc\. 

M 2 = m/C2 


(Ill) 

and the parameters. 




P,=a/i-m; 


ifM, <1,M2<1. 

(112) 

P,=^M,'-1 

’ P2=a/^2-1 

ifM,> 1,M2> 1. 

(113) 


K sx". 


( 110 ) 


we obtain the following partial differential equations 


P,' 


2 

dx'i^ 


if M, <1 




2 9V , SV ,2t7M 


i2 


dx[^ dx'i 


_^"V'> = 0 


if M 2 <\ 


and 


dx!^^ axi'^ 


if A/, >1 




2 aV , 9V ^ZryM 


-rv'> = 0 ifM2<l 


(114) 


(115) 


ajcJ^ dx'i 

In view of (1 10) and (1 1 1), the stresses and couple-stresses given by (107) reduce 


into 




‘■12 




'axj 


dx^dx', p 


dx'i 


= M. 


2^ 

dx'i 


-2^.2 


dx'i 


\in=2x]Y' 


^ d^l'^ 
\^dx[j 


dx'idx'2 ’ 


P23 — 2 tiV 


axi'axj 

p ax,ax2 '' dx'i 

aij; 

\^'2J 


( 116 ) 
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The boundary conditions (109) become, at jc' = 0 


.2^ = . 
ac;^ax; 


£ 5 (,) 






= 0 


(117) 


Now, we have to determine the functions (j) and \\i which satisfy the differential 
equations (1 14) or (1 15), the boundary conditions on the free surface (117) and the 
appropriate (1 14) or (1 15), the boundary conditions on the free surface (117) and the 
appropriate radiation and finiteness conditions at infinity. 

The nature of the solution depends on the Mach numbers Mu M 2 . Three cases can 
be distinguished : 

(a) M 2 > M\> 1, (Supersonic) 

(b) M 2 > Ml > 1, (Transonic) 

(c) 1> M 2 > Ml , (Subsonic) 

Since C 2 > Ci, so that M 2 > Mi, the three above exhaust all possibilities. 

Supersonic Case [Mi > 1, A/^ > 1] : 


In this case, we must solve the equations (115), we take the solutions as 

; v(; (118) 


AXx[)=Ae-’"^^ (119) 

A,{x[) = Be-''^’^^ (120) 


where 

Tl. =4'k"-p\ r[,=yJX^-q" 

«’=5^!{i+4(i+K)<"x“}Ui] 


( 121 ) 
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where X, A, & and C are arbitral^ constants. Hence, the solution are 

= = ( 122 ) 

which satisfy the radiation condition at infinity only backward running waves are 
admitted. The other possible solutions, of the form are rejected on the basis of 

the radiation condition because they represent disturbance which originate at infinity 
and coverage towards the load. We assume a general in the form 


vp (;c,' ,x\) = — ]{5(?i) + C(;t) }e'^ dk 

2tt_t 

The boundary conditions (117) give 

(a)^ (mI - 2)^(X)- 2a{T]2B(?t)+ Ti3C(>.)} = ?(?.) 
2r\^A{x)^[Ml -2){5(X)+C(X)}X = 0 


(123) 


(124) 


where ?(x) (= -P/p) is the Fourier transform of-(P/p)8(ar,'') . 
Hence, from (124) we get on simplification 


a(\) = (m]-2) 

where A‘ = {mI - 2)' 


1 + ^^ 
11. 9 j 

/ 


XP‘[X)IK ; Bi^) = -2\P'{X)IK (125) 


Tl. ? . 






rri +T1, 


.2^ 


1+^ 

V rj 


(126) 


In view of (125) and (126) becomes 
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271 r\,q ) A J 

M; {x[, x[) = ~ j _ mL X W. 

27t ^ A r\,a 


Since the parameter of the couple-stress is small, we may neglect its cubes and 
higher powers. Then the approximate forms of (j) and \j/ are follows. 

47tA -1 iX ■ J I 


=-^ t-— + 2ex^A,^—~ + ^ 

2nA .1 ir ' iX^ 2p 


a(xi-Mi) ('l + S^'l 


XP‘{x)dX 


where A = (m,^ - 2)+4p.p, ; ^ ; 2p.(l-fp0 . . 

' P,A ^ 


"■ " 2p,p,/f(x; - PX)}+ 

A. (2 - M.Oi' - px ) - |2PM, - (1 + (^; - px ) - 
p.(i+p:^;|7^.(^:-px) 


«. = ~{p,{2 -px;)-2p,//(x: -px;)}--/^ 


'{2-m:)a,fXx: - PX)+ 4A.i^,{x; -p,x;)- 

P 2 ^1 
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(2 - Ml ){2+{xi }5(^; ) + ^ 5(x: - >/ 


+ 


A 




ax! 


2 + -M,^^ 
. 1 ^ 


\{1-Ml)-A^,^^ |2A, —FM-^2A) + 


2 P 2 ax] j 2 P 2 ax, 


'" = A 


fs / 002 , ^ J ^2 


( 2 -M ,0 2 P^-M.^ 


8 (x:-p,x;)- 4 p.p. 5 (x;-px) 


+L^ 

ji A 


A,(2-Mn£7^(^.'-P')(2P^)[2p^ +^M,^jp: + 


+4 P.P J 2 A. £7 F, (x; - p,x; )| + x; F, (x; - p,x; ) - 


4p.(up;y a 

2 Pj ax; 


7 F,(x:-p,x;) 


T„=-^p,(2-M,^)[5(x;-p,x;)-5(x;-p;)]-^.:^[p.(2-Ar;)A,F,(x:-p,x;)+ 
A |i. A 

+ P,|2A, Air^(x; -P 2 xi) + ^^x^ -P2^2)|(2 -^ 2 ^)] 


where' 


60 


F, = 5">(x: -p,x;) = -(x; -px;)-'5(x; -p.x;) 
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= 6">(x: -PX) = 2(;c: -P.xD 

dx. 


^F, = -p.x;) = -:(x; -p,x,)-8(x: -p,x;),(/ = 1,2) 

ox, 


In case of the above relations (128-132), there is a term containing the parameter 
of couple-stress with the classical part, hence the displacements and stresses are 
slightly changed due to the effect of couple stresses. 

The difference of shear stresses 

t,. = -2nV'V = -4^P,(l + ■^F,(x’ -p.;c;) 

A obc; 


Is significant along the line xi=P 2^2 due to influence of couple stress. In the similar 
way subsonic and transonic cases may be studied. 

Discussion 

The Rayleigh and Love type of surface waves are increased under the action of 
couple stress. As the couple stress parameters 1 is small, the change in the wave 
velocity C is also small. The waves in the layer are also affected by couple stress as 
studied in this review paper. The stress distribution due to the presence of moving 
load in the solid elastic semispace may be considered worthy for numerical 
calculations as the improper integrals have been evaluated and they are written in 
terms of different functions. 

But it is to be mentioned that, this asymmetric theory does not have yet the 
complete or through experimental verification. We know merely the order of 
magnitudes and mutual relation between the material constants in studying most of the 
dynamic problem. However the complete correspondence of the experiment and 
theory exists in the case of discrete media (spatial grillages) where all the material 
constants may be uniquely determined. When passing from such a grillage to a 
continuous medium we obtain exactly Cosserat's continuum. The couple-stress theory 
of elasticity, which is now even considered as an Utopian theoiy by some sections of 
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research workers, will enhance the speed of development in science and technology in 
the future. We must remember that the phenomenon 'Utopian to reality' occurs as a 
repeated feature in the history of science. 
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Abstract 

The Extended Huckel Molecuhir Orbital (EHMO) calculation and the Normal Coordinate Analysis 
(NCA) have been carried out on pentacoordinated nitrosyl adducts of iron and cobalt 
diethyldithiocarbamates with the linear and bent metal-N-0 groups respectively. It is observed that 
the geometry of the metal-N-0 group depends on the symmetries of the HOMOs and LUMOs of NO 
and the parent dithiocarbamate fragments, the energies of metal d orbitals relative to and n* of 
NO as well as the electron density on the metla atom. This observation has been substantiated by the 
force constants of N-0, Co-N, Fe-N, Co-S, Fe-S and S-C bonds evaluated by NCA. 

(Keywords : EHMO/NCA/nitrosobis(diethyldithiocarbamato) iron(II)/nitrosobis(diethyldithiocar- 
bamato) cobalt(II) 


Introduction 

Bis(N, N“dimethyIdithiocarbamato)nitrosyI iron(II), [Fe(NO)(S2CNMe2)2]^ and its 
cobalt analogue, [Co(NO)(S 2 CNMe 2 ) 2 ]^ have the uncommon square pyramidal five 
coordinated configuration as predicted theoretically by Daudel and Bucher^. Although these 
two isomorphous"^ complexes have the same basal ligands yet their metal-(NO) bonding 
systems are found to be different. Unlike the linear Fe-N«0 system in the iron complex, 
the Co-N-O system is significantly bent (~127“)* with unsymmetrical Co(NO) n bond 
which is absent in the former. General discussions^'^ on pentacoordinated nitrosyls revealed 
that the bending of the metal-N-0 group depends upon the different coordination modes 
(NO'^, NO and NO") of NO, the energies of the metal d orbitals relative to nitrosyl and 

n*, electron donating capability of the basal ligands and the position of the metal in the 
periodic table. A comparative molecular orbital (MO) study of the bonding of NO groups 
in two nitrosyl complexes, [Fe(NO)(S 2 CNEt 2 ) 2 ] and [Co(NO)(S 2 CNEt 2 ) 2 ] will hence, be 
of heuristic value. It may be expected that the bonding and the structural differences in the 
M-NO groups (M = Fe or Co) of these complexes should be reflected on the symmetry 
properties of their frontier orbitals as well as the N-O, M-N and M-S force constants. 
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Qualitative MO calculations have been done by Extended Huckel Molecular Orbital 
(EHMO) method and the results thus obtained have been substantiaed by relevant force 
constants evaluated by Normal Coordinate Analysis (NCA). 

Materials and Method 

The complexes were prepared by the method reported earlier^^. The FTIR spectra 
(3000-400 cm“^) of the complexes as KBr pellets were recorded using a Perkin-Elmer 
321(K) FT-IR spectrophotometer. 

EHMO Calculation : EHMO calculation was done with ICON 8 computer program^ ^ 
The mean values of the structural parameters for the complexes (A&B) were taken from 
the published data^’"^ of their methyl analogues assuming that the alkyl substituents have 
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little effect on the distant M-N-0 moiety which is the main centre of interest. Hydrogen 
atoms were excluded to simplify the calculation. The orbital exponents (^), corresponding 
relative weights of d orbitals, the ionization potentials and the Huckel constant {K) 
were obtained from published works'^"''*. The energy levels of the molecular orbitals of 
main interest are furnished in Fig. (1 & 2). 

Normal Coordinate Analysis : The simplified models C and D of the complexes A and 
B were considered for the NCA assuming each ethyl group to behave as a single atom R. 
The 21 normal vibrations (3N-6; N = 9) of each of these 9 atomic models of Cg symmetry 
were found to be within the two symmetry species A' and A" and were classified as 


■ vib 


= 14A' + 7A" 
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The Wilson GF matrix method^^ was used to solve the secular equation using QCMP 
067 computer program^^^. The F matrix was constructed by using Urey-Bradley type force 
field^^. The initial sets of force constants were taken from earlier literature on similar 
system^^ which were refined by the procedure of Nelder and Mead^^ to get the final sets 
of force constants [Table 1(a) & (b)]. The observed and calculated frequencies along with 
the band assignments and potential energy distributions (FED) are furnished in Table 2. 

Results and Discussions 

EHMO : As both the complexes have the same basal dithiocarbamate ligand, discussion 
will cover mainly the bodning scheme of the M-N-0 groups. Cambi^^ reported the cobalt 
complex to be diamagnetic indicating a dsp^ type hybridization. The iron complex with 
one electron less is paramagnetic as confirmed by e.s.r. studies^. Depending upon the 
behaviour of the antibonding orbital (FI*), three different coordinating forms of nitric oxide, 
viz,, NO^, NO and NO” are observed in nitrosyl complexes'^. It is common to consider the 
linear nitrosyls as complexes of NO"^ and bent molecules as complexes of NO“^. Generally 
nitric oxide is like carbon monoxide in forming complexes with transition metals. The extra 
antibonding electron of NO is usually transferred to the metal atom to give the nitrosyl 
cation, Typical complexes are [Fe(CN) 5 NO]^” and [Fe(H20)5N0]^^ The spectral 
and magnetic properties of [Fe(NO)(S 2 CNMe 2 ) 2 ] revealed that nitrosyl bonds formally as 
NO"^ consistent with the observation of an infrared band at 1735 cm”^ usually assigned 
to coordinated NO'^ In their theoretical study, Alderman and Owston^'^ suggested the 
n* (NO) Co electron flow scheme. However, in some complexes, such as 
[Co(NH 3 ) 5 NO]^‘^, in presence of highly basic other Jigands, an electron may be transferred 
from metal to NO^^ The present EHMO calculation shows that the NO groups of 
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Fig. 2 - Interation diagram of [Co(NO)(S 2 CNEt 2 ) 2 ] 
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[Fe(NO)(S 2 CNEt 2 ) 2 ] and [Co(NO)(S 2 CNEt 2 ) 2 ] possess net positive charges of 0.50100 and 
0.36064 respectively indicating more electron donation to the metal atom in the latter 
complex. The fragment molecular orbital (FMO) analysis shows that the unpaired electron 
in the HJOMO (Jt* j) of NO flows to the LUMO of [Fe(S 2 CNEt 2 ) 2 ] fragment under a 
symmetrically and energetically favorable condition (Fig. 3a). Similar interaction is 
observed between the HOMO (7t* p of NO and the LUMO(dj^) of [Co(S 2 CNEt 2 ) 2 ] fragment 



Fig. 3 - HOMO-LUMO interaction in (a) [Fe(NO)(S2CNEt2)2] and (b) [Co(NO)(S2CNEt2)2]. 

in the cobalt complex (Fig. 3b). It is clear from the Fig. 3 (a & b) that the iron atom forms 
n bond with NO through the N atom only while in the cobalt analogue both the N and O 
atoms of the nitrosyl group are involved in tc bonding interaction with the metal orbital 
fulfilling the criteria for a bent Co-N-0 system. It was suggested^ that higher the energy 
of the metal d^^ (or dy^) and d^ relative to nitrosyl n* and a„ respectively greater is the 
M-N-0 bending. The present work has revealed that the metal d^^ (or dy^ and orbitals 
are more close in energy to n and a„ of NO in the cobalt complex than in the iron analogue 
(Fig. 1 & 2) suggesting a greater Co-N-0 bending than that of Fe-N-0 group. Mingos^ 
also reported that the bending of M-N-0 group in [Co(NO)(NH 3 ) 5 ]^'*' influences substantially 
the energies of mainly two molecular orbitals containing metal d^^ and and nitrosyl a„ 
and 7t*. 

Normal Coordinate Analysis : The IR spectra of [Fe(NO)(S 2 CNEt 2 ) 2 ] and 
[Co(NO)(S 2 CNEt 2 ) 2 ] were recorded in the range (30(X)-4()0 cm"'). The bands above 2000 
cm”' are not assigned as there were no significant absorptions except C-H stretching^^ and 
because hydrogen atoms were not included in the present calculation. The bands observed 
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Table 1 - Refined sets of force constants (mdyn/A) for (a) [Fe(NO)(S2CNEt2)2] and (b) [Co(NO)(S2CNEt2)2]. 


(a) fi(Fe-S) = 3.98 
f2(S-C) = 6.08 
f3(C-N) = 7.92 
f4(N-R) = 7.86 
f5(Fe-N) = 1.84 
f^(N-0)= 12.26 
fyCSFeS) = .54 
fgCFeSC) = .08 
f9(SCS) = .60 
flo(SCN) = .28 
fn(CNR) = .18 
fl2(RNR) = .38 
fl 3 (SFeN) = .72 


f j 4 (FeNO) = .68 
fl5(S ... S) = .80 
fl6(Fe...C) = .ll 
fl7(S...N) = .24 
fl8(C...R) = .50 
fl9(R ... R) = .40 
f2o(S ... N) = .03 
f2i(Fe ... 0) = .25 


(b) fi(Co-S) = 3.50 
f2(S-C) = 6.25 
f3(C-N) = 7.65 
f4(N-R) = 7.77 
f5(Co~.N) = 1.90 
f6(N-0)=12.2 
f 7 (SCoS) = .56 
fgCCoSC) = .08 
f9(SCS) = .55 
flo(SCN) = .45 
fn(CNR) = .15 
fl2(RNR) = .44 
fl 3 (SCoN) = .74 


fl 4 (CoNO) = .62 
fl5(S...S) = .78 
fl 5 (Co...C) = .i2 
fi7(S...N) = .27 
fl8(C ... R) = .43 
fl9(R ... R) = .33 
f 2 o(S...N) = .05 
f 2 i(Co...O) = .23 


Table 2 - Comparison of observed and calculated frequencies (cm‘*) and band assignments with potential energy 
distributions (FED). 


[Fe(NO)(S2CNEt2)2] 

[Fe(NO)(S2CNEt2)2] 

Obs 

Cal 

Assignment 

FED (%) 

Obs 

Cal 

Assignment 

FED (%) 

1687 

1687 

v(N-O) 

95 

1 

1682 

1673 

v(N-O) 

97 

1550 

1502 

'^.sym (C-N) 

69 

1495 

1497 

Vsym(C-N) 

67 

1272 

1265 

'^a.sym (N”R) 

94 

1269 

1253 

'^asym (N—R) 

92 

1067 

1074 

Vasym (S-C) 

94 

1071 

1083 

'^asym 

92 

967 

960 


61,35 

997 

979 


57, 37 

— 

686 

< 

3 

68 

649 

647 

Vsym (Co-S) 

60 

555 

554 


41,26 

551 

550 

Vsym(S-C)+Vsy„(Co-S) 

31, 17 

— 

— 


— 

504 

508 

|^(Co— N) 

30, 29 

477 

477 

^a.sym 

99 

— 

437 

^asym (Co— S) 

99 


379 

'^sym 

82 

— 

403 

'^sym (Co— N) 

46 

— 

212 

5sym(FeNO)+NOrock 

23,65 

— 

190 

5syn,(CoNO)+NOrock 

35, 53 
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at 1349.6 cm“^ and 1197.1 cm""^ in the iron complex and those at 1350.2 cm“^ and 1210 
in the cobalt complex are probably due to C-C stretching of the ethyl groups and were 
not calculated here. The bands at 1450.2 cm"^ and 1376.5 cm-^ in the cobalt complex and 
at 1419.1 cm"^ and 1375 cm~^ in the iron complex are comparable to those reported 
for 8 j|qj^(CH 3) and were excluded from calculation due to the same reason mentioned 
above. Similarly the bands at 848.7 cm'^ in the cobalt complex and that at 843.8 cm*"^ in 
the iron complex might be due to the C-C-C bending coupled with the C»C and C-CH3 
stretching which were not included in the present calculation. Eleven vibrational frequencies 
for each complex were calculated. The agreement between the observed and the calculated 
frequencies are quite satisfactory (Table 2); the root mean square errors being 4.7 and 10.1 
cm"^ for the iron and the cobalt complexes respectively. The IR assignments and the force 
constants for the basal ligands are consistant with those reported earlier for dithiocarbamate 
complexes The bands at 1687 cm"** and 1682 cm“^ are assigned to for 

[Fe(NO)(S2CNEt2)2] and [Co(NO)(S2CNEt2)2] respectively and are comparable to their 
corresponding methyl analogues^. However, these values are much less than the vibrational 
frequency of free NO which is about 1876 cm“^ indicating a decrease in the N-0 bond 
order in these two complexes. It is reasonable to expect a 7i-back bonding from and 

dy^ orbitals of the metal to the n* (NO)^^ reducing the N«0 bond order. Such a considerable 
7S-back bonding in most nitrosyl complexes is evidenced by the relatively shorter M-N 
bonds in these systems^^. The flow of electrons from 0^^ orbitals of NO to metal orbital 
also contributes to the decrease of N>0 bond order. The combined effects of sigma mode 
of complexing of the NO"^ ion as a Lewis base and the it-back bonding predominate the 
transfer of electrons between metal and NO molecules resulting in a net decrease of N-O 
bond order. As the iron atom has one valence electron less than the cobalt atom so the 
extent of electron flow from 71* of NO to the iron d orbital would be more and hence the 
N-0 bond order in the iron complex becomes slightly greater than that in the cobalt complex. 
It is reflected on the greater (1687 cm"^) and N-0 stretching force contant = 
12.26 mdyn/A) in [Fe(N0)(S2CNEt2)2] than the corresponding values (v^„^ = 1682 cm~^ 

12.20 mdyn/A) in [Co(NO)(S2CNEt2)2]. Eamshaw etaL^ reported that the electron 
removal from the cobalt atom by electron withdrawing substituents decreases the dn-pn 
bonding between cobalt and nitric oxide and increases the N-0 bond order. It suggests that 
with the same electron donating alkyl groups in the iron and the cobalt complexes the 
excess electron in cobalt will lead to a greater Co-NO dn-pn bonding than that between 
Fe and NO with the expected force constant order . Present observation 
substantiates this with respective values of/^^^= 1.90 mdyn/A and/^^^= 1.84 mdyn/A. 
The bands at 1495 cm*^ and 1500 cm'* have been assigned to partially double bonded C-N 
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stretches and those observed at 1269 cm"' and 1272 cm"' are assigned to Vjy_^ which are 
well within the range (1480-1550 cm"' for ^nd 1250-1350 cm-1 for observed 
by Chatt et As suggested by Cotton and McCleverty^^, the relatively higher 

(1500 cm"') in the iron complex than that (v^./v = 1495 cm"') in the cobalt complex might 
be the result of greater electron donation by the sulfur atoms to the iron atom which is 
ahead of cobalt in the position of the periodic table. The above explanation is substantiated 
by the following order of C-N and S-C stretching force constants : 

Fc_yv(Fe) = 7.92 mdyn/A > fc_N(Co) = 7.65 mdyn/A and 
fj_^Fe) = 6.08 mdyn/A < fs_<<Co) = 6.25 mdyn/A 

Conclusion 

The reason for linear and bent structure for nitrosobis (diethyldithiocarbamato) iron 
(II) and nitrosobis (diethyldithiocarbamato) cobalt (II) has been proposed on the basis of 
HOMO/LUMO interaction of NO and metal dithiocarbamates, derived from EHMO 
calculations. The NCA of the above systems corroborate the observations obtained from 
EHMO analysis. 
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Abstract 

Potentiometric equilibrium measurements have been made at 30“C for the formation of ternary 
complexes of 5-hydroxy-2-(hydroxy methyl)-4(H)pyran-4-one (known as kojic acid) or 3- 
hydroxy-2-methyl-4(H)-pyran-4-one (known as maltol) and Co(II), Ni(ll), Cu(II) and Zn(II) 
with imidazole, 1 -methyl imidazole, 2-methyl imidazole, 2-ethyl imidazole and benzimidazole 
in 1:1:1 ratio. The comparison of AlogAT values reveals that the order of stability of ternary 
complexes is: 2-methyl imidazole S 2-ethyl imidazole > imidazole a 1 -methyl imidazole > 
benzimidazole. 

(Keywords : potentiometry/mixed ligand complexes/substituted imidazoles/stability constants). 

Introduction 

The binding of imidazole with transition metals is of interest because of this 
ligand’s close relationship with biological systems of a more complex nature 
involving histidine residues. The literature contains free energy data obtained by 
potentiometric studies on the systems.''^ These data indicate that imidazole forms 
some of the most stable complexes of all heterocyclic -N ligands. This is corroborated 
by the high acid dissociation constants of imidazoles. It is also very well established 
that in the combination of cadmium and copper ion with serum albumin the principal 
binding sites are the imidazole groups of the histidine residues of the protein 
molecule. Several studies have been made on the imidazole ring to determine which 
nitrogen is involved in the binding to metals. Li et af determined the constants of 
imidazole and N-methyl imidazole with Cu(II). By their similarity in formation 
constants the binding site was assumed to be the pyridine nitrogen. 

Owing to the biological importance of imidazole, we have undertaken the study of 
substituted imidazoles with transition metals. In continuation of our earlier studies,*’’ 
we now report the stabilities associated with the metal complexes of substituted 
imidazoles with bivalent metals in the presence of oxygen donors such as 5-hydroxy- 
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2-(hydroxy methyl)-4(H)pyran-4-one (known as kojic acid) and 3-hydroxy-2-methyl- 
4(H)-pyran-4-one (known as maltol). 

Materials and Method 

The ligands, kojic acid and maltol were obtained from Lancaster (UK). Other 
ligands, imidazole and substituted imidazoles were obtained from Fluka 
(Switzerland). The metal salts of Co(II), Ni(II), Cu(II) and Zn(II) were of Anala R 
grade, and the metals were standardized volumetrically by titration with the sodium 
salt of EDTA.® Potentiometric measurements were carried out at 30°C with O.IM 
KNO3 as background electrolyte in aqueous medium using a control dynamics pH 
meter. The pH meter readings were plotted against a (moles of base added per mole of 
ligand) or m (moles of base added per mole of metal ion). Calculations were made 
with the help of BEST computer programme.® 

Results and Discussion 

Acid dissociation constants : The acid dissociation constants of kojic acid, maltol, 
imidazole 1 -methyl imidazole, 2-methyl imidazole, 2-ethyl imidazole and 
benzimidazole were remeasured and the values obtained are in good agreement with 
the values found in literature. 

Binary Constants : The formation constants for 1:1 [M-L] system where M = 
Ni(n), Co(II), Cu(n), Zn(II) and L = maltol, kojic acid, imidazole, 1-me imidazole, 2- 
me-imidazole, 2-ethyl imidazole and benzimidazole were measured in aqueous 
medium at 30°C and p, = O.IM KNO3. The values obtained (Table 1) are in good 
agreement with literature values.®’’ 

Table 1-Acid dissociation constants of ligands and their stability constants with bivalent transition metals. 


{T = 30°C; 7= O.IM KNO 3 ; aqueous medium] 


Ligands 

pK. 

Co(n) 

Ni(II) 

Cu(0) 

Zn(n) 

Imidazole 

7.01 

2.56 

3.21 

4.28 

2.64 

1-Me-imi 

7.19 

2.61 

3.32 

4.25 

2.71 

2-Me-imi 

7.93 

2.71 

3.38 

4.38 

2.93 

2-Et-imi 

7,95 

2.82 

3.41 

4.41 

2.97 

Bzm 

5.37 

2.25 

2.65 

3.28 

%35 

Maltol 

8.55 

5.35 

5.57 

7.61 

5.11 

Kojic acid 

7.62 

4.75 

4.93 

6.57 

5.14 


Standard deviations (a fit) are within ± 0.04 log units. 
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Ternary Constants : The mixed ligand titration curves of all the systems showed a 
feeble inflection at w = 2. Mixed ligand titration curve of Ni(II) with maltol in 
presence of imidazole is given in Fig. 1 along with other curves (free ligand and 
binary). These curves are related to the usual equilibrium. 



Fig. i-Potentiometric titration curves of the ternary complexes of [maltol-Ni-imidazole] A = free ligand 
(imidazole), B = Ni-imidazole, C = free ligand (maltol), D = Ni-maltol, E = 1 : 1 : 1 
(maltol-Ni-imidazole) ternary complexes. 

Calculations were made with the help of BEST computer programme. Out of the 
various models, one which gave the best statistical fit was selected. The binary 
stability constants were used as fixed parameters for the refinement of the ternary 
formation constants. The values thus obtained are presented in Table 2. 

The stability of the various ternary metal complexes and corresponding binary 
metal complexes have been quantitatively compared in terms of the parameter, A log K 
which is expressed as 
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Alog^ = log^l-log^“ (2) 

The AlogA:^ values for the different ternary complexes have been calculated and 
are listed in Table 2 along with ternary constants for different systems. In order to 
obtain precise values of Alog^, it is essential that the binary and ternary constants 
should refer to identical experimental conditions. Therefore acid dissociation 
constants of all the ligands under study and their binary constants have been 
remeasured. 


Table 2- Logarithms of ternary stability constants of bivalent metals with substituted imidazoles in 
* presejvce of maltol or kojic acid and their corresponding AlogK values. 

[T = 30®C; /= O.IM KNO 3 ; aqueous medium] 



Co(II) 

Ni(n) 

Maltol 

Cu(II) 

Zn(n) 

Co(II) 

Kojic Acid 

Ni(II) Cu(II) 

Zn(II) 





Ternary constants 





Imidazole 

7.99 

8.88 

12.07 

8.50 

7.29 

8.18 

10.92 

7.80 

1-Me-imi 

8.02 

9.01 

12.03 

8.58 

7.30 

8.31 

10.95 

7.88 

2-Me-imi 

8.14 

9.10 

12.23 

8.88 

7.42 

8.42 

11.08 

8.20 

2-Et-imi 

8.27 

9.16 

12.27 

8.93 

7.59 

8.42 

11.15 

8.22 

Bzm 

7.57 

8.28 

10.93 

8.14 

6.92 

7.54 

9.90 

7.43 





AlogA: 





Imidazole 

+0.08 

+ 0.10 

+0.18 

+0.09 

- 0.02 

+0.04 

+0.07 

+ 0.02 

1-Me-imi 

+0.06 

+ 0.12 

+0.17 

+ 0.10 

-0.06 

+0.06 

+0.13 

+0.03 

2-Me-imi 

+0.08 

+0.15 

+0.24 

+0.18 

-0.04 

+ 0.11 

+0.13 

+0.13 

2-Et-imi 

+ 0.10 

+0.18 

+0.25 

+0.19 

+ 0.02 

+0.08 

+0.17 

+ 0.11 

Bzm 

-0.03 

+0.06 

+0.04 

+ 0.02 

-0.08 

-0.04 

+0.05 

-0.06 


Standard deviations are within ± 0.05 log units. 
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Perusal of the AlogA^ values (Table 2 ) shows that the stability order with respect to 
the secondary ligands in the ternary complex is 2-methyl imidazole > 2-ethyl 
imidazole > imidazole > 1 -methyl imidazole > benzimidazole. This sequence may be 
attributed to the decrease in the basicities of the ligands. For steric consideration, the 
AlogK values in the case of complexes with substituted imidazoles should have been 
smaller than that for imidazole complexes. The increase in basicity in substituted 
imidazoles, thus compensates for the negative contribution of the steric effect and 
further enhances the values of Alog/C in all the cases. 

Amongst the substituted imidazoles, there is not much difference either in the 
ligand basicities or in the steric effect so that both 2-methyl imidazole and 2-ethyl 
imidazole gave almost similar Alog^ values. The AlogAT values of imidazole and 1 - 
methyl imidazole are also comparable, may be due to their similarities in the basicity. 
Replacement of hydrogen by CH3 group at one position in 1 -methyl imidazole has 
little effect. The lower stabilities of ternary complexes involving benzimidazole is in 
line with basicity order. 

The positive AlogX values of these ternary systems may be attributed to the dn-pn 
back bonding. Like 2 - 2 ’-bipyridyl, 1 , 10 -phenanthroline, imidazole or substituted 
imidazoles are also reasonably good n acceptors. Therefore, imidazole or substituted 
imidazoles would prefer binding with 0 - 0 " donors rather than N-N donors.^ Maltol 
or kojic acid being 0-0 donors prefer binding with imidazole or substituted 
imidazoles, hence positive or less negative AlogK values are obtained. 

In other words, the extra stabilization of ternary complexes can be explained by 
considering that there is a M(II)->L, n interactions and hence the electron density on 
the metal ion is reduced and the effective electronegativity of [M-L]"^ complex is 
almost the same as that of [M(H20)6]^'^. Thus a bonding tendency of imidazole or 
substituted imidazoles to combine with [M-L]"" and [M(H20)6]^'^ is not significantly 
different. Hence ternary values are more than binary values resulted in positive A\o%K 
values. The order of stability in terms of O-O" donors is as maltol > kojic acid. This 
order is in line with basicity order of these two ligands. The order of stability of 
ternary complexes in terms of metal ions is in conformity with Irving Williams natural 
order of stabilities.'® 
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Abstract 

Stability constants of ternary complexes of inorganic phosphates (viz., orthophosphate (OP), 
pyrophosphate (PP) and tripolyphosphate (TPP) with bivalent metal ions in presence of 
benzimidazole (Bzm) and substituted benzimidazoles (viz. 2-Methyl benzimidazole (Me-Bzm), 
2-Propyl benzimidazole (Pro-Bzm), 5,6-Dimethyl benzimidazole (5,6-Me-Bzm), 2- 
Hydroxymethyl benzimidazole (HMB), 2-Hydroxethyl benzimidazole (HEB)) have been 
determined potentiometrically at 30"C and p~0.1M KNO 3 in 20% methanol medium. Ternary 
stability constants of the various inorganic phosphates were compared. It was found that the 
stability of the ternary complexes in terms of inorganic phosphates decrease in order : TPP > 

PP > OP. The ternary complexes of TPP are more stable than PP due to more electrostatic 
attraction between TPP'^ and positive charge on coordination sphere of the metal ion. The 
ligand OP has comparatively less negative charges than PP (L^‘) or TPP (L^'), hence it 
forms least stable complexes. Further, OP may act as monodentate ligand, even if it behaves as 
bidentate, the resultant 4-membered chelate would be less stable. 

(Key words : pH metry/mixed ligand complexes/bihary complexes/temary complexes/stability 
constants). 


Introduction 

The interplay between nucleotides and metal ions is now well established for 
many enzymatic systems^ Consequently, metal ion-nucleotide interactions^ have been 
studied in the solid state^ and in solution"* and our knowledge on binary metal 
complexes of nucleotides 5'-monophosphates or 5'-triphosphates is considerable. The 
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stability^ of these complexes and their structures'* in solution are well characterized, 
and the tendency for macrochelate formation in monomeric metal ion complexes of 
nucleoside 5 '-triphosphates and the promotion of stacking by metal ions have been 
quantified*-’-®. In continuation of our earlier studies’- on stabilities of the ternary 
complexes involving nucleic acid constituents, we now report the stabilities associated 
with the ternary complexes involving benzimidazoles and inorganic phosphates and 
bivalent metals. 


Materials and Methods 

The ligands Bzm and 5,6-Me-bzm were obtained from Sigma Chemical 
Company(USA) and John Baker Inc., Colorado, USA respectively. The other ligands 
under study were prepared in this laboratory according to the procedures available". 
The other secondary ligands, orthophosphate (OP), tetra sodium pyrophosphate (PP) 
and pentasodium tripolyphospyhate (TPP) were obtained from Aldrich Chemical 
Company Ltd., (USA). 

All bivalent metal ions, EDTA, potassium hydrogen pthalate, KNO3, NaOH were 
of BDH Analar grade. A stock solution of bivalent metal ions were prepared and 
standardized by the various methods available in literature"- Carbonate-free sodium 
hydroxide was prepared and standardized by titrating with potassium hydrogen 
pthalate. Acid dissociation constants of the free ligands and formation constants of 
metal complexes were determined by potentiometric titrations of the various ligands 
and metal-ligands mixtures (in ratio 1:1 and 1:1:1) with standard carbonate free 
NaOH. The pH measurements were carried out at 30 ®C with O.IM KNO3 as 
background electrolyte in 20% (v/v) methanol-water medium using a control dynamic 
pH meter fitted with gel filled combined electrode assembly. The pH meter readings 
were plotted against a (moles of base added per mole of ligand) or m (moles of base 
added per mole of metal ion). Calculations were carried with the help of computer 
programe - BEST". 


Results and Discussion 


Acid dissociation constants 

(a) Benzimidazole and substituted benzimidazoles 

The free ligand titration curves for the monoprotonated (added one mole of extra 
acid) form of all the ligands under study, show a steep inflection at a = 1 (where a = 
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moles of base added per mole of ligand). The dissociation constants (pK^) were 
calculated in the region a = 0 to a = 1 . For ligands HEB and HMB apart from pK^, a 
second dissociation constants pK2a were also calculated in the region a = 1 to a = 2. 
The dissociation constants for all the ligands were computed from the potentiometric 
pH profiles of their solutions in the absence of metal ions. The negative logarithms of 
dissociation constants are listed in Table 1 . 


'able 1- Acid dissociation constants of different ligands under study (T = SO^C ; g = 0.1 M KNO3 ; 20% 
methanol) 


Ligands 

pKa 

PK2a 

Ligands 

pKa 

PK2a 

Imidazole 

6.94 


Benzimidazole 

5.37 


Me-Bzm 

5.95 


Pro-Bzm 

6.02 

- 

5,6-Me-Bzm 

6.12 


HMB 

5.40 

11.21 

HEB 

5.47 

12.00 

OP 

6.71 

10.62 

PP 

6.32 

8.50 

TIP 

5.61 

8.10 


The basicity order of ligands, listed in Table 1 , increase in the order Bzm ~ HMB 
< HEB < Me-Bzm < Pro-Bzm < 5 , 6 -Me-Bzm < Imi. The basicities of Bzm, HMB 
and HEB are almost similar. The presence of hydroxyl methyl or hydroxy ethyl 
groups at position 2 in the benzimidazole moiety has no effect on basicity of tertiary 
nitrogen of benzimidazole. The pK2a (dissociation of proton from hydroxyl group) of 
HMB, HEB is found to be above 11 . The pK^ values obtained in the present 
investigation agree closely with those recorded in the literature^ Introduction of 
electron donating groups at position 2 (e.g. 2 -Me-Bzm, 2-Et-Bzm, etc) or at positions 
5 and 6 (e.g. 5 , 6 -dimethyl benzimidazole) makes the tertiary nitrogen of 
benzimidazole more basic. 

(b) Inorganic Phosphates 

The free ligand (KH2PO4) titration curve of orthophosphate (OP) show an 
inflection at a = 1 , followed by a buffer region in higher pH. The pk^ and pk2a are 
calculated in the region a = 0toa=l&a=ltoa = 2 respectively. The free ligand 
curve of pyrophosphate (PP) taken in diprotonated form (Na2H2P207) shows a small 
inflection at a = 1 followed by a buffer region. Similarly, tripolyphosphate (TPP) in 
triprotonated form (Na3H2P30]Q) shows two moderate inflection at a = 1 and a = 2. 
Acid dissociation constants are calculated using computer program PKS and the 
values thus obtained are in good agreement with literature values 
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Binary Metal Complexes 

(a) Benzimidazole and substituted benzimidazoles 

The potentiometric titration curves for the binary 1:1 ligand systems shov^^ a 
single steep inflection at m = 1 (where m = moles of base added per mole of metal 
ion). The formation constants (log for the binary complexes formed were 

calculated and the values thus obtained are given in Table 2. The potentiometric 
titration curves for the binary 1:1 [M^''-HMB] or systems show a single 

inflection at m = 1, indicating a formation of protonated complexes. 


Table 2 - Logarithms of stability constants of binary complexes of substituted benzimidazoles [T = 30^0; 
ji = 0.1 M KNO3; 20% methanol] 


Ligands 





Imidazole 

2.56 

3.21 

4.28 

2.64 

Bzm 

2.25 

2.65 

3.28 

2.35 

Me-Bzm 

2.35 

2.73 

3.46 

2.51 

5,6-Me-Bzm 

2.45 

2.82 

3.63 

2.87 

Pro-Bzm 

2.19 

2.56 

2.89 

2.27 

HMB 

2.20 

2.23 

3.18 

2.15 

HEB 

2.25 

2.17 

3.25 

2.15 


The constants are accurate to ± 0.05 logK units 

The values of stability constants of HMB and HEB are log K 


The results of the present study (Table 2) show that bzm and its alkylated 
derivatives form less stable coordination compounds than those of the imidazoles. 
This is the result of two factors (i) the pK^ of the imidazole is greater than those of the 
benzimidazole / substituted benzimidazoles under study, (ii) The bulky benzene ring 
attached to the imidazole system in the benzimidazoles probably decreases their 
stabilities. Freiser^^ has also suggested that the benzene portion of the benzimidazole 
or 2-substitued benzimidazoles interfere sterically with chelation. 

(b) Benzimidazole vs. substituted benzimidazoles : 

The stability constants listed in Table 2 increase in the order: Pro»Bzm < HMB - 
HEB < Bzm < Me-Bzm < 5,6-Me-Bzm < Imi for Cu^^ metal ion. This sequence is in 
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accordance with the basicity order of these ligands except 2-propyI benzimidazole, 
which show lower stability though it has higher basicity comparatively. This anomaly 
is due to the steric hinderance caused by the propyl group present at position 2 in 2- 
propyl benzimidazole. Generally, an alkyl group in sterically blocking positions only 
decreases the stability of the chelates^^. Lane^^ showed that benzene portion of the 
benzimidazole does not cause any significant steric hinderance in planar chelates 
provided that the coordinating group in the 2-position is small e.g. -CH 2 NH 2 ,- 
CH20H,-C00H etc. However, if a bulky group is coordinating, such as the methyl 
amino methyl, or the o-hydroxy phenyl, this group may be sufficiently large to 
overlap with the hydrogen in the 4-position and thus cause steric interference. 

(c) Inorganic Phosphates 

(i) M^*-Orthophosphate: Potentiometric titration curves of M (Il)-orthophosphate 
(diprotonated) in 1 : 1 ratio shows a inflection at m = 1 followed by a buffer region. 
The separation of ML curve from the free ligand curve is small suggesting complex 
formation with low stability constants. The stability constants obtained are presented 
in Table 3 and are in good agreement with literature values’’. 

Orthophosphate can bind to the metal ion with two of its oxygen atoms to form 
four-member ring leaving another oxygen unbound (whose pK^ is more than 10). But 
the experimentally observed values (table 3), which are very low, suggest 
monodentate type of binding using the oxygen (whose pK^ is 6.71), leaving the 
oxygen (pK^ = less than 2) and another oxygen (pK^ > 10) unbound. 

(ii) M’'"-Pyrophosphate: The ligand pyrophosphate has four oxygen’s, which can 
coordinate to metal ion. Theoretically it is possible to form two four membered and 
one six membered rings using four oxygen’s simultaneously, or with two middle 
oxygen atoms, it can form a single six membered ring, leaving other two oxygen’s 
uncoordinated to metal ion. 

Potentiometric titration curves of M’* -pyrophosphate in 1 : 1 ratio shows a single 
inflection at m = 2. pH of metal-ligand titration curve is much lower than the free 
ligand curve indicating formation of strong complex with high formation constants. 
Potentiometric titration curves for M’^-pyrophosphate in 1:2 ratio gave a feeble 
inflection at m = 3 and another one at m = 4. The stability constants of ML 2 are given 
in Table 3. 

(iii) M’"'-Tripolyphosphate: Potentiometric titration curve of M2+-tripolyphosphate 
(in triprotonated form) in 1:1 ratio shows an inflection at m = 3. Metal ligand titration 
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curve is much more lower than the free ligand curve indicating strong complex 
formation. The ML complexes are quite stable with large formation constants. 
However, it is observed that the stability constants of [M-TPP] are lower than the 
stability constants of [M-PP], though the fact that TPP is a tridentate and would expect 
to form much more stable complexes. 

The stability constants for the mono binary complexes of OP, PP, TPP show that 
the stability decreases in the order: PP > TPP > OP. 


Table 3- Logarithms of stability constants of binary complexes of inorganic phosphates [T = 30 “C, p = 
0.1 M KNO 3 ; 20% methanol] 


Metal ions 

OP 

MHL 

MHL 

PvroDhosohate 

ML 

ML 2 

MHL 

TrioolvDhosDhate 

ML 

ML 2 

Co 2 + 

2.31 

3.50 

6.45 

2.50 

4.92 

6.05 

2.71 

Ni 2 + 

2.32 

3.83 

6.98 

1.60 

4.83 

6.64 

2.60 

Cu 2 + 

2.96 

5.38 

8.60 

4.86 

5.96 

7.94 

2.80 

Zn 2 + 

2.42 

3.80 

ppt 

ppt 

5.10 

7.67 

2.52 


The order of stabilities of 1:1 complexes of these ligands with respect to metal 
ions [Cu^"^ > NP > Zn^"^ > Co^^] is in agreement with Irving-William’s order^''. 


Ternary Metal Complexes 

(a) + P04' + Substituted benzimidazoles: The potentiometric titration curves 
of the teraaiy complexes involving P04' and substituted benzimidazoles show a 
single inflection at m = 2 indicating the simultaneous coordination of P04‘ and 
substituted benzimidazoles to the metal ion. The ternary constants are calculated in 
the buffer region m = 0 to m = 2 using the pKa of P04’ and pKa of substituted 
benzimidazoles. The values thus obtained are presented in Table 4 . 

(b) + P207^ + Substituted benzimidazoles: The potentiometric titration curves 
of the ternary complexes in 1 : 1:1 ratio shows an inflection at m = 2 indicating the 
simultaneous coordination (formation of protonated complexes). The constants are 
calculated in the buffer region m = 0 to m = 2 , using the pKa of P207^' and pKa of 
substituted benzimidazoles. The values then obtained are presented in Table 4 . 
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(c) + PjOio^ + Substituted benzimidazoles: The potentiometric titration 

curves of the ternary complexes in 1:1:1 ratio shows an inflection at m = 3 indicating 
the simultaneous coordination. The constants are calculated in the buffer region m = 
0 to m = 3. The values thus obtained are presented in Table 4. 

Table 4- Logarithms of stability constants 6 '^mal ternary complexes of bivalent metal ions 
[T = 30®C: p = 0. 1 M KNOj; 20% methanol medium] 


Systems 

Co^'^ 





L 

= Orthophosphate (OP) 



Bzm 

4.48 ± 0.05 

4.99 ± 0.03 

6.26 ±0.03 

4.69 ±0.07 

Me-Bzm 

4.61 ±0.04 

5.07 ± 0.04 

6.43 ± 0.06 

4.85 ± 0.03 

5,6-Me-Bzm 

4.74 ± 0.08 

5.22 ± 0.06 

6.64 ± 0.08 

5.27 ±0.08 

Pro-Bzm 

4.41 ±0.02 

4.84 ± 0.05 

5.80 ± 0.06 

4.57 ± 0.09 

HMB 

4.41 ±0.08 

4.44 ± 0.02 

5.99 ± 0.05 

4.44 ±0.01 

HEB 

4.44 ± 0.09 

4.36 ± 0.04 

6.04 ± 0.06 

4.45 ± 0.02 


L = 

= Pryophosphate (PP) 



Bzm 

5.77 ±0.04 

6.52 ± 0.06 

8.72 ± 0.03 

6.17 ±0.08 

Me-Bzm 

5.91 ±0.03 

6.62 ± 0.03 

8.92 ± 0.06 

6.37 ±0.06 

5,6-Me-Bzm 

6.03 ± 0.04 

6.76 ± 0.06 

9.13 ±0.05 

6.74 ± 0.05 

Pro-Bzm 

5.71 ±0.08 

6.40 ±0.01 

8.37 ± 0.09 

6.11 ±0.03 

HMB 

5.65 ±0.01 

6.09 ±0.05 

8.62 ± 0.06 

5.92 ± 0.04 

HEB 

5.69 ± 0.06 

6.02 ± 0.06 

8.67 ±0.08 

5.89 ± 0.02 


L = 

Tripolyphosphate (TPP) 



Bzm 

7.21 ±0.05 

7.54 ±0.08 

9.32 ± 0.04 

7.49 ± 0.05 

Me-Bzm 

7.35 ± 0.04 

7.64 ± 0.06 

9.54 ± 0.05 

7.70 ± 0.06 

5,6-Me-Bzm 

7.49 ±0.03 

7.80 ± 0.03 

9.80 ± 0.08 

8.10 ±0.08 

Pro-Bzm 

7.16 ±0.02 

7.43 ± 0.02 

8.93 ± 0.06 

7.43 ± 0.03 

HMB 

7.15 ±0.09 

7.11 ±0.08 

9.21 ±0.03 

7.29 ±0.01 

HEB 

7.19 ±0.08 

7.04 ± 0.06 

9.27 ± 0.02 

7.30 ±0.04 
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The relative stability of the ternary metal complexes (MLA) as compared to that 
of the correspondent binary complexes have been quantitatively expressed in terms 
the parameter AlogK given by the following equation. 

AlogK = [logK V a - (logK V + logK “ml)] 

The obtained values of AlogK are reported in Table - 5 


Table 5- Comparison of AlogK values associated with ternary complexes 
[T = 30^C : p = 0.1 M KNO^; 20% methanol medium] 


Systems 

Co^^ 


Cu^^ 

Zn^-*- 



= Orthophosphate (OP) 



Bzm 

-0.08 

+0.02 

+0.02 

-0.08 

Me-Bzm 

-0.05 

+0.02 

+0.01 

-0.08 

5,6-Me-B2m 

-0.02 

+0.08 

+0.05 

-0.02 

Pro-Bzm 

-0.09 

-0.04 

-0.05 

-0.12 

HMB 

-0.10 

-0.11 

-0.15 

-0.13 

HEB 

-0.12 

-0.13 

-0.17 

-0.12 


L 

= Pryophosphate (PP) 



Bzm 

+0.02 

+0.04 

+0.06 

+0.02 

Me-Bzm 

+0.06 

+0.06 

+0.08 

+0.06 

5,6-Me-Bzm 

+0.08 

+0.11 

+0.12 

+0.07 

Pro-Bzm 

+0.02 

+0.01 

+0.10 

+0.04 

HMB 

-0.05 

+0.03 

+0.06 

-0.03 

HEB 

-0.06 

+0.02 

+0.04 

-0.06 


L = 

Tripolyphosphate (TPP) 



Bzm 

+0.04 

+0.06 

+0.08 

+0.04 

Me-Bzm 

+0.08 

+0.08 

+0.12 

+0.09 

5,6-Me-Bzm 

+0.12 

+0.15 

+0.21 

+0.13 

Pro-Bzm 

+0.05 

+0.04 

+0.08 

+0.06 

HMB 

+0.03 

+0.05 

+0.07 

+0.04 

HEB 

+0.02 

+0.04 

+0.06 

+0.05 
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The AlogK values for ternary complexes containing benzimidazole or substituted 
benzimidazoles are less negative or positive than the value (-0.2) expected on the 
statistical grounds reflecting the extensive stabilization of these complexes. This 
stabilization is due to electrostatic attraction between the negatively charged oxygen’s 
of phosphate oxygen of OP / PP / TPP in the coordination sphere of the metal and 
positive charge on metal ion. 

The stability of the ternary complexes in terms of inorganic phosphate decrease in 
the order: TPP > PP > OP. The ternary complexes of TPP are more stable than PP 
due to more electrostatic attraction between TTP^' and positive charge on 
coordination sphere of the metal ion. Further TPP, though, potentially tridentate may 
be binding to the metal ion in a tridentate fashion. This is another reason that TPP 
forms more stable ternary complexes than PP. The ligand, OP has comparatively less 
negative charges L^' than PP (L'*') or TPP (L’‘), hence, it forms least stable complexes. 

The sequence of stability of ternary complexes with respect to other ligands 
decrease, in the order: 5,6 Me-Bzm > Me-Bzm > Bzm ~ Pro-Bzm > HMB > HEB. 
This trend is in accordance with basicity of ligands and steric hinderance caused by 
substituents present in position 2. 

The stability constants of ternary complexes of bivalent metal ions is in the order: 
Cu^"^ > Ni^* > Co^"^ < Zn^"^. This order is in conformity with Irving - William’s natural 
order^'* of stabilities. 

Distribution Diagrams : 

From a knowledge of proton-ligand and metal-ligand formation constants, the 
distribution of the total metal among the various metal-ligand species as a function of 
pH has been calculated by a computer program BEST*'’. As a representative case, the 
distribution diagrams for [Bzm-Cu-OP] and [PP-Ni-Bzm] systems are given in ;Fig.l 
and 2. 

[Bzm-Cu-OP] 

The maximum percentage of the Cu-Bzm or ML is 61.2 at pH 6.1, as the pH 
increases the percentage of ML reaches minimum. The other binary species [Cu-OP] 
has attained maximum percentage of 5.1 at pH 6.6. After pH 5 the formation of 
ternary complex has started (appearance of curve LMA) and its percentage goes on 
increasing with increase in pH, with simultaneous decrease in percentage of both 
binary species (i.e. MA and ML). The ternary species is the only predominant species 
in higher pH. 
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Fig. 1 : Species Distribution Diagram of Bzm-Cu-OP (1:1:1) 


[PP-Ni-Bzm] 

It is seen that binary complex ML or Ni-Bzm is formed to the extent of 69.5% at 
pH 6.9. It indicates complexation is taking place in stepwise manner. The maximum 
formation of binary species, MA is to the extent of 3.6% at pH 7.6 after this pH, the 
ternary complex is the only species present. It is concluded that up to pH 6.5 there 
exists only binary species MA (3.6%) and another binary species ML(69.5%). After 
pH 6.5 the MAL is the only species found. 



Fig. 2: Species Distribution Diagram of B 2 m-Ni-PP (1:1:1) 
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Abstract 

The oxidation kinetics of D-galactose with quinolinium chlorochromate, C 9 H 7 NHCr 03 Cl, have 
been studied in aqueous acetic acid medium 50% (v/v) in the presence of perchloric acid at 
constant ionic strength. The reaction under pseudo-first order condition, is first order with 
respect to both the [oxidant] and the [substrate]. The reaction is markedly catalysed by [H^ and 
the effect of [NaC 104 ] is negligible. Increasing the dielectric constant of medium decreases the 
reaction rate. A 1:1 stoichiometry is observed and the reaction did not induce the 
polymerisation of acrylonitrile. The reaction has been ceirried out at different temperatures and 
the variation parameters have been computed. A probable mechanism of the reaction is 
suggested. 

(Keywords : QCC/QCC-galactose intermediate/oxidation/hydride ion transfer) 

Introduction 

The kinetics of oxidation of reducing sugars which contain a large number of 
functional groups is of considerable interest from mechanistic point of view and have 
been fairly well studied^'^^. Quinolinium chlorochromate^^ adds to the select list of 
new Cr(VI) reagents introduced recently as an oxidising agent for the effective and 
selective oxidation of organic substrate under mild conditions. It has been found quite 
useful as an oxidant synthesis*^’^^ and is being increasingly used as an oxidant in 
kinetic investigation of many compounds*"^’*^. D-Galactose plays an important role in 
carbohydrate chemistry. A survey of literature, however, revealed that no work has 
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been carried out on the kinetics of oxidation of D-galactose with quinolinium 
chlorochromate. The study reported in this paper deals with the oxidation kinetics of 
D-galactose, by quinolinium chlorochromate in aqueous acetic acid 50% (v/v) 
medium with a view to throw light on the mechanistic aspects of the oxidation 
process. 


Materials and Method 

Materials : Quinolinium chlorochromate was safely prepared’^ by a careful 
addition of quinoline (Loba-chemie) to a solution of chromium trioxide in 6M HCl 
followed by the filtration of yellow orange crystals. Its purity was checked by 
estimating Cr(VI) iodometrically. Solution of D-galactose was always freshly 
prepared in aqueous acetic acid 50% (v/v). The ionic strength was maintained constant 
with the use of concentrated solution of NaC104 (B.D.H). Perchloric acid (E-Merck) 
and all other chemicals were used as such without further purification. 

Kinetic measurements : The reaction were performed under pseudo-first order 
conditions by keeping a large excess the of the [D-galactose] with respect to [QCC]. 
The medium of the reaction was always 1 : 1 (v/v) acetic acid-water in the presence of 
perchloric acid. Kinetic measurements were made on a Systronics 1 06 
Spectrophotometer at 440 nm. The optical density was measured at various intervals 
of time. Computations of rate constants were made from the plot of log [QCC] against 
time. 

Stoichiometry and product analysis : The stoichiometry of the reaction was 
performed by conducting the oxidation of D-galactose under the conditions of a 
known excess of quinolinium chlorochromate. A mixture of D-galactose (0.010 mol 
dm'^) and QCC (0.10 mol dm’^) was kept for several hours at 30°C for the reaction to 
go to completion. The unconsumed oxidant was estimated iodometrically at the end of 
the reaction. The stoichiometry was found to be 1 : 1 consistent with the following 
equation. 

C6Hi206+C9H7NHCr03Cl — C5Hio05+HCOOH+C9H7NHCl+Cr02 

Arabinose Formic acid 

For product analysis a reaction mixture containing D-galactose and quinolinium 

chlorochromate in the stoichiometric proportion of 1 ; 1 was left to equilibrate at 30°C 
for 24 hours. The reaction mixture was neutralised with NaHCOs, extracted with 
chloroform, washed with water and dried over anhydrous MgS04. 
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The formation of D-arabinose was confirmed by phenylhydrazone formation'®. 
The presence of formic acid was confirmed by spot test*’. 

Results and Discussion 

Empirical rate law : The pseudo-first order rate constants were determined at 
various initial concentrations of reactants. The results obtained are given in Table 1 . 

Table 1- Rate constants for oxidation of D-galactose by QCC at 30° C. Solvent : acetic acid-water (50- 
50% v/v). 


[QCClxlO^ 

(mol dm'^) 

[D-GaIactose]xl0^ 
(mol dm‘^) 

[H^]xl0 
(moi dm'^) 

kix lO" 

(sec"') 

(mol''dm’sec'‘) 

0.66 

1.66 

6.13 

6.83 

- 

1.33 

1.66 

6.13 

6.82 

- 

2.00 

1.66 

6.13 

6.84 

- 

2.66 

1.66 

6.13 

6.82 

- 

3.33 

1.66 

6.13 

6.82 

• - 

3.33 

0.33 

6.13 

1.36 

4.12 

3.33 

0.66 

6.13 

2.70 

4.09 

3.33 

1.00 

6.13 

4.12 

4.12 

3.33 

1.33 

6.13 

5.43 

4.09 

3.33 

1.66 

6.13 

6.82 

4.11 

3.33 

1.66 

3.1 

3.38 


3.33 

1.66 

6.1 

6.82 

- 

3.33 

1.66 

7.7 

8.47 

- 

3.33 

1.66 

9.2 

10,22 

- 

3.33 

1.66 

10.7 

11.92 

- 


Plots for different concentration QCC versus time were linear and the rate 
constants were independent of initial concentration of QCC, showing first order 
dependence of the rate on [QCC]. The reaction is first order with respect to [D- 
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galactose] also. A plot of log k, against log [D-galactose] was linear with a slope of 
unity thereby confirming first order dependence in [D-galactose] (Fig. 1). 



Fig. l-Plot log ki versus log [D-glactose] [QCC]=3.33xl0‘’ mol dm’’; 
[HCIO4]=6.13 x 10-‘ mol dm’’; iNaC 104 ]= 1 . 66 xlO ^ 

Rates of oxidation were found to increase with increase in [H^] and the slopes of 
the plots of log ki versus log [HCIO4] was approximately unity showing that reaction 
is acid catalysed and follows a first order dependence in [HCIO4] (Fig. 2). 



1+log [HCIOJ 

Fig. 2-Plot log ki versus log [tf] [QCC]=3.33xlO’^ mol dm"^; 
[D-galactose]=l .66x10’* mol dm’* [NaC 104 ]= 1 . 66 xlO’'; mol dm’*; 
Solvent ; acetic acid-water (50-50% (v/v/); temperature = 30°C. 
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Consequently the empirical rate law is described as follows : 

_ <^tQCC] ^ jD . Galactose][HC104 ] 
dt 

Effect of ionic strength : The reaction rate was not influenced by ionic strength 
when NaC104 was initially added to the reaction mixture over the range of 0.66xl0‘’ 
to 4.00x1 0'*mol dm'^. Similar observations were also reported in the oxidation of diols 
by PCC‘* and chromic acid*’. 

Effect of radical-forming agent : When the reaction was initiated by adding a 
radical scavanger, acrylonitrile, into a solution containing D-galactose and QCC, no 
retardation in the rate was observed. .No turbidity due to the polymerisation of 
acrylonitrile was observed. Thus the formation of radical intermediate is ruled out in 
the course of the reaction. 

Effect of solvent : The reaction has been studied under varying composition of 
acetic acid-water mixture. It has been observed that the reaction rate increases with 
the increase of solvent (CH3COOH) in acetic acid-water mixture. (Table 2 and Fig. 3). 

A linear plot between log kj and 1/D (inverse of dielectric constant)^® with a 
positive slope suggesting an interaction between an ion and a dipole^*. 



Fig. 3-Plot log kj versus log 1/D [QCC]=3.33xl0'^ mol dm 
[D-galactose]= 1.66x10'^ mol dm'^; [HC104]=^. 13x10''; mol dm 
[NaC104]=1.66xlO'' mol dm'^; temperature=25“C 
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Table 2- Dependence of rate on solvent composition 

IQCC] = 3.33 X 10'^ mol dm'^ [D-galactose] = 1.66 x 10'^ mol dm'^ [HCIO4] = 6. 13 x 10'' mol 
dm'^ [NaC 104 ] = 1.66 x 10'' mol dm'^ [NaClOJ = 1.66 x lO"' mol dm‘^ Temperature - 25‘'C. 


CH3COOH 

: H 2 O (v/v) 

1/D 

ki X 10"* (sec'^) 

20 

80 

0.016 

3.28 

30 

70 

0.018 

3.58 

40 

60 

0.020 

4.02 

50 

50 

0.024 

4.84 

60 

40 

0.028 

5.82 

70 

30 

0.036 

8.23 


Effect of temperature : The reaction rates at different temperatures were 
determined and the values of activation parameters were calculated from the slope of 
linear plot of Ig ki versus T''. The data are presented in Table 3. A perusal of data 
(Table 3) shows that these reactions are characterised by a high negative value of 
entropy of activation (AS*). It indicates that solvation effects are predominant in the 
reaction, which suggests the formation of a charged and rigid intermediate state. 
Further, the high positive value of energy of activation and enthalpy of activation 
indicate that intermediate is highly solvated. 

Table 3- Temperature dependence and activation parameters of the oxidation of D-galactose by QCC. 

[D-galactose] = 1.66 x lO'^ mol dm'’; [QCC] = 3.33 x 10'’ mol dm'^; [HCIO4] = 6.13 x 10''; 
[NaC 104 ] = 1.66 X lO’' mol dm"’; 

Solvent : acetic acid - water (50-50% v/v) 


Temperature 

(“C) 

kixlO" 

(sec"') 

Ea 

kJ moP^ 

AH' 

kJ mo!'^ 


-AS* 

JK-' mol ' 

20 

3.28 

- 

53.17 

91.28 

130.06 

25 

4.84 

54.26 

53.13 

91.31 

130.14 

35 

9.88 

58.10 

53.05 

93.25 

130.53 

40 

14.29 

- 

53.01 

93.85 

130.47 


Aldoses in aqueous acidic solution exist as equilibrium mixture of a - and |3- 
pyranose (Scheme 1) . The equilibrium mixture of D-galactose mainly consist of p- 
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anomer with equitorial orientation of the glycosidic hydroxyl group. The p-anomer is 
supposed to be more reactive species^^. 



H OH H H 

a - Galactose, 30% P - Galactose, 70% 


Scheme 1 

An increase in the oxidation rate with acidity suggests protonation of QCC 
species. This protonated QCC interacts with the substrate in the rate determining step. 
Involvement of such species is well known in chromic acid oxidation^^. 

The oxidation of D-galactose may be assumed to follow the direct hydride ion 
transfer of the axial anomeric proton of p-galactose as shown in Scheme 2. 

However, UV spectra recorded did not show the existence of intermediate 
complex formation. It indicates the instability of intermediate (Fig. 4). 



Wavelength, nm 

Fig. 4-Electronic spectra of QCC in aqueous acetic acid 50% (v/v) (A) and QCC 
in aqueous acetic acid 50% (v/v) with D-galactose. 
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Scheme 2 


Rate Law- 


K, 


K, 


QCC + 


QCCH^ 


( 1 ) 
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QCCH^ + S ■' - -■J - > I + Cr(IV) (2) 

I + QCC — > Products (3) 

Rate = -^t^^^SJa[QCCH"][S] = k,[QCCH"][S] 
dt 

Rate = k, [QCCH^] [S] (4) 

Applying law of mass action to equation (1) we get, 

K, 

[QCC][H^] 

[QCC H"] = K, [QCC] [H"] (5) 

substituting the value of [QCCHf"] from equation (5) into equation (4) 

Rate = k,K, [QCC] [H^] [S] 

= kobs [QCC] [H^] [S] 

Conclusion 

The oxidation of D-galactose by QCC in aqueous acetic acid medium proceeds 
via hydride ion transfer mechanism involving specific cleavage of C 1 -C 2 bond of the 
substrate to give the products. 
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Abstract 

The spce J is defined as a generalization of space \ Its properties including completeness and 

inclusion are studied. It is shown that the space is dense in G^’ The generalized Hankel Potential 
is defined and its properties are studied. The L^-space of all such Hankel potentials, denoted by 
, is defined. It is proved that the space is a Banach space. It is also proved that the generalized 
Hankel potential is an isometry of onto 

(Keyvi^ords ; Sobolev type space/generalized Hankel potential/Beurling type ultradistribution/ 
Pseudo-differential operator/Hankel transform). 

Introduction 

Let CO be a continuous real valued function defined on / = (0, <») satisfying the following 
conditions : 

(a) 0 < C0(^ + Tl) < t0(4) + Cl)(Tl) V t T1€ / 

(P) f 

0 

Sometimes we assume the following condition also : 

(Y) ©(^) > a + log (1 + ^), ^ e /, for some real number a, and b>0. 

The class of all (O satisfying (a) - (Y) will be denoted by M. 

Clearly, (o(^) = log (1 + 4), ^ > 0 and a)(^) = ^“ ^ > 0, 0 < a < 1, 
satisfy (a) - (y). 


From (y) it follows that 
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^ < exp {-alb + o )(^) ! b),\> 0. 

Let {h,. u) (£) = ui^) denote the Hankel transform of m e L (0, «>) defined by 

M (^) = {hf, m) (^) = ^ 

Then the space f^iF) is defined to be the set of all complex-valued infinitely 
differentiable functions ^onl such that 

^ (^) = sup I (x:"‘ didxf (j) (x)) | < «> 

’ xel 

and 

pj:, ($) = sup|e^^(r‘d/4)*r'^-''^$(^))i<oo. 

Clearly, = ^ - For (o(a:) = log(l+x), 0^(1) = the Zemanian - 

spaced The dual of well be denoted by (H^Y which consists of all Hankel transformable 
Beurling type ultradistributions. 

We shall make use of the function space (0, oo), i <p < oo, which consists of all 
those real measurable functions on (0, co) for which 

Up 

< 

where 

daix) = (2^‘r(p+l))-‘ x^t^ + ’nh:. 

A tempered distribution u belongs to the Sobolev space H' (R"), se R, if its Fourier 
transform corresponds to a locally square integrable function such that 


\\f\\p = \f^\mrdo(x) 
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It has been generalized to by Pahk and Kang^ using ultradistribution theory of 
Beurling and Bjorck. It has been further generalized to 1 < p < <», by Pathak^. The 
Sobolev type space R , \<p<°° was introduced by Pathak and Pandey'^ using 

Hankel transformations. 

In this paper a space is defined as a generalization of the space It is shown 
that for 5 e /f, p. > -1/2, 1 <p < «>, the space is complete. Moreover, the generalized 
Hankel Potential is defined and its properties are studied. The L^space of all such 
Hankel potentials, denoted by is defined. It is proved that the space ^ is a Banach 
space. It is also shown that the generalized Hankel potential is an isometry of 
onto KT- 

The Generalized Sobolev type space 

Definition 1 : For s e R, p e R , 1 Sp < », the generalized Sobolev type space 
is the set of all ultradistributions u e such that its Hankel transform corresponds 
to a locally integrable function over I = (0, °°) for which 

For = log (1 + 4), the space G^’ ^ reduces to the space G^’ 

Theorem 1 : The space Gfj 1 < p < «>, j > 0, is complete. 

Proof. Let {m^} be a Cauchy sequence in Then (^)} is a Cauchy 

sequence in D’. Since is complete, therefore there exists a function/ 6 2/ 

such that 



f IF z&k ^ «>. 

(1) 

Set 

= ^>0 

(2) 


Clearly, g& e (7^)'. Therefore its inverse Hankel Pansform exists. Now (2) gives 
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Therefore, by Definition I, ge. 

Let g = u. Then 

(h^ u) (^) = {h^ (A-' (5) = g{^) =A^) 


Now it follows from (1) that 


Uj, U in 
This proves the theorem. 

Theorem 2 : For t > r, c and the inclusion map is continuous. 
Proof. For u e and r > r, we have r = r + 8, 8 > 0. 

Now, 


Since 


r i ih^ « (^) f 4? < f r I (h, u (0 f dl 


we have 


II » lie;; £ II “lb;; v»ecs,; 

Hence c;;sg;-;. 

{m^.} -> 0 in ^ then 

{Uf.} -> 0 in . 


(3) 


Again let 
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This proves the continuity of inclusion. 

Theorem 3 : The space is dense in 

Proof. Let u € Then u) (^) e L^. Since D is dense in L’’, there exists a 

sequence {(j)^ ) of functions belonging to D such that 

{h^u){^) in L^ 

Define \\ij (x) = /i'* (j)^.)(.t) 

Since (J)^- € D, therefore D-,D being a subspace of (j)^- also belongs 

to Hence its inverse Hankel transform exists; so that >j/y e H^. Now, 

II - u llc^.* = II (h^ (V; - «)) il) \\p 

= II <1>; - m) (0 ip -> 0 as J -4 oo ; 

so that 


y,. -> M in 

The following lemma will be used in the sequel. 

Lemma 1 : For tliree real numbers , s and S 2 satisfying jj < 5 < J 2 Ve > 0 there exists 
C (e) > 0 such that 

(4) 


Proof. For p = and C (e) = £ 


as in Chazaran and Piriou^, we have 




< ep^2 + e 





80 


R.S. PATHAK and K.K. SHRESTHA 


i.e. 1 < B + B ^ (5) 

1 

Put k = B*^^ Then, by the inequality (5) we get 
1 < + (kpfr^ 


If /:p > 1, then we have (/:p /2 ^ > i. Hence the inequality is satisfied. 

If /:p < 1, then kp^~^ > 1 so that the inequality is true in this case also 
This proves the Lemma. 

Theorem 4 : Let s and ^2 be three real numbers satisfying Si<s< ^ 2 . 
Then Vb > 0, there is a positive number C(£) such that 

II « lie"'" ^ e II « llc/’-’i + C (8) II M ||(/^, Vm e 

Proof. If w G then by Theorem 2, we have 
ue GP-^andHe 

Now, 




By Lemma 1, Ve > 0, there exists a C (e) > 0 such that for all p > 1 we have 


I “ ll(f-^ ^ e 


I (h^ u) (^) I T + C (^) fj^ I (h^ u) & r 4 


so that 
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ii « ^ e II « ll(/''2 + Qe) II « llcf’^i 

0341 Oi,n (a,n 

Theorem 5 : For s <0, (1) c 

Proof. Let ue (I), Then from Macauley-Owen^, we have 

II V w \y = Ii u \y . 

Now, 

for some constant K; so that u e 

For s=0, p=2, the space G^’° is identical with (I). 

Example 1 : For ae R, let be defined by 

<8„,<t>> = (l)(a) V<t)G M. 

Then 

8a e Ga/n for 5 < - 1 / (pb). 

Indeed, we have 

so that 



II (\ 5.) ® It = II ^ J, m II, 
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^ „ M/p 

V ) 

Since | V|a (^) j is bounded by the last integral is bounded by 

- *' * ^ by property (y) 

S (1 + ^r* 4 

The ^-integral is convergent for 5 < - 1 / (pi,). 

The Generalized Hankel potential 
Let a(x) 9t 0 be a multiplier in such that 

V ^ (0* ®°) V m = 0, 1, 2 (6) 

Then, for f € we define the products/a"* and/a"'" 6 by the following relations. 

</a'", 0 >=</, a™ 0 > V 0 € 
and 

</a“'", o'" 0 >=</, 0 > V 0 £ 

where m is a non-negative integer. 

Now, the generalized Hankel potential is defined by 

H™ = /j~‘ (fl-™ (A^ u)) V « 6 (i^)' 


(9) 
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Since the generalized Hankel transformation /i^' is a continuous linear mapping of 
onto the same being true for we conclude that the generalized Hankel potential 

defined by (9) is a continuous linear mapping of onto for all m 6 Z. Also 

(H; u) W = J (xci) (^) u) (^) 

This shows that the generalized Hankel potential is a pseudo-differential operator with 
symbol a~”'. For a(x) = (1 + he R, (9) reduces to 

(H;, u) ix} = (^(1 + (A^„) (^) for 5 = mA 

which is the Hankel potential investigated by Pathak and Pandey**. 

Theorem 6 : Let « e (H^'. Then for m,leZ'we have 

(i) H” Hi, « = HjT' « 

(ii) Hj « = «. 

Proof. For u e we have 

H“ H;, H = A;;‘(a-”A^(Hi,«)) 

» A-‘(a-" A^(A-‘(a-' /^i,«)) 

- A"‘(a~<'"-^A^«) 

= 

This proves the first part. For the second part, we have 

M = Aj^* (a® Aj, «) = u. 
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.THE SPACE 

For m e Z and 1 </?<«, the space (0, »=) is the set of all those ultradistributions 
u e (Ft^Y for which 

The norm of this space is given by 

Theorem 7 : The space (0, <») is a Banach space with respect to the norm 

II llnj7,co,n • 

Proof. Let {m^.} be a Cauchy sequence in (0, °°). Then 

is a Cauchy sequence in L^. Since is complete, there exists a function 
« e such that 

jj-m M in as — > o« . (10) 

Let V = u. Then 

V ^ ^n+i/2 „ or H;;'" V = «. 

Hence v e . 

From (10) it follows that v in as ^ > «» . 

Theorem 8 : The generalized Hankel potenial is an isometry of 



GENERALIZED SOBOLEV TYPE SPACES 
(0, °°) onto (0, «); and we have 

II ^ llm+<,p,(0,|i “ II *t* llm,p,(0,n 

Proof. For (j) G (0, oo), we have 

“ II *1^ llm^,to,n' 

Again let v € . Then 

so that H;;' 4) e and 

H[,H-'v = v. 

That is, for each v in there exists an v G such that 

V = V. Hence, is onto. 

Theorem 9 : For 1 <p < <», m = 0, 1, 2, 

I! H™ (t) II < c II lip V G . 


In the sequel we shall need the following lemmas. 
Lemma 2 ; If /, g G (0, <»), then 
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(f#gf(x)=fix)g(x), 0<x<o°, ( 11 ) 

where 

hx) = fj{xt)nt)da(t) ( 12 ) 

and 

j{x) = 2^‘r(n+l)^-'‘7j,(x). 

The poof of this lemma is given by Haimo^. 

Lemma 3 : L&tf^x) € L}^ (0, «>) and g(jt) e (0, «>), l<p < » and let 

f* ix, y) = Az) D{x, y, z) dc (z). (13) 

Then the Hankel convolution 

(f# g) (x) = {X, y) g(y) dc(y) (14) 

satisfies the norm inequality 

ll(r#g)Wlip^lUA:)||ilUWllp. (15) 

The proof is given by Pathak and Pandey^. 

Proof of Theorem 9. Proceeding as in the proof of the Theorem 4.3 of Pathakand Pandey^ 
we can show that 

[v «)) (^) = K Q . ftp g(t)) (^) . ( 16 ) 

Now assume that 


t>(0 = 


(17) 
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a-"'& = 

Then 

M = h-' a-" (^)) (/) (18) 

^ ,-^-l/2 r ^ J 2^^ r (|X + 1) ^"2^*"' do (4) 

0 ^ 

Using (13) we obtain 

M = «■'" (^) 2^^ r (II + 1) (^)-^‘ at) da a) 

= «■'" (%)j (^) da a) = (a-”) (0 e L* (0, =o) by (6). 

Now, using (16), (17) and (18) we obtain 

(f# g)) a) = V Am (^))] Am ^ 

= a-“ (|)(/»^,<t))(^); 

so that 

(T# g) = A-‘ (a-'" h^<l>) = H;<I>. 

f#g = H” (1>. 


Therefore, 
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Hence, 

Using Lemma 3 we get 

ll,---'" h; 4,11, S ll/ll, Hr-'-''"* II,; 

SO that 

This completes the proof of the theorem. 

Theorem 10 : For 1 < p < <» and m < I, 

Vi/'’’ c VT'’’ 

and we have 

li<t>iU<0,(X ^ C lid) 

Proof. Let d) e Then 

■ H-' d)6 


Also 

H-'"d> = d). 

Now, 

II<1>1U»., = 11 H-d) II, 


= ||y-M/2H^«H-'d)||, 


< C||y-M/2H-'d,||^ 
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Theorem 11 : For 1 <p < «« and any symbol o'" , me Z, define the pseudo-differential 
operator 

(A (D) u) (x) = h~^ [n'" u] (x) . 

Then 

A(D) : is a bounded linear operator. 

Proof. Consider the following operators ; 

-> 

A (D) h;;' : ^ 

Tj/-m . ^/0,p . 

* *^.0),^ ^ 0),[l 

The operator is bounded. Indeed, for (j) G we have 

II ♦11;,.., = II II, < ~ 

and 

l|H-'(>|lo.p.a>., = <oo. 

The second operator is also bounded. For, let (j) G . Then 

Now, 


II A (D) h; 4, Ik,,.,, = II J-'-'” .1 (0) h;: ♦ II, 
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For the third operator we see that 

= II II, <-; 

SO that third operator is also bounded. Hence the product A (D) is a bounded 
linear operator forna into 

u/'P ► 

co,^ 

H m-l 

n 

WZ-^ A(D) 

Moreover, by Theorem 8, the operators and in Fig. 1 are isometric and onto. 
A generalization of the above theorem is given below. 

Theorem 12 : For Z and symbol a^\ the pseudo-differential operator 
A{D) = h~J a'” h,, : 

is a bounded linear operator. 

Proof. Let M € Then 

iiHU^ = ii3'''^"''^H-;ii<oo. 
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II ^ (D) u = 11 H;r' MD) u \\p . 

Here A(D) is a pseudo-differential operator with symbol a‘. For, 
AiD)u = h-^[a'^^‘hJ,A{D)u]\ 

= ;i-‘ [a-"' (a- V»)]] 

Then (19) gives 


( 19 ) 


» <0) “ ' II V f"' 'V “ II, 

= II Vii, < ~. 

This completes the proof of the theorem. 
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Abstract 


In this paper, the thermosolutal instability of Rivlin-Ericksen visco-elastic fluid in the 
presence of rotation in porous medium has been considered to include the effects of magnetic 
filed and variable gravity field. For the stationary convection, the stable solute gradient has a 
stabilizing effect, whereas rotation has stabilizing/destabilizing effects when gravity 
increases/decreases upwards. The magnetic field has stabilizing effect in the absence of 
rotation, whereas in the presence of rotation, magnetic field has stabilizing effect when 


< 1 + X H ^ I [(l 4- ] * The mediiun permeability has destabilizing effect in the 

absence of rotation, where as in the presence of rotation medium permeability has stabilizing 


effect when T, 

A 


r 


<e 


PQ: 


l + X + ^ 
\ e 




increases upwards (i.e. X > 0). The case 


of overstability is also considered wherein the sufficient conditions for the existence of 
overstability are derived. The effect of stable solute gradient, rotation, ma^etic field and 
medium permeability have also been shown graphically. 


(Keywords : thermosolutal instability/Rivlin-Ericksen elastico-viscous fluid/rotation/magnetic 
fieid/variable gravity.) 


Introduction 

The theoretical and experimental results of the onset of thermal Instability 
(Benard convection) in a fluid layer under varying assumptions of hydrodynamics and 
hydromagnetics has been treated in detail by Chandrasekhar'. Veronis^ has 
investigated the problem of thermohaline convection in a layer of fluid heated from 
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below and subject to a stable solute gradient whereas the problem of thermohaline 
convection in a horizontal layer of viscous fluid heated from below and salted from 
above has been studied by Nield^ 

With the growing importance of non-Newtonian fluids in modem technology and 
industries, the investigations on such fluids are desirable. There are many viscoelastic 
fluids that cannot be characterized by Maxwell's constitutive relations or Oldroyd’s 
constitutive relations. One such class of viscoelastic fluids is the R.ivlin-Ericksen fluid. 
Rivlin and Ericksen^ have proposed a theoretical model for such viscoelastic fluid. 
This and other class of polymers is used in the manufacture of parts of space-crafo, 
aeroplanes, tyres, belt conveyers, ropes cushions, seats, foams, plastic, engineering 
equipments etc. Recently, polymers are also used in agriculture, communication 
appliances and in biomedical applications. When fluid permeates a porous material, 
the gross effect is represented by the Darcy's law. As a result of this macroscopic law, 
the usual viscous terms in the equations of Rivlin-Ericksen viscoelastic fluid motion is 


replaced by 


— ^(p + p')— ^ 

Jt, ^5/ 


, where p and p' are the viscosity and viscoelasticity of 


the thermosolutal instability of Rivlin-Ericksen fluid, k\ is the medium permeability 
and q is the Darcian (filer) velocity of the fluid. 

The problem of thermosolutal convection of rotating fluid in a porous medium is 
of importance in geophysics, soil sciences, ground-water hydrology and astrophysics. 
The development of geothermal power resources hold increased general interest in the 
study of the properties of convection in porous media. The scientific importance of the 
field has also increased because hydrothermal circulation is the dominant heat transfer 
mechanism in the development of yoimg oceanic crust Lister^. Srivastava and Singh 
have studied the unsteady flow of a dusty elastico-viscous Rivlin-Ericksen fluid in the 
presence of time-dependent pressure gradient. Recently, Sharma and Rana have 
studied the instability of streaming Rivlin-Ericksen fluid in porous medium in 


hydromagnetics. 

Keeping in mind the importance in ground-water hydrology, soil sciences, 
geophysics and astrophysics, the effect of magnetic field on thermosolutal convection 
in ^Wvlin-Ericksen viscoelastic rotating fluid in porous medium in the presence of 
variable gravity field has been considered in this paper. 


Formulation of the problem and perturbation equations : 

Consider an infinite, horizontal, incompressible Rivlin-Ericksen viscoelastic fluid 
of thickness d, bounded by the planes 2 = 0 and 2 = in an isotropic and homogenous 
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porous medium of porosity e and medium permeability k\, which is acted upon by a 
uniform rotation 0(0,0,0), magnetic field H{0,0,H) and variable gravity g(0,0,-g). 
The layer is heated and soluted from below such that a uniform temperature gradient 



dz j 


and uniform solute gradient 


P' 


dz j 


are maintained. 


Let p, p, T, C, a, a', q(u, v, w) u, o' and Pe denote, respectively, pressure, 
density, temperature, solute concentration, thermal coefficient of expansion, an 
analogous solvent coefficient of expansion, fluid velocity (initially zero), kinematic 
viscosity, kinematic viscoelasticity and magnetic permeability. Then equations of 
motion, continuity, heat conduction, solute concentration through porous medium and 
equation of state of Rivlin-Ericksen fluid are 


]_ 

e 


dt 6^ ^ 


'-1 

,P«; 


V/j+lf 1+^ 



4^Po^ ’ ’ 


V.q = 0, 


( 1 ) 

( 2 ) 


£^ + (f.V)r = KVT, 


( 3 ) 


E'— + {IV)C = kV^C, (4) 

and p = po - To) +a'(C - Co)] (5) 

where the suffix zero refers to values at the reference level z = 0. The kinematic 
viscosity u, kinematic viscoelasticity u', the thermal diffiisivity k and solute difiusivity 
k' are all assumed to be constant E=e+(l-6)(pjC^/PoC^) is constant and is a 

constant analogous to E, but corresponding to solute rather than heat, Pf, c^, po and Cf 
denote the density and heat capacity of solid (porous) matrix and fluid respectively. 

Maxwell's equations and the modified Ohm's law yields 
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e — = Vxfex.^)+eTiV^^ 
dt ^ ' 

(6) 

V.F=0, 

(7) 


where t) is the electrical resistivity. 

The steady state solution is 

q =(0,0,0), r=-p 2 +ro,C = -P'2 + Co,p = po(l+apz-a’p'z). (8) 


Consider a sniall perturbation on the steady state solution and let 5p, 5/>, 0 , y 
q(u, V, w) and h{K, hy, denote, respectively, the perturbations in density p 
pressure p, temperature T, solute concentration C, velocity q = (0,0,0) and magnetic 
field H = (0,0,H). 

The change in density 5p caused by perturbation 0 and y in temperature and solute 
concentration, is given by 

5p = -po(a0-aV). (9) 

Then the linearized perturbation equations relevant to the problem are 


1 dq 
Tdt 



( 10 ) 


V.?=0, 


( 11 ) 


= Ph’+kV^0, 
dt 


( 12 ) 



P'w+K'V^y, 


(13) 
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e — = Vx(^Xj^)+e r| , 


V.h=0, 

In the Cartesian form, equations (10) - (15) with the 
expressed as 


(14) 

(15) 

help of equation (9) can be 


1 _ 1 dhp 1 

€. dt Po Sx: A:, 

1 _ 1 dbp 1 ^ 

€dt p„ dy k, 


r 


u + u 






o + u 


dt 


u + 


47Tpo 


dh^ dh 
dz dx 


\ 


+ -Qv, 
€ 


dh^ 




dt) 47tp„ 


dz dy 


■Qu , 


\ dw _ 1 ddp 1 

e dt Po dz 


r 


u + u 


V 


dt 


w + g(a0 - a'y), 


du dv dw . 

— + — + — = 0 , 

dx dy dz 

dt dz 


dh dv 2 

e — - = H — +6 TiV A , 
dt dz ' ' 


dh dw / 

dt dz 


— K + — h^-¥ — K =0 

dx ^ dy ^ dz ' 


E — = Pw + kV^0, 
dt 


(16) 

(17) 

(18) 

(19) 

( 20 ) 

( 21 ) 

( 22 ) 

(23) 

(24) 
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£'^ = pV+KV^y, (25) 

dt 

Operating equations (16) and (17) by — and respectively, adding and using 

cbc oy 

equation (19), we get 


i = _Lfv^ - — ] 5;) - -f o + u'— 1 w + - 20^ (26) 

p„v Szv 47tp„ 


, 3 3 3 

where = — + — + — , 

dx By dz 


and ^ ^ is the z-component of vorticity. 


Operating equations (18) and (26) by 


5 1 d 

and — respectively and adding, to eliminate 5p between 

dz J dz 


equations (18) and (26), we get 


1 3 /„j \ (3^ 5 1/ Q f \ if < ^ /r 72 7 t 

(V •H')=g — - + — - (a0-ay) u + u — V w+2-s (V /z J-2Q— 

edt^ \^3x dy J k\ dt) 47Tp„ 3/'' ' dz 

(27) 


Operating equations (16) and (17) by 

d d 

and — respectively and adding and using equations (19), we get 
dy dx 


13C If ,3^ di, ^^dw 
e 3t k\ dt) 4irp„ dt dz 
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dh dh 

where ^ = — is the z-component of current density. 

ox ay 

Operating equations (20) and (21) by 

3 3 

-—and — respectively and adding, we get 
oy ox 


^ = /f-^+GTlV^4, 

edt dt 

The dispersion relation : 


(29) 


Analyzing the disturbances into normal modes, we assume that the perturbation 
quantities are of the form 


[w, e,y, K C, = [W{z), ©(z), r(z), Z(r), X(z)] exp [ik^ + ikyy + nt] 
where k^, ky are the wave numbers along the x andy directions respectively, 

* = fc+^,v)^ 

is the resultant wave number and n is the growth rate. 

Using equation (30), equations (27) - (29), (22), (24) and (25) become 


"z = 


^ A 

ydz^ 


|jT = -gA^^(a0 - aT) - -^(u + - k 


+ 


1^. d 


”z = -l 


f 




o + u — 
dt 


Z + 


4np, 


4np^ dz 
■DX + 2Q. 


/ 


dz^ 


K-2Q. 


W 

dz 




dz 


„ „dZ I 
enX = H — +en| 
dz 


f ^2 




~-k 
dz ) 


(30) 


(31) 

(32) 

(33) 
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enK = H +en — --k K 

dz \dz J 


€n© = pfr + K -rj-k © 
[dz 


( d^ \ 

enr = ^'W + K' T 

[dz J 


Expressing the coordinates (x,y,z) in the new unit of length d and letting a = kd. 


G = , F = -^ and D' =d— and the superscripts * is suppressed. 

V) d dz 


Equations (31) - (36), in non-dimensional form, become 


- + — (l + aF) (£>^ © 

e ’ 47tp„u^ ^ u 

-^^r + ^Z)Z = 0, 

U 60 


- + — (I + gF) z = dx + ^^DW , 

€ /), 47tp„o e u 


[d^ -FjSig]© = - W 




[d^ -a^ - P,g]x = -\ — \dZ . 



thermosolutal instability of rivlin-ericksen 


101 








DW, 


( 42 ) 


where d , = — , is the thermal Prandtl number, p = — is the magnetic Prandtl number, 
K r| 

o ^ 

q = — is the Schmidt number, and p, = — 7 is the dimensionless medium 

k' d 

permeability. 

Applying the operator {D^ - cP' - pjo) to equation (38), to eliminate Y between 
equations (38) and (41), we get 


— + — (1 + aF)l (D^ - a^p^o) + — D^ 
e A J e 


Z = ^ {d^ -a^- pA DW. (43) 

€ U 


Eliminating K, ©, T and Z between equations (37) - (43), we obtain 


a 1 


+ — (l + aF)|(i)' -a^ - p^aJ-h^D^ 


X |— + — (l + aF)| 


(d^ -a^- Ep,o) {D^ -a^- E'qa) { d^ -a^- p,o) { d^ - a^) + ^ 
{d^ - a^ -p^o) + Sa"X{D^ -a" -Ep,a){D^ -a^ - p^o)]w 


I± 


D^{d^ -a^ -P,g){D^ -a^ -E'q^)[D^ -a^ -E'q(5)]W = Q (44) 


where Q = 


LI I 

, is the Chandrasekhar number, 

^47rpoUiij 


R-. 


J g^M. 




OK ) 


, is the thermal Rayleigh number. 
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c ~ 

goOC'P'^ 

u ■— 

. UK' 

nr 

(2Qd^\ 

T, = 

.. 


, is the analogous solute Rayleigh number, 


\ J 


, is the Taylor number. 


Here we assume that the temperatures and concentrations at the boundaries are 
kept fixed, the fluid layer is confined between two boundaries and the adjoining 
medium is electrically non-conducting. The boundary conditions appropriate to the 
problem, using equation (16) are 

W = D^W^X=DZ = @ = T = 0 atz = 0andl, (45) 

and the components of h are continuous. Since the components of the magnetic field 
are continuous and the tangential components are zero outside the fluid, we have 


DK=0, 


on the boundaries. Using the boundary conditions (45) and (46), we can show that all 
the even order derivatives of W must vanish on the boundaries and hence, the proper 
solutions of equation (44) characterizing the lowest mode is 


W= iT'osinj:z_ (47) 

where Wo is a constant. Substituting (47) in (44) and letting x = — , 

n 

R S O cs T 

■^=—5',=—, Q=- 7 ,/a=-r, P = 7C^P, andr^ =-Yand keeping in mind that o 
HZ n n ' K 

can be complex, we obtain the dispersion relation 


R. 


^16 Pi e ?a \ (\ + x + pJg^) J 


+ 5 + 

' (l+x + E'qia,) 


(l + x + Pj/a, ) (l X + Epja^ ) 


h. (l-i-/a,7t^F)l / , Q, 


^x + iEp,a, 


( 48 ) 
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Equation (48) is the required dispersion relation studying the effects of magnetic 
field, stable solute gradient and rotation on the thermosolutal instability of Rivlin- 
Ericksen fluid in porous medium. 

The stationary convection : 


When the instability sets in as stationary convection, the marginal state will be 
characterized by a = 0, the dispersion relation (48) reduces to 




1 + x 


he 


l + x , + Q 


1 + X + 


V 


IQi 
€ J 


+ s, , 


(49) 


which expresses the modified Rayleigh number Ri as a function of the 
dimensionless wave number x and the parameters Si, Qi, and P. Rivlin-Ericksen 

viscoelastic fluid behave like an ordinary Nevrtonian fluid since viscoelastic 
parameter F. vanishes with a. 

To study the effects of stable solute gradient, rotation, magnetic filed and medium 

permeability, we examine the behaviour of .^^and^^ analytically. 

dS/dT/dQ, dP 

Equation (49) yields. 


dS, 


= 1 , 


(50) 


which implies that the stable solute gradient has a stabilizing effect on the 
thermosolutal convection. The adverse solute gradient has destabilizing effect on the 

system since then becomes negative. 

dS, 
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Also from equation (49), we have 




{l + xfP 


hce^ 


1 + a: + 



l e J 


(51) 


which implies that rotation stabilizes the system when gravity is increasing upwards 
X > 0). 


dQ, he 


e 






PQ. 


1 + JC + 

V e ; 


(52) 


which implies that magnetic field stabilizes/destabilizes the system and gravity 
increases upwards (i.e., A, > 0) when 







(l + jc)P^ 

Also from (49), we have 


dPi _ 
dP ~ 


(l + x)P^ 


Xx 


(l + x) 


f 


l + ;c + 


PQ. 


(53) 


The medium permeability stabilizes/destabilizes the system when 


{Ux) 


{Ux) 




( PO^ 

1 + X+ 

V e J 

' / ' 

\ ^ J 


and X. > 0. In the absence 
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of rotation and for constant gravity field is always negative implying thereby the 

dP 

destabilizing effect of medium permeability. 


The dispersion relation (49) is analyzed numerically. Graphs have been plotted by 
giving some numerical values to the parameters, to depict the stability characteristics. 



<-«—x=0.S ~X—X=0.8 


Fig. 1 - Variation of Rayleigh Number (J^i) with stable solute gradient {S\) for A. ~ 2, 6 - 0.5, P - 0.2, 
Qi ==3, = 20 for fixed wave numbers x = 0.5 and x - 0.8. 



— ♦-^xaO.5 — X— x=0.8 


Fig. 2 “ Variation of Rayleigh Number (i?i) with rotation ) for A. ~ 2, 6 - 0.5, P - 0.2, Qi - 3, Sj 
10 Tor fixed wave numbers x = 0.5 and x = 0.8. 
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Fig. 3 - Variation of Rayleigh Number (Ri) with magnetic field (Qi) for X = 2, € = 0.5, P = 0.2, 
=20, Si = 10 for fixed wave numbers x = 0.5 and jc = 0.8. 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 
P 

— x»0.5 

Fig. 4 - Variation of Rayleigh Number (Ri) with medium permeability (P) for X = 2, € = 0.5, P = 0.2, 
Qi = 3, =20, Si = 10 for fixed wave numbers x = 0.5 and x = 0.8. 

In Figure, 1, Ri is plotted against Sj for X = 2, s = 0.5, P = 0.2, =20, = 3 for 

fixed wave numbers x = 0.5 and x = 0.8. This shows that stable solute gradient has 
stabilizing effect. In figure 2, Ri is plotted against for X. = 2, e = 0.5, P = 0.2, Q\ 
= 3, Si = 10 for fixed wave numbers x = 0.5 and x = 0.8. This shows that the rotation 
has stabilizing effect. In figure 3, R] is plotted against Qi for X = 2, e = 0.5, P = 0.2, 

=20, Si = 10 for fixed wave numbers x = 0.5 and x = 0.8. This shows that magnetic 
field has destabilizing effect, for Qi “ 1 to 2 and has stabilizing effect for Q\ = 2to 20. 

In Figure 4, R, is plotted against R for X = 2, e = 0.5, Si = 10 =20, gi = 3 for 

fixed wave numbers x = 0.5 and x = 0.8, which shows that medium permeability has 
destabilizing effect, for R = 0.1 to 0.3 and has stabilizing effect for R = 0.3 to 1 . 
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stability of the system and oscillatory modes : 

We discuss oscillatory modes, if any, on the onset of stability problem due to the 
presence of magnetic field and rotation. Multiplying equation (37) by W*, the complex 
conjugate of fV, integrating over the range of z and making use of (38) - (42) together 
with the boundary conditions (45), we obtain 


- +— (l + cF) 
6 P/ 


47tpoU up 


up i_ e ^ 

+ ^ -[^9 + P2^^W ] = 

4rtpoU 


(54) 


where I, = \{dW^ +a^\W'\)dz, 


I, = J +aY\" + 2a^\DKf)dz, 


+a^\K\)dz, 

0 

+a^\@f)dz, 

0 


/.=j(jDr|’+o-|rr);ir, 

0 
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0 

0 

I, = \{dX^ +a^\x{)dz, 

0 


which are all positive definite. Substituting a = + ia,-, in equation (54) then equating 

the real and imaginary parts, we obtain 




i L 

i''’pj ' 47CpoU' 

^ 1 T?\ .. ^2 

d^L 


1 F 

■ + ■ 


up 


up' 


Pi 


ye Pi) 47tpoU 
p, e'n^ gaicfl^ ^ ^ga'K.'a'^ 

2 _ o ^ 


47tPoU up 


up' 




47rpoU 


(56) 


and c, 


^1 


Ir KS-n „ , gaica gaKfl , 

Yi --. — M ;v~^ 

' 47tP(,U up up 




‘J 

1 


' y 


+hi^pi 

‘t'P.u " 


= 0 


(57) 


Equation (38) yields that CTr may be positive or negative. The system is, therefore, 
stable or unstable. It is clear fi’om equation (57) that ct, may be zero or non-zero, 
meaning thereby the modes may be non-oscillatory. The oscillatory modes are 
introduced due to the presence of a magnetic field, rotation, viscoelasticity and stable 
solute gradient, which were non-existent in their absence. 
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The case of overstability : 

Here we discuss the possibility of whether instability may occur as overstability 
Since for overstability we wish to determine the critical Rayleigh number for the onset 
of instability via a state of pure oscillations, it suffices to find conditions for which 
(48) will admit solutions with cj real. 

Equating real and imaginary parts of (48) and eliminating Rj between them, we 
obtain 


+M+A^ = 0 

where we have put and 


(58) 


[e P J p 


6 P 


Epfi'Y, 


A. 


_ (l+y)' ( 1 F% 




+ 


J 


(l + x) 


, a ri .FTt^.ri , 


- + 

P’ eP 


+ (l + x)‘ 


s P 


J 


9LA 

r)3 


eP 


— + ■ 


O, 


4" • 
e P 


eP 




+ 


+ 




eP 


[ 


(63) 


Since ai is real for overstability, the four values of ci are positive. The product of 

the roots of equation (58) is y-, and if this is to be positive then > 0 since Jfrom 

A 4 

equation firom equation (59) A 4 > 0 . 

Equation (59) shows that this is clearly possible i.e. Ao is always positive if 

Epi > E'q , Epx > P2 and ^ ^ 

P e 


( 60 ) 
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which implies that 

E'k > Ek', k < jBt) and V) > 

Therefore, 

E'k > Ek', K<Er\ and u > 


r4fiV 

U.e J 

r4Qy 

li? J 

'UKPo] 

J 


are the sufficient conditions for the existence of overstability and the principle of 
exchange of stability is not valid. 


Conclusion 


The stable solute gradient is found to have stabilizing effect whereas rotation has 
stabilizing effect for X > 0 and destabilizing effect for ^ < 0 which is in contrast to the 
Newtonian fluids. The magnetic field has a stabilizing effect for 




1 + X + 


Pa 



— and as gravity increase upwards (i.e., X > 0) and .has 

PQiJ 

— and as gravity increase upwards 


1 + X + - 


' (l + x)?' 

(i.e., X. > 0) whereas medium permeability has stabilizing/destabilizing effect on the 


system for — < 





\ + x 




■ for X, > 0 and the same 


behaviour have also been shown graphically. The case of overstability is also 
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and u>| 

r4oy 



[hj 

U.sJ 


t] are the 


sufficient conditions for the existence of overstability. 
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Abstract 

The linear theory of generalized diermoelasticity without energy dissipation for homogeneous 
and isohropic microstretch elastic materials is employed to study the disturbance due to 
mechamcal and thermal sources in a half-space by the use of Laplace-Hankel transform 
echmques. The integral transforms have been inverted by using a numerical technique. The 
anmytical expressions of the displacement components, force stress, couple stress, microstress 
an temperature field are obtained. Normal force stress and tangential couple stress for 
mstantaneous and continuous mechanical and thermal sources have been computed numerically 
and illustrated graphically. Stretch effect on various expressions has also been obtained 
analytically and depicted graphically. 

(Keywords : Microstretch, Thermoelasticity, Sources, Integral transforms) 

Introduction 

Thermoelasticity theories which admit a finite speed for thermal signals (second 
sound) have aroused much interest in the last three decades. Recently, relevant 
theoretical developments on the subject of finite velocity of heat propagation are due 
to Green and Naghdi‘'^ which provide sufficient basic modifications in the 
constitutive equations that permit treatment of a much wider class of heat flow 
problems. An important feature of this theory, which is not present in other 

ermoelasticity theories, is that this theory does not accommodate dissipation of 
thermal energy. 

Micropolar elasticity” termed by Eringen'^ is used to describe deformation of 
eastic media with oriented particles. A micropolar continuum is a collection of 
interconnected particles in the form of small rigid bodies undergoing both 
translational and rotational motions. The linear theory of micropolar thermoelasticity 
was developed by extending the theory of micropolar continua to include thermal 
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effects by Eringen^ and Nowacki.* The theory of mi.crostretch elastic solids is a 
generalization of the micropolar theory and was introduced by Eringen.’’* The theory 
of thermo-microstretch elastic solids was developed by Eringen.® Recently many 
authors'®"'* discussed different type of problems in microstretch and thermo- 
microstretch solids. Following®’* we study the disturbance due to mechanical and 
thermal sources acting at a point on the surface z = 0 in the generalized thermo- 
microstretch elastic half space by applying integral transform techniques. 

Formulation of the problem and solution 

A homogeneous, isotropic microstretch generalized thermo elastic solid 
occupying the semi-infinite region z > 0, is considered in an undisturbed state and 
initially at uniform temperature Tq. We take the origin of the cylindrical co-ordinate 
system {r, 0, z) on the surface z = 0 and z-axis pointing vertically into the medium. Let 
J^r, z, /) be the change in temperature of the medium at any time. A mechanical 
source or thermal source is assumed to be acting at the origin of the cylindrical co- 
ordinate system having isothermal boundary. 

For two dimensional problem, we assume 

u=(u„0,u,), $=(0,<t)e,0) (1) 

and introducing dimensionless quantities defined by the expressions 



in the field equations derived by Eringen® and Green and Naghdi* and then inserting 
the potentials O, \|; and F defined by the relations 
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5<D aV 
dr^drdz’ 




dz 


( 


dz^ 




-ar 

dr ’ 


<I>6=-F-. (3) 


in the resulting equations, we obtain 


V^-a. 


'’a/' 


p - OgT + = 0 , 


(4) 


72 5 


2\ 




dt 


vi;+a4r = 0 , 


(5) 


a2 ^ 


-2a,- a 




r-a,vV = 0 , 


( 6 ) 


^72 


ar 


o2 

jj - e-^V^O = 0 , 

^ a/^ 


(7) 


a9r-agV^O + 


p)2 ^ 

r72 ^ 

a.V 


(t) =0, 


( 8 ) 


where 


_^Kcl 

y0 


^1 = 


_P-/C2 


a,= 


PC2 




K 


)i + K 


a< =- 


A# 4* 2fx + 


“.4 


«7- — 
0 


0 


2viC^ 

^9=— V 


2 . 




y Gq 2 > 

9v(0^j 


2 ^ + 2jj. + iiT 


c, = 


p^* ’ ar^ r dr dz^ 


(9) 
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Applying Laplace and Hankel transforms defined by 

0 

and 

f{q,z,p)=\f{r,z,p) rJ„{rq)dr , (10) 

0 

on equations (4)-(8) and then eliminating f,^' andf from the resulting expressions, 
we obtain 

[v'« + ^V"'+fiV'^+c][6]= 0 (11) 

and 

( 12 ) 

where 


A = -p^{aQ+d^+e a^) - dj +a^d^ , 

B= p^ + Oo + i/y )- e + 6 a„d^ - ] , 

C = - a^p^d^dj , E=-2q^ - 2a^- p^ [a^ + a^)+a^ai , 

F = [2a, + (aj + ^3 ) - 0,04 Ya^p^ (2a, + Oj ) , 






( 13 ) 
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The solutions of equations ( 11 ) and ( 12 ), satisfying the radiation conditions at 


infinity are 

6 = exp(- ^,2)+ A2 exp(- 42 z)+ 4 exp(- ^32) , ( 14 ) 

f =blAi exp (-4,2)+ 62^2 exp(-^22)+^>34(-^3^) . (15) 

I* = al Aj exp(- §,2)+ 0*2 A2 exp(- ^2^)+ 0*3 Aj exp(- §32) , ( 16 ) 

H/ = ^4exp(-^42)+^5exp(-^52), ( 17 ) 

f = 04^4 exp(- ^42)+ a’, As exp(- ^52) , ( 18 ) 


where 


t2 2 , ^ 

Sn ■^n 2 : 


= 2^-H cos 


(|) + 2rc(«-l) 


(«= 1 , 2 , 3 ) 


;2 

s4,5 ■ 


-e±4e'^~af 


V 3 

h, ^=tm-'{j\A\/-G), A=G^+4H\ 


G=^-^ + C, H=j-^, a*={^j -g^-p^a,)/(a,N,-as) 




HO* + -q^ - P^aX 

Oo 





d,\ 



(/ = 1 , 2 , 3 ) 


♦ 




-a^p^ 


1 


and ^,(j = 1 , 5) are supposed to have positive real parts. 


( 19 ) 
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Case I : Instantaneous Mechanical Source acting on the Surface 
The boundary conditions in this case are 

tzz = . t,,=ma=T=X^=0 atz = 0, (20) 

Inr 

where t^, t„, and have been borrowed from the constitutive relations derived by 
Eringen’ and Green and NaghdP, 5( ) denotes the Dirac delta function and P is the 
magnitude of the force. Making use of equations (2) and (3) in the boundary 
conditions (20); applying the transforms defined by equation (10) and then 
substituting the equations (14) - (18) in the resulting expressions, \ve obtain the 
expressions for displacement components, stresses, temperature field and microstress 
as 

L\o 

u, , ( 22 ) 

LAq 

C + /sAse'^^^] , (24) 

^ze =^^fe4^lA4e'^'‘" (25) 

LXq 

f = +blA,e-^^^ +fc3*A3e-^3^] , (26) 




+^2^2A2e”^^^+^3a3A3e ^^*], 


(27) 
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where 


-44a>,)+|3a3(^4a>3 -^5a5e4)]+ 

^2 kl«i* - ^5«5«2 )+ ^3«3 - ^4«Ie5 )] + 

*3 ks^l^A - 44«>3 )+ ^ 2«2 ^4^^!% “ ^ 5 « 5«6 )] » 

=Poei^iblal > ^2 =Poenk\^'xbl -libla\), 

= Po^i ^iK <h - ^ibW\ ) ’ ^4 = Po^s4b2^6 ’ 

Aj ^-P^^^alb^A^, Aj = - al)+ -aj )-bl%^%-^ (al -a,‘ ), 


fi={b\+bj^] -b^q^ -b]a„+b^a] , /= 1,2,3 ; 

fj =-b2^] +^4^ -Q^ybiOj , j = 4,5 ; P„ =%„ , 
= 73/5 +9^*2^3^5 . H = 73/4 +^f^*2^3^4 > «3 = fifs ^ blq^i^^ks » 
^4 =7274 +?^*2 ^2^4 > ^5 =7l75 + ?^*2 ^1^5 > ^6 =7l74 +9%^^1^4 , 


- 75^4^4 "■ /a^S^S > bo - °2 

PC2 


bi=-y 


2p + K 


■'2 ~ 2 
PC2 


A - 

O4 - j- 


^ ..LI. A Y ® A — — 

Ot = r- , 03 = 04 + 05 , O 5 = 7 , 07 — 2 

' PC^ PC2n^2j PC2^2 
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Particular Case I : 

By taking = X., = 0 in equations (21) - (27), we obtain the expressions 

for displacement components, stresses and temperature field in a micropolar 


generalized thermoelastic medium as 

IXo 

+A'e-^5^)], (30) 

CXq 

(33) 

f (34) 


where 


A1 = - Po^z U^jOs - / 5 ^ 4«4 ) . Aj = P„q^ - 75^404 ) , 

a; = PcksoXi^i^z -^ki) > A'j -^192) . 
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53=/.'/5+*2V|;^5. 






-A' ±4^^^^']/^, A'=-p\d,+a,+ ea,)-2g\ 


B' = q*+ q^p^{d^ + a„+ e aj+ p*a„d^ , 9 , 2 = — - 9 ^ -p^aj , 

«o 

//=(*!+ *2 , / = 1,2. (35) 


Case n : Continuous Mechanical Source acting on the Surface 
The boundary conditions, now take the form 




zMM) 

2 itr 


t^,=m 


•zO 


= r = A,j=0 arz = 0, 


(36) 


where H{t) is Heaviside unit step function. The expressions for displacement 
components, stresses, temperature field and microstress in this case are again given by 
equations (21) - (27), with Po in right hand side of A, (/=!,.... 5) be replaced by 

fp\ 

p:= ^ . 

KP) 

Particular Case 11 : 

Neglecting stretch effect, the expressions for displacement components, stresses 
and temperature field are given by equations (29)-(34) with Po in right hand side of 
A' (i = l,2,4,5)be replaced by jP* . 

Case in : Instantaneous Thermal Source acting on the Surface 


The boundary conditions in this case are 


= , T = 


-P&41 


2nr 


atz = Q , 


(37) 
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where P is the constant temperature applied on the boundary and the boundary is kept 
sfress free. With the help of these boundary conditions (37) and following the 
procedure adopted in Case I, the expressions for displacement components, force 
stresses, couple stress, microstress and temperature field are again given by equations 
(21)-(27) with A,</=1,....6) replaced by A* , 

where 

Aj = —Pg [^2^2 )"*■ ^3^3 (^5^5^4 ~ ^4*^4% )] ’ 

A*2 = Po - ^5«i«2 )+ ^3^*3 - lAes )] , 

A*3 = -Pg [^lOl (^4ale3 - ^ 5 ^ 5*64 )+ ^ 2^2 - ^4«I% )] > 

A*4 =-Po|jajZ)2^*6 > ^*5 " ■^o^4‘^4^2'^*6 ’ 

^6 if 2^3 ~ f3^2)~^2^2i/l^3 ~ ~ 

Particular Case in : 

Neglecting stretch effect (a^ =Xo =A,, =0), we obtain the expressions for 
displacement components, stresses and temperature field in a micropolar generalized 
thermoelastic medium as given by equations (29)-(34), only with A' (/ = 1,2,4,5) 

replaced by A° , 
where 

~~P}^5^5^2 “^4^4‘^l]» ^2 “~-^[^4^4‘^3 “45^5‘^4]5 

^4 ~ ~Po^3^5^2if\%2 “/2^l) > ^5 ~ Po^A^A^2if f)2 “ /2^l) • 

Case IV : Continuous Thermal Source acting on the Surface : 

The boundary conditions in this case, take the form 
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^zz ^zr 




atz = 0 , 


(40) 


The expressions for displacement components, stresses, temperature field and 
microstress are obtained in the similar manner as obtained in case II and again are 

given by equations (21) - (27) with A, (/ = 1, ,6) replaced by A* given by equations 

(38) and Po in right hand side of A* (/ = 1 S) be replaced by P/ . 

Particular Case IV 

If we neglect stretch effect, the expressions for displacement components, stresses 
and temperature field are given by equations (29)-(34) with the replacement of 
A' = 1,2,4, 5 ) by A® given by equation (39) and with Pq in right hand side of 
A" (/ = 1, 2 , 4 , 5 ) be replaced by P^* . 

Following Honig and Hirdes’’ and Press et al.^®, the transforms in equations (21)- 
(27) and (29) - (34) are inverted for all the four cases to obtain the solution of the 
problem. 

Discussion of Numerical Results 

We take the case of magnesium crystal^’ like material subject to thermal 
disturbance for numerical calculations. The physical constants used by us are 

p=1.74gm/cm\ 7■=0.2xl0■'’cm^ ?u=9.4xl0"dyne/cm^ 

p=4.0xl0"dyne/cm^, Ar=1.0xl0"dyne/cm^, Y=0.779xl0‘‘*dyne, 

A,o=0.5xl0’'dyne/cm^, A,i=0.5xl0'’dyne/cm^, ao=0. 779x1 0*'* dyne, 

.Sr*=1.0005xl0"cal dyne/gm°C cm^, C’=0.23caI/gm°C, /o=6.13 lxl0‘’^sec., 

ti=8.765xl0-‘^sec., e=0.073, ei=0.069, To=23°C, I = \ cm. 

The variations of normal force stress T^{=t^lP) and tangential couple stress 
M^{= / P) with radial distance ‘r’ at time t = .1 sec. have been shown by (a) solid 

line for cases I and HI, and dashed line for cases 11 and IV in microstretch generalized 
thermoelastic (MSGTE) medium and by (b) solid line with centered symbols for 
cases I and HI and dashed line with centered symbols for cases H and IV in 



NORMAL FORCE STRESS 
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micropolar generalized thermo-elastic (MGTE) medium. The variations for the above 
said cases (I-FV) are shown in Figs. 1-4. 



Fig. 1 -Variations of normal force stress 
^z 2 (“ ^zz ^ radial distance r (Case I, II). 


Fig. 2-Variations of tangential couple stress 
(= / P) with radial distance r (Case 1, 11). 




Fig. 3-Variations of normal force stress 
(= / P) with radial distance r (Case III, IV). 


Fig. 4-Variations of tangential couple stress 

M^(= / p) with radial distance r (Case in, rV). 
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The values of normal force stress in case I, vary in a similar manner in two 
media (MSGTE and MGTE) and in case II, values of are large in MSGTE medium 
in comparison to MGTE medium in the range 0 < /■ ^ 2.0 but are small in the further 
range (2.0 < r ^6.0) as shown in Fig. 1. The variations of tangential couple stress (Me) 
lie in the range -0.02 < M^< 0.08 in both the media and for both cases I and II and 
the variations have oscillatory behaviour as depicted in Fig. 2. Fig. 3 illustrates the 
variations of Tiz for the cases III and IV in both the media. Due to stretch effect the 

values of in MSGTE medium are large in the range 0 < r < 1.8 and 3.5 ^ r < 5.0, 

small in the range 1.8 < r < 3.5 and 5.0 < /• < 6.0 in comparison to MGTE medium for 
case in, whereas for Case IV, values of Tz- in MSGTE medium are small in the range 
0 < r < 2.55, in comparison to MGTE medium and thereafter attains almost same 
values for both the media. The slight oscillatory behaviour of Mze is shown in Fig. 4 
for the cases III and IV and for both the media. Thus, from the above numerical 
results, we conclude that stretch has a significant effect on normal force stress and 
tangential couple stress in all the four cases. 
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Abstract 

An analytical solution to the steady state three dimensional atmospheric diffusion equation with 
multiple sources has been developed. The modelling is done using Non-Gaussian model, where the 
wind speed and eddy diffusivities are assumed to vary with the power of vertical height which is 
uniformly assumed to be 1/2. The removal mechanism is also incorporated in the model and Green’s 
function concept has been applied for the solution of the problem. This model is useful in the prediction 
of pollution concentration not only from a single point source, but also from multiple point soruces. 

It also helps in studying the effect of sink mechanism in decreasing the concentration of pollutants 
from the atmosphere. 

(Keywords : Multiple sources/variable wind velocity/variable diffusivities/removal rate/Green’s 
function) 


Introduction 

In past, many steps have been raised to understand the dispersion of air pollutant in 
the atmosphere under various meteorological conditions. An atmospheric dispersion model 
is needed in order to calculate the spatial distribution of pollution concentration due to the 
emission source. Many investigators^’^ have tried to develop such models considering 
different meteorological conditions, and adopting different methods of solution. Generally, 
the Gaussian^ and the non-Gaussian type models are in practice. When compared to Gaussian 
models, the non-Gaussian models have wider applicability. 

Analytical solutions of the atmospheric diffusion equation using non-Gaussian type 
models have been found in literature since the 1950s. Earlier, the solutions of the 
non-Gaussian type models were found out, considering only the two dimensional cases, or 
neglecting the inversion effect. By applying the Green’s function method, solutions for 
three dimensional diffusion equation under both bounded and unbounded boundary 
conditions have been obtained^. But in the above models, the presence of several removal 
mechanisms in the atmosphere have been not taken ito account. In the environment various 
mechanisms such as rain out, wash out, deposition oh vegetative canopies^ etc. are known 
to exist which are effective in removal of pollutants from the atmosphere. 
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In the light of above, we have analysed the problem of steady state transport of a 
contaminant, emitted from multiple sources, under Dirichlet’s boundary conditions, 
considering the effect of removal mechanism present in the atmosphere and using the 
non-Gaussian plume equation, to describe dispersion. 

In order to simplify the problem, we have considered the case of equal power profiles 
of wind speed and diffusivities which is assumed to be 1/2. The solution of this problem 
is obtained utilizing the Green’s function concept. 

Diffusion Equation 

The partial differential equation describing the steady state dispersion of a non-reactive 
contaminant released continuously from ‘n’ point sources located at (x* , y \ , z]), 
(Xj , y\ , Zs) — (^ . y" » z") in a Cartesian co-ordinate system, where the removal of the 
contaminant which is assumed by some removal mechanism in the region of interest is 
given by 


U(z) 


dC {x, y, z) 
dx 


d_ 

dy 


{x, z) 


dC jx, y, z) 
dy 


+ 


dz 






dz 


+ X e'5(^-4)8(y-y.)8(z-4)-/f,C(^,y,z) (1) 


where C (jc, y, z) is the ambient concentration of the contaminant. Q is the emission strength 
of /* source located at , 4)> 8 is the Dirac delta function* with the property 

J 6 (x - Xq) tir = 1 and K^ is the removal rate of the contaminant. 

R 


Here the wind is assumed to be blowing in the x-direction, the turbulent fluxes are 
approximated by gradient transport (K-theory) and the turbulent diffusion is neglected 
compared to advection (Slender plume approximation) since the scale of the turbulent 
transport is smaller than plume dimensions^. 

The wind speed U{z) and eddy diffusivities Ky (jc, z) and (z) are assumed to vary 
with height and are taken as the following power law profiles. 
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t/(z) = 


If- 


^ _ K„Kx)z^ 


K, 


if 


where J{x) is an integrable function of ;c and U[f and Kf^ are the wind speed and eddy 
diffusivity respectively at reference height. H. For simplicity in our paper, we have assumed 

a=p = Y= 1/2. 

Equation (1) is reduced to 


f = I 


dC 


H 

V J 


' 3y I 9z 






+ S i2'S(x-4)5(y->/)5(z-4) 

7=1 


(2a) 


Boundary Conditions 

It is assumed here that the pollutant is immediately removed on touching the boundaries 
(Dirichlet’s condition) 


C(x,y,z) = 0 

, X OO 

2(b) 

C(x,y,z) = 0 , 

, y —> ± oo 

2(c) 

C(x,y,z) = 0 , 

, z = 0 

2(d) 

C(x,y,z) = 0 , 

, z = H 

2(e) 


where H is the height of inversion layer. The dispersion equation and it’ s boundary conditions 
are made dimensionless by considering the following dimensionless quantities : 
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- Z - 

H'""' UffH^ 



U„H^C 

~Q 


, K = 


-KiH^ 

Kh 


U 

5(J-?,) = 6(x-4) 


^(y-t) = 8 (y-i) H, 


Z(z-t) = 5(z-4)H 


7® =m, e'' = -f 


On removing bars, dispersion equation (2) and its boundary conditions are changed to the 
form : 



dy 


Ax)z 


'4 K 
By 




Bz 


+ KC + X 2'8(j:-Jci)8(>’-ys)8(z-4)» 
;=i 

C (j;, y, z) = 0 , X -> oo 
C{x,y,z) = 0 , y ±oo 
C(x,y,z) = 0 , z = 0 
C (x, y, z) = 1 , z = 1 


(3) 

4(a) 

4(b) 

(5) 


Analytical Solution 

This boundary value problem is solved following the procedures of Yeh and Huang^ 
and Yeh^ for a single isolated source. The solution is comprised of three components, a 
source strength a vertical dispersion factor, the vertical sub-Green’s function 
(x, z, , 4) and a cross wind dispersion factor - the cross wind sub-Green’s function 

G^y(x,y,4,yi) . 

Thus solution of equation (3) is 

n 

C{x,y,z) = X 
y=i 


( 6 ) 
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where the two dispersion factors and are both equal to zero if x > . The product 

of these two factors (or sub Green’s functions) constitutes a Green’s function^^’^ ^ which 
can be interpreted as the concentration response due to unit disturbance (source) at 

(4 - - zi)- 

The total concentration at (x, y, z) is, therefore, the sum of the responses from ‘n’ 
various sources. 


Method of Solution 

For the concentration response due to a single source, equation (6) can be written for 
; = 1 as - 


Cix,y,z) = cc c Gz . Zo ; z) . Gy (xq ,yo-,x, y) 
• 8 (jCq - 5 CVo “ 3'i) 5 (zo - Zs) dx^dy^ dz^ , 


Cix, y, z) = QGj (a:, z ; , z,) Gy (x, y;x^, y,) 

where (aq , Zq ; x, z). Gy (jfg ^y^^x, y) is the Green’s function that satisfies the following 
pairs of adjoint partial differential equations'^’*^. 


dG, dG, 

Zo^ + Zo"^ ■^ + = ■'S(^o“^)5(zo-z),^o<^ 


(7) 


and Gj (Aq ',Zq\x,z) = 0 for Aq > a 


and 




8(ao - x)5(yQ-y) ,Xq<x 



and Gy (aq , Zq ; a, z) = 0 for Aq > a 

^0 > % ’ Zq used as independent variables instead of a, y, z so that finally we get 
C(a, y, z) after integration. For any location Aq , equation (7) can be expressed as 
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dzo 




so. 


lA 


lim G, = S(zo - z). 


( 9 ) 

( 10 ) 


The method of separation of variables is used and on setting 
Gj (xq ,Zo) = M(zo).N (xq) , we have 

Sm I dM 


d zl 2j:o dzQ 


+ (Kz^'^ + t]^)M ^0 


( 11 ) 


and 


1 dN 2 

Nd7r^ 


( 12 ) 


where is the separation constant. The method of power series is employed to solve 
equation (11). 

M{zq) = al { 1 - 0.6667 - 0.333q^z^ + 0.0899 + 0.0952 Kr\^ ^ 


- 0.0952 - 4.938 x 10~^ - 8.1834 x 10"^ K^z^ - 4.329 x lO"-" 

+ a\ {2z^ - 0.6661 Kzl - 0.4 + 0.0635K^^ + 0.0162Ky\^z^ 


fA-3 


+ 0.0222q'‘ - 2.822 x 10”^ K^4 ~ 5-0794 x 10"^ 


+ 0{z^) 


(13) 


The solution of equation (12) has been found as 


Nixo) = C, 


(14) 
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(^0 ’ Zo) = exp {- Ti^ (x -Xq)} [qo (1 - 0.6667 ) 

+ a, (2^;^ - 0.6667 )] (15) 

We have applied the conditions (5) and (10) to equation (15) and replacing variables 
to get the vertical dispersion factor due to the source at (^ , }/ , z^) as 

Mj (z) M,. (4) exp {- til (x - 4) } 

1‘z'^ {M,{z)fdz 

where r\^ is given by the eigenvalue equation 

2 - 0.6667 K - 0.4 Ti^ + 0.0635 + 0.0762 Kx]^ + 0.0222 

- 2.822 X 10"^ - 5.0794 x 10”^ r\^ « = 0 (17) 

and 

M.{z) = 2z'^ - 0.6661 Kz^ - 0.4 r\jz'^ 

+ 0.0635 K^z^ + 0.0762 A-n/ + 0.0222y]f z"'^ 

- 2.822 X lO'^A^z^ - 5.0794 X 10"^ A^Tifz"'^ 

+ 0(z®) (18) 

Now, equation (8) can be written as 



G,(x,z;x',4)= S 

i=i 


( 19 ) 
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lim Gy (xq , ; -c. j) = 8 CVo “ 3') (20) 


The Fourier transform is employed to solve equation (19) and (20), to get 


Gy(:ro,yo;-*:.y) 


exp 


-(y-yo>^ 1 

4 

jf /(<t>o)#o- 1 

' /(<t>o) #o| 

V4^ 1 1 

r 

r /(<i)o) #o| 


( 21 ) 


where is some reference location in Cartesian co-ordinates. To find an explicit jf(;c), we 
have used the relation (as is frequently done in dispersion modelling) as 


i.e. 




fix) 


1 da^ 

2 dx 


(22) 


where cj is the mean square particle displacement. 

For multiple source modelling, we have assumed thatj((x) depends on the relative distance 
between source and the location of interest. 

We have put (22) in (21) and replaced variables to get the cross wind dispersion factor 
due to the/'' source at (4 , , 4) as 


Gj,(x,y;x' ,)/) = 


1 

V2Jco^(x-xi) 


. exp 


-cv-y/ 

2oJ(x-x') 


(23) 


where Gy is standard deviation. 

We have finally obtained the total concentration response at (jc, y, z) due to n sources 
located in the region of interest as - 
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C(x, y, z) 


ll 


s 


oo 


G'- 


s 






M,.(4)exp {-Tif (x-4)) 



2o^(x-4) 




exp 

\2na^(x-x'^) 



Discussion 


(24) 


The case of a single elevated point source at dimensionless parameters (0, 0, 0.2) is 
discussed here. The dimensionless removal rate constant K = - 2.77, standard deviation 
Oy = 1.878 , and = 1. 

To illustrate the behaviour of this three-dimensional atmospheric dispersion model, 
concentration profiles have been calculated for a variety of conditions using equation (24). 



z 

Fig. I - Dimensionless concentration C(X, 0, Z) plotted against dimensionless vertical height Z. 


In Fig. 1, the concentration of pollutants is plotted with respect to vertical distance (0 
< 2 < 1) for different down wind distances (x = 0.02, 0.04, 0.06), keeping the crosswind 
distance constant. It is seen that for a particular downwind distance (x = 0.02), the 
concentration profile attains it’s peak at z = 0.45. But with increasing down wind distance 
(x = 0.04, 0.06), concentration decreases, and the concentration profiles approach uniform 
distribution. This results primarily from horizontal spreading. As a consequence, the 
concentration level of pollutant near the ground becomes high, which may prove dangerous. 
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Fig. 2 - Dimensionless concentration C(X, 0, Z) plotted against dimensionless vertical height Z. 

In Fig. 2, the concentration is plotted with respect to vertical distance for different 
down wind distances (jc = 0.2, 0.4, 0.6). It is seen here that at low down wind distance (x 
= 0.2), the plume attains it’s peak, while for higher down wind distances (x = 0.4, 0.6), 
there is a uniform distribution. 



Fig. 3 - Dimensionless concentration C(X, 0, Z) plotted against dimensionless downwind distance X. 

In Fig. 3, concentration is plotted against the down wind distance (0.02 < x < 1) for 
different vertical distances (z = 0.1, 0.2, 0.3). The extensive horizontal spreading shown in 
this graph confirms the deductions made in (1) and (2). 
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Fig. 4 - Dimensionless concentration C(X, 0, Z) plotted against dimensionless downwind distance X, 

In Fig. 4, concentration is plotted agaisnt the cross wind distance (- 03<y< 0.3) for 
different values of down wind distances (x = 0.02, 0.04, 0.06). It is seen here that in the 
cross-wind direction, the concentration profiles are symmetric with peak concentrations 
along the centre line of the plume. It is also seen that concentration levels decrease and 
there is marked lateral spreading with increasing down wind distance (x > 0.4) from the 



Fig. 5 ” Dimensionless concentration C(X, Y, 0.2; plotted against dimensionless crosswind distance Y. 

In Fig. 5, concentration is plotted with respect to down wind distance (0,02 < x < 0.1) 
for different cross wind distances (y = 0.1, 0.2, 0.3). For a particular cross wind distance 
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(y = 0.1), concentration is seen to decrease regularly with the down wind distance. For the 
cross wind distance (y = 0.2), concentration profile is almost uniform for down wind 
distance (0.02 < .r < 0.04) and then decreases regularly for x > 0.04. 

For the cross wind distance (y = 0.3), it is seen that the concentration profile is a curve 
which attains it’s peak at (x ~ 0.045) with extended horizontal spreading in the increasing 
down wind direction. 


Summary 

A steady state three dimensionless non-Gaussian model with multiple sources is 
developed, taking into account the various sink mechanisms present in the atmospheric 
which are effective in the removal of the pollutants. An analytical solution to this problem 
is obtained using the concept of Green’s function, he diffusion equation is broken into three 
components - a source strength, a cross-wind dispersion factor and a vertical dispersion 
factor. 

The model is useful in the prediction of the harmful effects of pollutant from multiple 
sources and its removal by the natural or artificial sinks present in the atmosphere. 
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Abstract 

An interdigitai microstrip line has been analysed and its image impedance and phase 
constant 0 determined as a flmction of frequency, for the given microstrip dimensions and 
substrate dielectric constant, for both even and odd numbers of strips. An open circuit Z-matrix 
was obtained between the voltages and currents of the input and output ports in terms of the 
coupled microstrip even-and odd-mode characteristic impedances Zoe, Zoo, and velocities Ve, vq, 

Zit, and 0 were then determined from the matrix in each case. The results were compared with 
those obtained for meander lines, complementary to the corresponding interdigitai lines. The 
curves of-Z//, and 0 as a function of frequency were also plotted and compared. 

(Keywords : MICs/microstrips/interdigital line/periodic structures). 


Introduction 

In a previous paper by Agrawal’ the analysis of n-coupled microstrips was carried 
out taking all the couplings into account. This was in contrast with the 'imit cell' 
approximation used by other authors^'^ wherein the coupling between the unit cell and 
other adjacent strips were not considered, thus giving erroneous results. 

In this paper, the 'interdigitai' line (which is complementary to 'meander' line) has 
been analysed for « = 3, 4 and 5 coupled microstrips. The iterative impedance Zit and, 
phase shift 0 of the interdigitai line as a function of the even-and-odd-mode 
characteristic impedances Zoe and Zoo, phase velocities of propagation v* and vo, and 
firequency /have been derived for each value of n. The Z,f and 0 of the interdigitai line 
have been compared with those of the 'corresponding' complementary meanders 
previously obtained by Agrawal^ 
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Open Circuit Impedance 2^Matrix for a Coupled Line Section 


Fig. 1 shows a coupled line section. Assumed fictitious current sources drive the 
line in the even-and odd-modes, and are denoted by small letters. The port voltages 
and currents are also defined in the figure. By superposition, the total port currents, I, 
can be expressed in terms of the even-and odd-mode currents as 


■“ he ■“ ^30 

/3= ho 

he^ ho 

h ~ he ^*40 


‘ yCvI 

va ■*'V2 j 

* 

he — 

12-^ 2/V- 

— I2e 

*' V 

Zoe.Zo* '4Q^7 


l1e--^/3 — 13 

T 

1 

1 

1 

ize 

7^ 

Vb +V4/ 


1 



0 

1 



( 1 ) 


Fig. 1-A parallel coupled line section with even and odd-mode fictitious current sources. 

The terminal voltages and currents on the strips are shown by capital letters. 

First we consider the line as being driven in the even-mode by the iu current 
sources. If the other ports are open circuited, the impedance seen at ports 1 or 2 is 

z: = -/z„.cotp; (2) 

where % is the phase constant in the even-mode. 

The voltage on either conductor can be expressed as 

f;(z)= f;(z)= 2k;cosp.(i -z) O) 


so the voltage at port 1 or 3 is 


^;(o)=f^;(o)=2F;cosp.i=/,.z: 


(4) 
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Using (2) and (4), V* is substituted in (3) to get 


-yZp, Cos /,, 

Sin 


Similarly the voltages due to the current sources Z 2 e are 


F/(z)= F/(z) = -yZoc 


Cos Q.z . 

Sinp,l' 


(5) 


( 6 ) 


And the odd-mode voltages due to the current sources ho and /40 are 


V,\z)=V,\z) = -jZoo 


Cos Po(l-z) . 
Sin Pol 


V^\z)=V,\z) = -jZoo 


Cos Po^ 
Sin Pol 


*40 


(7) 

( 8 ) 


The total voltage on a strip a or 6 is equal to the sum of voltages due to all the 
current sources on that strip. 

From (1), he fictitious current sources can now be replaced by the terminal 
currents as 


h=ii,+h)/2 h,= (i2+m 

ho = {h - /i)/2 /40 = (/4 - h)l2 (9) 

The total voltage at port 1 is 


F, = FJ(0)+F/(0)-hFJ(0)-hF;(0) 
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= - jZShfiot p,l + /,Csc p.l) + ;Z,(/;„Cot P„1 + z;„Csc p„i) 


Using (9) we get 

2F, = -jZoe [(/; + h) Cotp.l + ih + h) Cscp.l] + 

yZoo [(/3 - h) CotPol - (/4 - h) Cscpol] 

or V, = -;72[/i {Zq^ Cotp.l + Zoo Cotpol } + h{Zoe CscP.l + Zoo 
CotPol} +/3{Zoe CotPel - Zoo CotPol} +/4{Zoe CscPel - Zoo 
CotPol} (10) 

This result yields the top row of the open-circuit impedance matrix [Z] that 
describes the coupled line section. From symmetry all other matrix elements can be 


found once the first row s known. The matrix elements are then 

Z\\ = Z22 — Z33 = Z44 = (ZoeCotpel + ZooCotPol) (1 Is) 

Z12 — Z21 — Z34 = Z43 = ~y 72 (ZoeCscPel + ZooCscPol) (11b) 

Zi3 = Z31 = Z24 = Z42 = -y72 (ZoeCotpel — ZooCotPol) (11c) 

Zi4 = Z41 = Z23 = Z32 = "772 (ZoeCscPel — ZooCscpol) (1 Id) 


A two port network can be formed from the coupled line section by terminating 
two of the four ports in either open or short circuits depending on the boundary 
conditions. The image impedance Z„ and dispersion 0 of the two port network can 
then be found out. 


A. Coupled Microstrip Interdigital Lines 

3-Coupled Microstrip Interdigital Line 

The current source excitations of 3-coupled microstrips a, b and c of an 
interdigital line are shown in Fig. 2. Proceeding as in the case of a pair of coupled 
microstrips, the terminal voltages can be written as — 
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INPUT 


c 




II 

Lli 


V2 


Vi 


T 

^4 


I3 


L^5 


* 


"6 


Vc 


OUTPUT 


Fig. 2~Interdigital line structure for 3-coupled microstrips. 

Fi = -yZofiCscpel (/uCosPel + lie) +yZooCscpol(/3oCosPol+i4o) (12a) 
V2 =-yZoeCsCpel (/le + i2eCosPel) -7Z00 QscPe 1 (/30 + /4oCosPel ) ( 1 2b) 
Vi = -yZoeCscPel (luCoSpel + he) 

-jZoo CscPol [(/30 - <5o)CosPel + (/40 - ^6o)] (12c) 

F4 = -jZoe Cscpj (he + he CosP.l) 

-jZoo CscPol [(/30 - ho) + (^40 - ho) CosPol] (12<1) 

Vs = -jZdeCsc^^l (/'uCosPel + he) -yZooCscPoK/soCosPo 1 +/6o)( 1 2e) 
= -jZoe Cscpel (ile+ /2eCosPel) -yZooCsCpol(/50 + /eoCoSpol) (12f) 


where 


/u= 1/3(/,+/3+/5) 
ho = 1/3 (h + 7s - 27i) 
/50= l/3(2/5-/l-l3) 


he- 1/3(12 + h+h) 

Uo= 1/3(74 + 76-2/2) 
ko= 1/3(276-72-/4) 


( 13 ) 
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The boundary conditions for the itnerdigital line of Fig. 2 are 


/3 = 0 = 741 V2 = Vs 

h = -h (14) 

Applying the boundary conditions in (12) and (13) and simplifying we obtain a 
[z] matrix between the input and output quantities as : 



(15a) 


where 


Zii = Z55 = -y/3Zo.Cotp,l + yZoo{Csc(2f3ol)-2/3Cotpol} 
Zxs = Z 51 = -y/3ZoeCotp,l -yZoo{Csc(2pol)-l/3Cotpol } 
The iterative impedance of the 3 strip interdigital line is 

^)t(3/) = “-^15 


and dispersion 

0(3/) = Cos'* 


^ 7 ^ 

±55 


(15b) 

(15c) 


(16a) 


(16b) 


4-Coupled Microstrip Interdigital Line 


The 4-coupled microstrip interdigital line with current source excitations is shown 
in Fig. 3 
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INPUT 



Fig. 3-Interdigital line structure for 4-coupled microstrips 
The terminal voltages are given below : 

Fi = -y'Zoe (/ieCot(3el+/2eCscpel)+yZoo(i3oCotPol+/4oCscpol) (17a) 

V 2 = -jZoe (/ieCscPel+/2eCotPel)+yZoo(/3oCscPol+/4oCotPol) (17b) 

F3 =-yZoe (/'leCotPel + /2eCscPel) 

~jZoo [(^30 ~ / 5 o)CotPol + (/40 — 4o)] CsC Pol (17c) 

F 4 = -_/Zoe (/’leCsCpel + hgCot^e^) 

-jZoo [(*30 - i 5 o)CscPol + (*40 - *6o)] Cot Pol (17d) 

F5 = -yZoe(*leCotpel + Z2eCscPel) 

—jZoo [(/50 — /7o)CotPol + (/’eo ~ /8o)Csc Pol] (17e) 

V6 = -JZoe (/’leCscPel + /2eCotPel) 


-jZoo [(/50 - / 7 o)Cscpol + (/‘eo - i&o)Cot Pol 


(17f) 
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F7 =-_/Zoe(/leCotPel+i2eCsc|3el)-yZoo(2loCotPol + /goCsCpfll) (17g) 
Fs =-yZo« (/leCscPel+/2eCotpel)-yZoo (iToCsCpol + /goCotPol) (17h) 

where 


he = (/,+/2 + / 5 +/ 7 )/ 4 ; 
ho= (h + h+h-31i)/4-, 
ho= {h + Ii-h-h)l2-, 
/70= (3/7-/i-/3-/5)/4; 

The boundary conditions for a 4-i 


he- {h + h + h+ h) 14 
Uo= {h + h + h-3h)l4 
ieo - {h + h-h- U) /2 
780 = {3h-h-U-h)l4 

interdigital line are : 


h=-h; V2 = Ve 

h=-h; F3 = F7 

/4=-/5 =0 


(18) 


On solving the above equations we get an open circuit (Z) matrix between the 
input and output ports ; 



(20a) 


where 


Zn = Zgg = -774[ZoeCotpel + Zoo{4Csc(2pol)-3CotPol] (20b) 

Zig = Zgi = — y74{ZoeCscPel + ZooCotpol} (20c) 


and 



(21a) 


and dispersion 
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0(4/) 


= Cos ' 


^88 


(21b) 


5-Coupled Microstrip Interdigital Line 

Fig. 4 shows the interdigital line structure made up of 5-coupled microstrips. Here 
also we apply the current sources from both sides of microsrips. The terminal voltage 
can be written as ; 



Fig. 4-Interdigital line structure for 5-coupled microstrips 
Vi = -yZoe(/iieCotPel+/2.Cscp,l)+yZoo(/3oCotPol+/4oCscPol) ( 22 a) 
F2 = -yZoe (Z,eCscPel+/2.Cotpel)+yZoo(/3oCscpoH-f4oCotpol) ( 22 b) 
Fs =-yZoe (/leCotPel + /2eCscPel) 

+ y'Zoo [(/'so ~ i3o)CotPol + (ieo ~ /40) CscPol] (22c) 

F4 =-yZoe (/‘leCscPel + /2eCotPel) 

+ yZoo [(/'so ~ i3o)CscPol + (ieo ~ /40) CotPol] ( 22 d) 


Vs = -yZo« (ileCotpel + i2eCsCpel) 
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+ y'Zoo [(^70 “ /5o)CotPol (^80 ~ ^6o)Cscpol] (22e) 

F6 = -7Zoe(/leCscPel + /jeCotpel) 

+7^00 [(^70 ~ /5o)CscPol + (/’so ~ *6o)CotPol] (22f) 

K7 =-yZoe (/leCotP^l + /2eCscPel) 

■’■y^’oo [(^90 “ /7o)CotPol + (/'lOo “ i8o)CscPol] (22g) 

K8 = -7Zoe(/leCotpel + laeCscPel) 

■*■7^00 [(*90 - *7o)CscPol + (/]0o - /8o)CotPol] (22h) 


F 9 = -yZoe(/le-CotPel+/ 2 «CsCpel)H 7 Zoo(/ 9 oCotPol+iOoCscPol) (22i) 
V\Q- -jZoe (/leCscPel+/ 2 eCotpel)+ 7 Zoo(/ 9 oCscPol+/lOoCotpol) (22j) 

where 

he = (-^1 +h + l5 + 19 ) /5; he = ih + 4 + /e + 4+ ho) /5 

ho ~ (-(3 + /s + -^7 + -^9 “ 4 / 1 ) /5; Uo = (4 + 4 + 4 + ho - 44) /5 

/50 = [2(4 + 4 + 4)-3(4 + 4)] /5; *60 = [2(4 + 4 + hoh ^{h + 4)] /5 

ho = [3(4 + 4) -2(4+ 4 + /i)]/5; *8o = [3(4 + 4o) -2 (4 + 4 +4 )]/5 

770= (44 -4-4-4 -/i)/5; /8o= (44o-4-4-4-4)/5 '(23) 

The boundary conditions for the 5-coupled strips interdigital line shown in Fig. 4 
are : 

4= -4; V2 = v, 

4= -4; F3 = F7 

4=-4o; F6 = F,o 

4 =4=4=0 


( 24 ) 
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Substituting (23) into (22) and using (24) we get an open circuit [Z] matrix 
between the input and output quantities as given below : 



(25a) 


where Zn = Z99 =? -^■Zoe/5[CotPel + Csc(2Pel)] 

+ yZoo/5{4Csc(2pol)Pol- 4Cotpol)} (25b) 

Zi9 = Z91 = (-yZo,/5) Cotpel) + (/Zoo/5 )Cotpol (25c) 


The interative impedance Z/, and dispersion 0 are then obtained from : 


Zit(5/) - — Z, 


6(5/) — Cos 


-I 


^ Z ^ 
^99 

V^19 J 


(26a) 

(26b) 


B. Coupled Microstrip Meander Lines 

Crampagne and Ahmadpanah^ had published the Zu and 0 curves for a 'unit cell' 
BCDEF of coupled microstrip meander as shown in Fig. 5. The meander BCDEF is a 
series connection of BCD and which are uncoupled. Agrawal' on the other hand 
showed the calculated curves for the meander ABCDE whose input and output ports 
are located at ;c = -LI2 instead of x= 0. These curves for a pair of coupled microstrip 
meanders are shown in Figs. 6(a) and 7(b). The curves obtained by Agarwal (solid 
lines) reveal that they are not independent of x, the location of the input and outpiit 
ports of the unit cell. According to Weiss^ and Crampagne and Ahmadpanah the 
phase 0 had been assumed to be the same for all x, i.e., 

V,(x) = Vi{x)e>° 


(27) 
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Unit cell 


Fig. 5-A unit cell of the periodic array of meanders. 

Agrawal’ showed that their assumption was invalid (by taking all the couplings 
between the meanders into account and) obtained the Zu and 9 curves for meanders 
consisting of both the even n and odd (n-1) numbers of lines. 

In this paper the same line dimensions and the same variations of Zoe, Zoo. Ve, vo, 
Pe, Po with respect to frequency as assumed by Crampagne and Ahmadpanah^ have 
been taken which were plotted with respect to frequency as shown in Fig. 8 by 
Agrawal’. The strip length 1 was 0.01836m. The Zunnd 0 curves as obtained from the 
above analysis of intergdigital lines have been compared with the corresponding 
complementary meanders of Agrawal’. The results are discussed below. 

Results 


3<Microstrip Interdigital Line 

Fig. 6 shows the Zt curves of 3 strip interdigital line and 2 (solid curves) and 4 
strip meanders, and Fig. 7 shows the 0 curves. The 3-strip interdigital line resembles 
the 2-strip meander except that the pass band of interdigital line is less broad and cut 
off is sharp. The pass band of 4-strip meander is smallest, and in this case there are 
two pass bands, within the frequencies of interest. Also, the Z„ curves of 4 strip 
meander are in the opposite sense. 




0 0-4 08 1.2 1-6 2.0 2.4 2.8 

f (GHz) 


Fig. 6(c) 

(Fig. 6b same as Fig. 1 1(a)) 
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Fig. 7(a) 



Fig. 7(c) 

Figs. 6 & 7-Iterative impedance Zu vs. frequency plots of 3-strip interdigital line compared with 2- and 4- 
strip meanders; and corresponding dispersion 0 vs. frequency plots. 

(Fig. 7b same as Fig. 12(a)) 
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Fig. 8-The characteristic impedance, plmse velocity and phase shift constant 
of the unit cell vs. frequency. 

4-Microstrip Interdigital Line 

Figs. 9 & 10 show the comparison of 4-strip interdigital line with 3 and 5 strip 
meanders. The interdigital line has 2 very narrow pass bands. The 3-strip meander has 
very wide pass bands. The curves for meanders are very different from that of the 
interdigital line. The number of passbands increases with ninnber of strips. The 
interdigital lines have sharp cutoff. 
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Here also the interdigital line has a very large stop band and its curve does not 
resemble those for the complementary meanders, being of opposite sense. The 6-strip 
meanders has the widest pass band. The interdigital line has the narrowest passband 
and very sharp cutoff characteristics. 

The aforesaid results can be summarised as follows : 

As compared to the corresponding complementary meanders the curves for 
interdigital lines are of opposite sense. By increasing the number of strips the munbers 
of pass and stop bands increase in both the meanders and interdigital lines. The latter 
has less number of pass bands in comparison. Odd pairs (2,6,10) of meanders have 
wider pass bands and even numbers of interdigital lines have narrow pass bands and 
wide stop bands. These results can be used in the design of filters and other periodic 
structures. Typical Z,f values as a fimction of frequency are given in the table below — 


Number of Zi,(ohms) 


Strips 


Interdigital line 



Meander line 


Zo 

4 

AGHz) 

3 

4 

5 

2 

3 

4 

5 

6 



0.5 

71.5 

00 

9.5 

79 

88 

82 

94 

00 

57 

132 

1.0 

68.5 

10 

1.4 

72.5 

114 

00 

44 

56 

58 

109 

1.5 

62.25 

CO 

00 

67 

52 

62 

72 

63 

63 

83 

2.0 

75 

00 

00 

64 

104 

88 

71 

86 

87 

62 
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Abstract 

The peculiar changes - characteristics of mesomorphic behaviour which occur at phase 
transitions, are primarily governed by the intermolecular interactions acting between sides, 
planes and ends of a pair of liquid crystalline molecules. In view of this fact, intermolecular 
interactions between a pair of MBA molecules have been evaluated using modified second 
order perturbation treatment along with multicentred-multipole expansion technique. An all 
valence electron method, CNDO/2, has been employed to compute the net atomic charge and 
corresponding dipole moments located at each atomic centre of the liquid crystal molecule. 
Using the results of stacking and in-plane interaction energy studies, probability calculations at 
varied angular and positional configurations in a molecular pair of MBA, have been carried out 
using Maxwell-Boltzmann formula. An attempt has been made to elucidate the nematogenic 
behaviour of 4'^-methoxyben2ylidene ~ 4 *- acetoxyaniline (MBA) molecules. 

(Keywords : nematic liquid crystal / molecular ordering / mesogen / CNDO / multicentred - 
multipole) 


Introduction 

Liquid crystalline phases are stable condensed phases in which molecules pack 
together with an order that is intermediate between the three-dimensional order of a 
crystalline solid and the disorder of an isotropic (ordinary) liquid. Liquid crystalline 
phases always have partial orientational order of the molecules and some liquid 
crystalline phases also have partial positional order of the molecules. The partial 
molecular ordering that is a characteristic of liquid crystallinity occurs frequently in 
natural and synthetic materials. Thus, liquid crystals are of considerable basic and 
applied interest'’^. 

The peculiar changes-characteristics of mesomorphic behaviour which occur at 
phase transitions, are primarily governed by the nature and strength of various types 
of intermolecular interactions acting between sides, planes and ends of a pair of 
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molecules^. In view of the key role of molecular interactions in mesogenic 
compounds, semi-emperical calculations have been emphasized by several workers to 
explain liquid crystallinity'*. Perrin and Berges have employed PCILO, INDO, CNDO 
etc. methods to analyse (i) the internal rotations, (ii) possibilities of motion in 
aromatic core as well as in the terminal chains, and (iii) the influence of conjunction 
between the oxygen lone pairs and the benzene ring on the internal rotations in several 
mesomorphic compounds^’ Tokita et af used Lennard - Jones potential to evaluate 
intermolecular interactions between a couple of pure nematogens and attempted to 
correlate their results with those of the molecular field theory* . However, it has been 
observed that ‘6-exp’ type of potential functions are found to be more effective in 
explaining the molecular packing instead of ‘Lennard - Jones’ potentials’. 

Using Buckingham ‘6-exp’ type of potential function, intermolecular interaction 
energy studies have been carried out in case of some mesogens in our laboratory*®'*'*. 
It has been observed that the pair potential between such molecules is anisotropic in 
nature which is in general regarded as the prime requirement for the mesophase 
formation in thermotropic liquid crystals* ^ The minimum energy configurations 
obtained through such theoretical calculations correspond to the crystallographic 
structure of molecules in the solid state. Similar observations have been reported by 
others also*®. Although considerable progress has been made in understanding the 
liquid ciystalline phases using computer simulation techniques*’’^®, it seems pertinent 
to study simple models in detail with an aim to elucidate the nature and strength of 
molecular forces, molecular ordering, orientational freedom etc. Since mesogenic 
properties are related to molecular aggregation in a specific manner, probability 
distribution calculations based on interaction energy results are expected to provide 
information about most probable molecular aggregation and tendency to retain 
translational/orientational order at different transition temperatures. 

In the light of the above facts, relative probabilities of various configurations 
obtained through interaction energy calculations have been examined in case of some 
mesogens’*'’*. The present paper embodies the results of probability studies carried 
out in case of 4^-methoxybenzylidene - 4 - acetoxyaniline (MBA). The 
thermodynamic parameters reveal that MBA shows crystal to nematic transition at 
358.2 K and nematic to isotropic melt at 382 K’'*. 

Method of Calculation 

The molecular geometry of MBA has been constructed using crystallographic data 
from literature and standard values of bond lengths and bond angles’*. Net charge and 
corresponding dipole moment components at each of the atomic centres of the 



STATISTICAL STUDY OF MOLECULAR ORDERING IN MBA 


161 


molecule have been computed by CNDO/2 method^. Modified Rayleigh-Schrodinger 
second order perturbation theory alongwith multicentred-multipole expansion 
technique as developed by Claverie and coworkers has been used to calculate 
intermolecular interaction energy between a couple of MBA molecules. According to 
this method, the total interaction energy (^tot) between two molecules is expressed 


E-roj - .Eel+-Epol+£^disp+-^rep (1) 

where Eel, Epol, Emsp and -£rep represent electrostatic, polarization, dispersion and 
repulsion energy components respectively. The formulae for various energy terms are 
given as under. 

Electrostatic energy : 

According to the multicentred-multipole expansion method^*, the electrostatic 
energy term is expressed as : 

.^EL ~ •£qQ+^QMI+-£mIMI (2) 

where £qq, ^qmi and £^mimi etc. represent monopole-monopole, monopole-dipole, 
dipole-dipole and interaction energy terms consisting of multipoles of higher orders 
respectively. However, consideration up to the first three terms has been found to be 
sufficient for most of the molecular interaction problems^*’^’. Hence, only three terms 
have been included in the present calculation. 

Polarization Energy : 

The polarization energy of some molecule (say, s) is obtained as a sum of the 
polarization energies for the various bonds : 

■EpoL=C(.l/2)|Elf'.Af4*' (3) 


(0 

where e'/' 

tflts X 


(4) 
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is the electric field created at the bond u by all surrounding molecules and ^ u is the 
polarizability tensor of this bond. Rxu is the vector joining the atom % in molecule (t) 
to the ‘centre of polarizable charge’ on bond u of molecule (s). 


Dispersion and Repulsion Energy: 

Dispersion and repulsion terms are calculated together using Kitaigorodskii type 
of formula as given below 

i^DisP + £rep = 1 1'^ I (5) 

where 


E{X,\) = Kx K, (-A/z‘ + Be'^") 


( 6 ) 


and 


z=RxdRl=yf(mm (7) 

where R^ and R^ are the van der Waals radii of atoms k and v respectively. The 
parameters A, B and y do not depend on the atomic species: this necessary dependence 
is brought about by R°, and the factors Kx and which allow the energy minimum to 
have different values according to the atomic species involved^^. The values of these 
various parameters and of the van der Waals radii have been given by Caillet and 
Claverie^^’^'* . 

The energy minimization has been carried out for both stacking and in-plane (side 
to side and end to end) interactions separately. One of the interacting molecules is 
kept fixed throughout the process while both lateral and angular variations are 
introduced in the other in all respects relative to the fixed one. The first molecule has 
been assumed to be in the X-Y plane with X-axis lying along the long molecular axis 
while the origin is chosen approximately at the centre of mass of the molecule. The 
second molecule has been translated initially along the Z-axis (perpendicular to the 
molecular plane) and subsequently along X-and T-axes. Variation of interaction 
energy with respect to rotation about Z-axis has been examined in the range of ± 60°. 
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Accuracies up to 0.1 A° in sliding (translation) and 1° in rotation have been achieved. 
The details of the mathematical formalism and optimization process may be found in 
literature 

The calculation of probabilities have been carried out using Maxwell-Boltzmann 
formula as given below: 

Pi= exp (-Ps,)/i;exp(-PS() (8) 


where Pi stands for the probability corresponding to the configuration /; P= MkT, 
where k and T are Boltzmann constant and absolute temperature respectively and S/ 
represents the energy of configuration i relative to the minimum energy value in a 
particular set for which the probability distribution is being computed. Molecular 
charge distribution, stacking and in-plane (side to side) interaction energy results are 
published elsewhere". In this paper, results obtained from probability distribution 
calculations have been reported. 

Results and Discussion 

The molecular geometry of MBA has been shown in Fig. 1. The variation of 
probability with respect to translation of one of the stacked molecules along the 
Z-axis corresponding to four sets of axial rotation viz. X{0^Y{0\ A(0°)y(180°), 
A(180'’)y(0®) and Ar(180®)7(180°) has been depicted m Fig. 2. It is observed that 
configuration A!(0°) Y( 1 80°) and X{ 1 80°) F( 1 80°) show a sharp peak at the same inter- 
planar separation, 4.0A° with nearly equal probability i.e. 91%. Below 4.0A°, the 
probability of occurrence of these configurations reduces to zero. The configuration 
X(0°)Y(0°) and X(180°)Y(0°) show their maxima at 4.5A° with probability values 
nearly 70% and 60% respectively. An analysis of the relative probability of being at 
maximum point corresponding to four temperatures viz. 300K (room temperature), 
358.2K (solid to nematic), 382.0K(nematic to isotropic) and 500K (above room 
temperature), is shown in Table 1. As evident from this table, the configuration 
Y(180°)y(180°) has almost 66% probability of occurrence at nematic to isotropic 
transition temperature while configurations Y(0°) 7(0°) and X(0°) 7(180°) have nearly 
probabilities 7% and 27% respectively. The remaining one configuration namely, 
X( 180°) 7(0°) possesses negligible probability. This indicates that in a stacked pair of 
MBA molecules, the probability of the configuration in which one of the molecules is 
rotated by 180° both about AT and 7-axes relative to the other, is maximum. 
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Table 1- Comparative probability of the minimum energy stacked complexes of MBA molecules with 
respect to translation along Z-axis corresponding to various rotational sets 


Configuration 

Separation 

(A“) 

Energy 

(kcal/mole) 


Probability (%) 



7= 300K 

r= 358.2K 

7= 382K 

r=500K 

x(0)r(0) 

4.0 

-6.90 

3.98 

5.95 

6.74 

10.30 

A'( 180)7(0) 

4.5 

4.71 

0.10 

0.27 

0.37 

1.13 

X (0)7(1 80) 

4,0 

-7.97 

23.64 

26.42 

27.28 

29.97 

X( 180)7(1 80) 

4.0 

-8.63 

72.27 

67.34 

65.59 

58.50 



Fig. 1 -Molecular geometry of 4^-methoxybenzylidene-4-acetoxyamline (MBA). 
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translation ( A°) — »- 



Fig. 2-Variation of probability with respect to translation of one of the stacked molecules along 
Z-axis corresponding to four distinct rotational sets namely X(0°)Y(0°), 

Y(0“)r(180“), XIlSO^IFfO®) and Y(180“)>'(180“) at 300K. 

Fig. 3 represents the variation of probability with respect to translation of one of 
the stacked molecules along the long molecular axis (X-axis) for configuration 
X(180‘’) 7(180°) in the range of ± 10A° corresponding to three selective rotations 
about Z-axis by 0°, 90° and 270° .The inter-planar separation between the molecules is 
4.0A°. The energy obtained in case of Z(1 80°) is always repulsive. Hence, probability 
is found to be zero at each point in this case. It may be observed that configuration 
Z(0°) shows a relatively sharp peak with a maximum probability of 52% at room 
temperature while Z(90°) and Z(270°) curves exhibit lesser values. Comparative 
probability of minimum energy stacked complexes with respect to translation along X- 
axis, is shown in Table 2 which shows that configuration Z(0°) has maximum 
probability of occurrence, 97% at nematic to isotropic transition temperature. 

A graphical representation of the probability distribution with respect to 
translation of one of the stacked molecules along an axis perpendicular to the long 
molecular axis and lying in the molecular plane (F-axis) in the range of ±2.0A° has 
been shown in Fig. 4 which indicates an overlapped structure as most probable one. 
Since the peak is not very sharp, translation along this axis in a very small range 
(<0.5A°) is probable at increased thermal agitations. Fig. 5 exhibits the variation of 
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probability with respect to rotation about Z-axis in the range of ±60° at 300K. It is 
evident from this figure that maximum probability (approximately 61%) lies at 
perfectly aligned structure. 

TRANSLATION (A°) ► 


-10.0 -5.0 0 5.0 10.0 



Fig. 3-Variation of probability with respect to translation along the long molecular axis (X-axis) 
in the range of ± lO.OA^ corresponding to three selective rotations about Z-axis 
viz. Z{Q% Z{9(f) and Z(270°). The inter-planar separation between 
the molecules is 4.0A^ at 300K. 

Table 2- Comparative probability of the minimum energy stacked complexes of MBA molecules with 
respect to translation along the long molecular axis corresponding to four selective rotations 
about Z-axis 


Rotation 

about Z-axis 

Energy 
( kcal/mole) 


Probability (%) 


r=300K 

r=358.2K 

T=382K 

J'=500K 

0° 

-8.63 

99.07 

97.81 

97.13 

92.64 

90® 

-5.58 

0.58 

1.33 

1.73 

4.28 

o 

O 

OO 

- 

- 


- 

- 

270® 

-5.24 

0.33 

0.84 

1.12 

3.07 
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Fig 4-A graphical representation of the ibility distribution with respect to translation 
ot one of the stacked molecules alo .g i axis perpendicular to long molecular 
axis (F-axis) in the : .i _■ of + 2.0A° at 300K. 



Fig. 5-Variation of probability with respect to rotation auuut Z-axis in the range of ± 60^. 
One of the stacked molecules has Z(180®) F(180°) Z(0^) configuration with respect 
to other and is separated at an inter-planar distance of 4.0A^ at 300K. 
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The variation of probability at 300K with respect to in-plane (side to side) 
translation of one of the interacting molecules along Z-axis has been plotted in Fig. 6 . 
The intermolecular separation along T-axis is 7.0A** and one of the molecules 
possesses X(0®)T(180®)Z(0®) configuration. As evident from this figure, maximum 
probability (nearly 30%) occurs when one of the interacting molecules is displaced by 
6.0A° along ^-axis with a simultaneous rotation of 180° about f-axis. Fig. 7 shows a 
plot of probability distribution for translation of one of the interacting molecules along 
f-axis during in-plane interactions at an interval of 0.2A°. This figure shows a sharp 
peak at an intermolecular separation of 7.0A°. Further, it may be observed here that 
like the case of stacking interactions, distinct maximum exists during in-plane 
interactions in both the cases (Fig. 6 & Fig. 7). 


TRANSLATION ( A° ) ► 



Fig. 6-A plot of probability distribution for in-plane (side to side) translation of one of the 
interacting molecules along Z-axis. The intermolecular separation (along F-axis) is 7.0 A° 
and one of the molecules possesses A'(0'’)F(180®) Z(0°) configuration at 300K. 

The variation of probability with respect to change of inter-group separation 
during terminal interactions for the energetically most stable interacting terminal 

groups viz. (-OCOCH3 H3COCC)-) configuration has been shown in Fig. 8 . It is 

observed from this figure that probability of any stable end to end configuration is 
very small (< 7 %) even at room temperature. The optimum end to end separation in 

this case (-OCOCH3 H3COCO-), is 3.03A° whereas the centre of mass 

positions of interacting MBA molecules are separated by 18.8A°. Further other 

possible interacting groups i.e. (-OCOCH 3 H 3 CO-) and (-OCH 3 H 3 CO-) 

show negligible probability of occurrence. 
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Fig. 7-A plot of probability distribution for translation of one of the interacting molecules along 
T*axis during in-plane (side to side) interactions at an interval of 0.2 A®. The initial 
value of intermolecular separation is 7.0 A® at 300K. 



TRANSLATION (A) 

Fig. 8~A plot of probability distribution for translation of one of the molecules along X-axis during 

in-plane (end to end) interactions corresponding to -OCOCH3 H3COCO- interacting 

groups at an interval of 0.1 A®. The initial value of intermolecular separation is 
19.0A° while optimum end to end distance is 3.03 A*^. The calculation 
has been carried out at 300K. 
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Table 3- Translational and rotational rigidities corresponding to most probable configuration during 
stacking, in - plane (side to side) and terminal (end to end) interactions in case of a pair of 
MBA molecules at various temperatures. 


Temperature 

(K) 



Translational Rigidity 



Rotational 
Rigidity about 
Z-axis (±10°) 



Stacking 


In-plane 
(side to side) 

Terminal 



X-axis 

T-axis 

Z-axis 

X-axis 

T-axis 

X-axis 



(+2.0A‘’)* 

(±0.4A'') 

(+0.5A®)*» 

(± 2.0A“) 

(± 0.2A°) 

(± 0.2A®) 


300.0 

1.26 

0.74 

18.71 

1.38 

0.96 

0.59 

1.82 

358.2 

1.21 

0.69 

11.63 

1.24 

0.88 

0.57 

1.48 

382.0 

1.20 

0.68 

9.97 

1.22 

0.85 

0.56 

1.38 

500.0 

1.15 

0.64 

5.80 

1.04 

0.76 

0.55 

1.09 


*At -2.0A® interaction energy comes out to be repulsive due to violation of van der WaaTs contact 
distance and hence probability of the configuration vanishes as shown in Fig.3. 

** At 3.5A° interaction energy comes out to be repulsive due to violation of van der Waal’s contact 
distance and hence probability of the configuration vanishes as shown in Fig.2 

Translational and rotational rigidities corresponding to most probable 
configurations during stacking, in-plane and terminal interactions at different 
temperatures, have been shown in Table 3. It should be mentioned here that rotational 
rigidity has been defined as the probability ratio of being at maximum to that 
corresponding to ±10° of rotations. Translational rigidity has been defined as the 
probability ratio of being at maximum to that corresponding to ±2.0A°, ±0.4A° and 
±0.5A° of translations of one of the stacked molecules along X-, Y- and Z-axes 
respectively during stacking interactions. During in-plane (side to side) interactions 
translational rigidity has been calculated in a range of ±2.0A°forZ-sliding and ±0.2A° 
for 7-sliding whereas the same has been considered in the range of ±0.2A° for 
translation along A'-axis during in-plane (end to end) interactions. As evident from 



STATISTICAL STUDY OF MOLECULAR ORDERING IN MBA 


171 


Table 3, both rotational and translational rigidities decrease with increase in temp- 
erature in each case. The probability of any stable configuration during terminal 
interactions is poor. Similarly during stacking and in-plane (side to side) interactions, 
there is a lesser probability of existence of stable configurations corresponding to the 
translation along an axis perpendicular to the long molecular axis (7-axis). Further, 
during translation of one of the stacked molecules along X- and Z-axes, interaction 
energy comes out to be repulsive at the intermolecular separation of -2.0 A® and 
3. 5 A® respectively due to the violation of van der Waal contact distances in both the 
cases (Fig. 3 & Fig. 2). However, these results favour nematic behaviour of the 
system. Similar results have been obtained in case of para-azoxyanisole^’ . 

Conclusion 

It seems that MBA exhibits a strong preference for aligned structure at transition 
temperature. In a stacked molecular pair both orientational flexibility and transitional 
freedom are found to be small corresponding to minimum energy configuration. Other 
configurations show greater translational flexibility alongwith their intrinsic 
preference for aligned structure which accounts for the nematic behaviour of the 
molecule 
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